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Abstract 


A  general  predictive  methodology  for  determining  residual  strength  in 
impact  damaged  composite  laminates  has  been  developed  and  incorporated  into 
a  computer  code  designated  RESTRAN  (REsidual  STRength  ANalysis). 
RESTRAN  is  a  finite  element  based  design  tool  that  can  analyze  composite 
structures  with  arbitrary  three-dimensional  (3-D)  geometry,  loading  and  support 
conditions,  material  properties,  and  initial  material  and  delamination  damage. 
Material  failure  modes  are  predicted  using  a  robust  suite  of  failure  criteria  and 
damage  laws.  Structural  failure  due  to  sequential  sublaminate  buckling  of 
delaminated  layers  is  also  accounted  for.  A  progressive  failure  analysis  is 
performed  until  ultimate  structural  failure  is  predicted,  thereby  yielding  an 
estimate  of  the  residual  strength.  This  report  contains  a  user's  manual  for  the 
RESTRAN  program  with  complete  descriptions  of  input  statements  and  program 
output.  Several  examples  are  shown  to  illustrate  the  use  of  the  RESTRAN 
computer  code. 
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1  Introduction 


The  development  of  the  RESTRAN  (REsidual  5ri2eng:th  Analysis)  computer  program  represents 
a  significant  contribution  as  a  design  tool  to  predict  the  reduction  in  structural  strength  due  to 
impact  damage.  The  numerical  finite  element  basis  provides  the  greatest  versatility  for  modeling 
structures  of  arbitrary  geometry,  support,  and  load  conditions  while  accounting  for  both  material 
and  structural  failure  modes.  Although  designed  specifically  for  the  analysis  of  advanced  fibrous 
composites,  the  support  of  a  layered  media  permits  additional  material  systems  to  be  modeled  such 
as  sandwich-type  constructions,  piecewise  linear  approximations  of  homogeneous  materials  with 
varying  properties  along  a  thickness  direction,  and  isotropic  materials  modeled  as  a  single  layer. 
For  composites,  material  failure  may  be  differentiated  into  fiber  and  matrix  modes,  while  structural 
failure  involves  buckling  of  local  sublaminates  due  to  the  presence  of  delamination  damage.  In 
addition  to  built-in  features,  a  suite  of  user-defined  subroutine  interfaces  are  available  for  problem- 
specific  tailoring  of  the  analysis.  Optional  graphics  may  be  requested  and  output  as  input  files 
for  display  using  MATHEMATICA  [1],  TECPLOT  [2],  or  in  a  user-defined  format.  The  analysis 
capabilities  of  RESTRAN  are  thoroughly  discussed  in  the  Theoretical  Manual  [3].  The  current 
report  details  the  use  of  the  RESTRAN  code.  Complete  descriptions  are  given  of  the  program 
input  statements  to  define  the  structural  geometry,  support  and  load  conditions,  material  behavior, 
and  solution  procedures.  Finally,  a  series  of  demonstration  problems  showing  program  input  and 
output  are  provided  to  illustrate  the  use  and  capabilities  of  RESTRAN. 

2  RESTRAN  Input  Statements 

The  following  sections  detail  the  syntax  and  basic  format  of  input  statements  available  to  describe 
the  model  geometry  and  invoke  program  options.  Some  degree  of  flexibility  has  been  implemented 
in  the  parsing  of  these  statements  such  that  either  upper  or  lower  case  letters  can  be  used,  and 
arbitrary  space  or  comma  delimitators  can  be  used  in  listing  statement  options  and  data  items. 
The  basic  format  of  all  input  statements  is: 

♦STATEMENT,  OPTIONi  =  X,  OPTION2  =  Y,  ...  ,  OPTION^  =  Z 
DATAi,  DATA2,  DATA3,  ...  ,  DATA„ 

One  requirement,  however,  is  that  the  *STATEMENT  entry  must  begin  in  the  first  field  in 
the  record  or  line  of  the  input  file.  The  input  statements  may  generally  occur  in  any  order  and  are 
parsed  sequentially  as  contained  in  the  input  file.  Statements  may  be  repeated  to  supersede  prior 
instructions  or  to  append  additional  quantities  to  the  input  database  such  as  nodes  or  elements. 

In  developing  syntax  for  RESTRAN  input,  many  of  the  primary  statements  such  as  node,  element, 
boundary,  and  mesh  generation  axe  based  on  the  syntax  used  in  the  commercial  finite  element 
analysis  programs  ABAQUS  [4].  Other  statements  have  a  parallel  in  NASTRAN  [5].  This  was 
implemented  so  that  preprocessor  programs  such  as  PATRAN  [6],  which  facilitate  the  preparation 
of  input  for  these  commercial  codes,  could  be  used  to  create  the  majority  of  the  input  data  for 
RESTRAN.  The  addition  of  a  small  set  of  RESTRAN-specific  statements  completes  the  input 
description  to  perform  the  various  analysis  options.  The  input  statements  supported  are  listed  in 
Table  1. 
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Table  1:  RESTRAN  input  statements. 


STATEMENT 

CATEGORY 

PAGE 

*  PREPASS 

EXECUTIVE 

3 

*  SOLUTION 

3 

*  MEMORY  ALLOCATION 

SYSTEM  RESOURCES 

5 

*  NODE 

5 

*  NGEN 

6 

*  NFILL 

6 

*  ELEMENT 

7 

*  ELGEN 

8 

*  NSET 

8 

*  NSET2 

8 

*  ELSET 

9 

*  ELSET2 

9 

*  BOUNDARY 

9 

*  BOUNDARY2 

9 

*  CLOAD 

10 

*  PRESSURE 

MODEL  DESCRIPTION 

10 

*  MODEL  GENERATION 

11 

*  LAYER 

17 

*  MATERIAL 

17 

*  NONLINEAR  MATERIAL 

18 

*  FAILURE  CRITERION 

21 

*  DAMAGE  LAW 

24 

*  INITIAL  DAMAGE 

25 

+  DELAMINATION 

27 

*  POSTBUCKLED  MATERIAL 

27 

*  ORIENTATION 

28 

*  DEFORMED  GEOMETRY 

29 

*  EQUIVALENCE 

29 

*  INCLUDE 

37 

*  EXCLUDE  ELEMENT 

30 

*  ENDDATA 

30 

*  PARAMETER  USRBKL 

31 

*  PARAMETER  STATUS 

EXECUTION  MODIFICATION 

33 

*  PARAMETER  LAPFAIL 

33 

*  PARAMETER  DIRECTORY 

34 

*  NODE  PRINT 

34 

*  ELEMENT  PRINT 

34 

*  PLY  FAILURE  PRINT 

OUTPUT  CONTROL 

34 

*  GRAPHICS 

35 

*  ECHO 

37 

*  HEADING 

37 
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2.1  Executive  Control 


Two  statements  are  available  to  dictate  executive  control  over  the  execution  of  RESTRAN.  These 
statements  are  used  to  invoke  preliminary  model  checking  or  initiate  various  analysis  procedures. 


2.1.1  Prelimineiry  Model  Evaluation 

To  perform  a  detailed  check  on  the  inputted  model,  the  following  executive  control  statement  can 
be  included  anywhere  in  the  input  stream. 

*PREPASS 

This  statement  causes  RESTRAN  to  process  the  input  file,  perform  error  checking,  assess  memory 
requirements,  and  generate  optional  graphics  of  the  model.  This  option  overrides  any  other  speci¬ 
fied  solution  procedure  and  is  used  to  assess  the  validity  of  the  inputted  model  prior  to  performing 
any  further  analysis. 


2.1.2  Analysis  Procedures 

Several  solution  procedures  are  available  in  RESTRAN.  The  executive  control  statement  to  request 
specific  analysis  procedures  is  given  by 

♦SOLUTION,  METH  =  PROG 
LIMi,  LIM2,  ACCEL 
PROC,  ITERB,  CVAL,  MTOL,  RTOL 

This  statement  selects  the  analysis  procedure  to  be  invoked  by  RESTRAN  by  assigning  the  manda¬ 
tory  parameter  METH.  Two  linear  analyses  and  three  residual  strength  procedures  may  be  se¬ 
lected.  For  linear  static  analysis,  the  method  parameter  is  set  as 

METH  =  LSA 

This  directs  RESTRAN  to  perform  a  linear  static  analysis  to  determine  the  basic  elastic  response 
of  the  structure.  This  option  may  be  used  to  check  displacements,  stresses  and  strains,  and  an 
optional  graphical  depiction  of  deformation  state,  to  assess  applied  loads  and  support  conditions 
imposed  on  the  model.  In  addition,  the  inclusion  of  the  *PLY  FAILURE  PRINT  option  will  out¬ 
put  first-ply  failure  statistics.  The  data  list  following  the  *SOLUTION  statement  can  be  omitted. 

METH  =  LBA 

This  option  directs  RESTRAN  to  perform  a  linear  buckling  analysis.  This  procedure  may  be 
used  to  analyze  stability  in  a  multiply-delaminated  composite  with  all  delaminations  considered 
simultaneously.  The  predominant  buckling  mode  can  then  be  determined  such  that  only  a  sin¬ 
gle,  critical  delamination  may  be  input  for  a  subsequent  residual  strength  analysis.  This  would 
greatly  speed  execution  because  the  residual  strength  algorithm  will  consider  each  delamination 
in  sequence  to  automatically  choose  the  critical  delamination  involved  in  sublaminate  buckling. 
Optional  graphical  output  can  be  selected  to  view  the  buckled  mode  shape.  As  with  the  LSA 
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solution  procedure,  the  data  list  following  the  *SOLUTION  statement  can  be  omitted. 

To  invoke  an  incremental/iterative  residual  strength  analysis,  the  following  options  are  available. 


METH  =  < 


MTL 

BKL 

CMB 


(1) 


These  options  specify  the  range  of  failure  modes  to  be  accounted  in  the  analysis.  The  METH  = 
MTL  parameter  is  used  to  specify  that  only  material  failure  modes  are  to  be  processed.  METH 
=  BKL  causes  only  sublaminate  buckling  failures  to  be  analyzed,  and  METH  =  CMB  directs 
the  solution  algorithm  to  consider  both  material  and  buckling  failure  modes. 


The  control  parameters  set  basic  iteration  bounds  in  the  residual  strength  analysis.  LIMi  is 
the  maximum  number  of  sequential  load  increments  to  be  calculated  in  the  prediction  of  residual 
strength.  These  loads  are  automatically  determined  as  scalar  factors  applied  to  the  initial  set  of 
external  loads,  and  this  limit  may  be  set  arbitrarily  high.  In  certain  cases  where  user-defined  sub¬ 
routines  are  used  to  degrade  material  properties  undergoing  failure,  arbitrary  degradation  laws  may 
be  input  which  could,  in  some  situations,  cause  internal  tests  performed  by  RESTRAN  to  avoid  de¬ 
termining  catastrophic  failure.  An  upper  bound  on  the  number  of  total  analysis  cycles  is  therefore 
needed  to  stop  runaway  execution.  LIM2  is  thus  used  to  set  the  maximum  number  of  iterations  to 
converge  the  damage  state  at  the  cmrent  load  level.  For  each  iteration,  failures  are  accounted  and 
internal  loads  redistributed,  which  may  cause  additional  failures  to  occur.  Although  most  failure 
events  tend  to  converge  rapidly  to  an  equilibrium  damage  state  or  cascade  to  ultimate  structural 
collapse,  this  limit  sets  bounds  on  possible  slow  failure  cascades  which  may  warrant  changes  in 
other  input  parameters.  The  optional  ACCEL  parameter  is  used  to  accelerate  the  prediction  of 
failure  by  encompassing  larger  groups  of  local  failures  per  analysis  cycle.  RESTRAN  automatically 
determines  a  scale  factor  applied  to  the  vector  of  initial  input  loads  to  identically  satisfy  selected 
failure  criteria  for  a  critical  ply  in  a  particular  element.  This  set  of  loads  is  then  increased  by 
a  default  value  of  1%  to  cause  other  plies  to  fail  which  may  have  calculated  failure  indices  near 
unity.  The  value  of  the  ACCEL  parameter  overrides  the  default  value  such  that,  for  example,  an 
input  value  of  1.1  will  cause  all  plies  within  10%  of  the  critical  ply  load  to  fail  in  the  current  cycle 
and  thus  accelerate  the  prediction  of  local  failures.  Care  must  be  given  in  setting  the  value  of  this 
parameter,  as  it  directly  modifies  the  applied  loads  and  therefore  the  predicted  value  of  the  residual 
strength.  A  value  less  than  unity  would  preclude  material  failure  and  is  therefore  automatically 
reset  to  the  default  value.  A  large  value  would  cause  an  overprediction  of  the  residual  strength  and 
would  have  to  be  manually  factored  out  of  the  final  predicted  strength  measure  by  the  user. 

For  residual  strength  analysis  involving  sublaminate  buckling  failure,  PROC  may  be  set  equal 
to  ’SIMULTANEOUS’  to  analyze  all  delaminations  as  a  single  group  or  ’ITERATIVE’  to  analyze 
each  delamination  separately  and  to  iterate  ITERB  times  to  obtain  a  solution  to  the  constrained 
buckling  contact  problem  of  an  embedded  delamination.  To  approximately  solve  the  contact  con¬ 
straints  involved  with  sublaminate  buckling,  iterations  are  performed  to  minimize  the  degree  of 
interpenetration  of  a  buckling  mode  to  surrounding  layers.  This  is  performed  by  repeated  solution 
of  a  linear  buckling  analysis  of  a  particular  delamination.  The  complete  satisfaction  of  constraints 
is  impractical  and  unnecessary  to  obtain  reasonable  local  buckling  loads.  The  parameters  CVAL, 
MTOL,  and  RTOL  specify  specific  treatment  of  delamination  buckling  and  range  from  0  to  1  due 
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to  an  internal  normalization  of  the  buckling  mode  vector.  To  limit  the  number  of  iterations  and 
to  specify  a  criteria  for  accepting  a  particular  buckling  mode,  the  fraction  of  nodes  exhibiting  a 
physical  opening  mode  is  required.  The  default  value  of  CVAL  is  0.9,  indicating  that  a  mode  is 
accepted  if  only  10  %  of  nodes  in  the  delamination  plane  are  involved  in  nonphysical  interpenetrat¬ 
ing  displacements.  The  entry  MTOL  indicates  the  ratio  of  maximum  local  modal  displacements 
across  the  delamination  plane  with  the  maximum  overall  buckling  mode  displacement.  This  ratio 
is  used  to  predict  whether  a  mixed-mode  instability  is  occurring  wherein  the  delamination  being 
analyzed  may  be  undergoing  an  opening  mode;  however,  the  overall  structure  is  exhibiting  a  greater 
deflection  indicating  a  global  mode.  The  default  value  for  MTOL  is  0.5.  Setting  this  tolerance  to 
zero  will  preclude  a  mixed  global-local  buckling  mode  from  being  predicted.  Finally,  the  RTOL 
entry  indicates  the  threshold  below  which  the  relative  motion  of  opposing  delamination  surfaces  are 
to  be  considered  zero,  and  coincident  nodes  are  to  be  condensed  out  of  the  current  delamination 
plane.  Unless  specified,  the  default  value  for  RTOL  is  0.0. 


2.2  System  Resources 

A  single  large  array  is  used  within  RESTRAN  to  perform  numerical  operations  on  global-sized  ma¬ 
trices.  This  in-core  array  is  used  to  provide  working  storage  to  perform  fast  matrix  operations  with 
the  bulk  of  generated  internal  data  stored  in  external  scratch  files.  In  the  normal  default  execution 
mode,  RESTRAN  uses  a  half  bandwidth  storage  scheme  internally.  In  order  to  run  large  jobs,  an 
out-of-core  solution  algorithm  may  be  invoked.  This  causes  additional  processing  and  will  slow  down 
the  execution  speed.  The  internal  algorithm  for  memory  usage  is  specified  using  the  following  entry: 

♦MEMORY  ALLOCATION 
STYPE 

where  STYPE  =  ’BANDWIDTH’  for  half  bandwidth  storage  and  ’OUT-OF-CORE’  for  fully 
external  storage.  A  storage  format  of  ’FULL’  is  recognized  which  causes  all  operations  to  be 
performed  using  full  matrix  storage.  The  only  practical  use  of  this  feature  is  to  speed  execution 
of  small  problems  in  which  the  connectivity  unavoidably  yields  a  global  stiffness  matrix  that  is 
densely  populated  as  opposed  to  sparse.  If  this  statement  is  omitted,  the  default  storage  mode  is 
BANDWIDTH. 


2.3  Model  Description 

The  following  subsections  detail  the  RESTRAN  input  statements  used  to  define  a  structural  model. 


2.3.1  Node  Definition  and  Generation 
Nodes  are  input  using  the  following  statement: 

♦NODE,  SETID  =  NSD 
NODE  ID,  X,  Y,  Z 

where  the  node  number  and  position  in  global  x,  y,  z  coordinates  are  repeated  for  each  node 


5 


in  the  model.  The  optional  parameter  SETID  is  used  to  include  the  specified  nodes  into  a  new  or 
existing  set.  The  node  set  ID  must  be  numerical;  general  character  strings  are  not  supported. 

To  automatically  generate  nodes,  two  utility  routines  may  be  invoked  to  generate  a  linear  and 
planar  field  of  nodes.  These  are  specified  by  the  commands  *NGEN  and  *NFILL. 

To  generate  a  linear  array  of  nodes,  the  following  statement  is  used 

*NGEN,  SETID  =  NGI 
NODEi,  NODE2,  INCR 

As  shown  in  Figure  1,  *NGEN  generates  a  line  of  nodes  between  two  previously  defined  nodes, 
NODEi  and  NODE2.  The  increment  value,  INCR,  specifies  the  increment  in  the  node  numbers 
of  equally  spaced  nodes  to  be  generated  between  the  two  existing  nodes.  This  increment  thereby 
specifies  the  number  of  nodes  generated,  which  is  given  by 

NUMGEN  =  {\NODE2  -  NODEi \/INCR)  -  1  (2) 

The  generated  nodes  may  optionally  be  assigned  to  a  node  set  by  including  the  SETID  parameter. 
This  option  must  be  specified  if  the  *NFILL  command  is  to  be  used  with  these  generated  nodes. 


A  two-dimensional  plane  of  nodes  may  be  generated  using  the  following  statement: 

*NFILL,  SETID  =  NFD 
NSETi,  NSET2,  NUMGEN,  INCR 

As  shown  in  Figure  2,  two  lines  of  nodes  with  specified  set  identification  numbers  NSETi,  NSET2 
are  used  to  generate  a  plane  of  nodes.  In  generating  the  two-dimensional  field  of  nodes,  NUM¬ 
GEN  indicates  the  number  of  equally  spaced  nodes  to  be  generated  between  bounding  nodes  in 
each  set,  and  INCR  indicates  the  increment  in  node  numbering  between  each  set.  The  node  set 
parameter  is  optional. 
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2.3.2  Element  Definition  and  Generation 
Element  definitions  are  input  using  the  following  statement: 

♦ELEMENT,  LAYUP  =  L,  ORIENTATION  =  M,  ELSET  =  EID 
EID,  Ni,  N2,  N3,  N4,  N5,  Ne,  N7,  Ng 

The  data  list  gives  the  element  number  and  the  node  numbers  of  the  nodes  defining  the  ele¬ 
ment  geometry  as  shown  in  Figmre  3.  The  mandatory  parameter,  LAYUP  =  L,  associates  the 
elements  with  a  particular  input  sublaminate.  The  optional  orientation  parameter  designates  the 
identification  number  of  an  inputted  local  coordinate  system  to  be  assigned  to  these  elements.  The 
ELSET  parameter  is  also  optional  to  assign  each  inputted  element  to  an  element  set.  The  data 
list  is  repeated  for  all  elements  associated  with  the  specified  layup  and  local  coordinate  system. 
The  *ELEMENT  data  block  is  repeated  for  each  group  of  elements  with  a  unique  assigned  layup 
and  spatial  orientation. 


Figure  3.  Hexahedral  element  geometry. 
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Automated  generation  of  elements  is  supported  by  the  following  input  statement: 

*ELGEN,  ELSET  =  EID 

ELID,  NGNx,  NODa,  NIDa,  NGNy,  NODy,  NIDy,NGN2,  NODz,  NIDz 

ELID  indicates  the  starting  element  from  which  other  elements  are  to  be  generated.  The  pa¬ 
rameters  NGNA.y.orZ  indicate  the  number  of  elements  to  be  generated  along  the  global  X,  Y,  or 
Z  directions,  NODA,y,orZ  indicates  the  increment  in  node  numbers,  and  NID a, y,or2  indicates  the 
increment  in  element  numbers  along  the  X,  Y,  or  Z  directions,  respectively.  The  optional  parameter 
ELSET  may  be  used  to  add  the  generated  elements  to  a  new  or  existing  element  set.  Element 
generation  is  depicted  in  Figure  4. 


Figure  4.  Generation  of  a  one-,  two-  or  three-dimensional  field  of  elements. 

2.3.3  Node  and  Element  Sets 

Node  and  element  sets  are  used  to  define  delamination  planes,  initial  damage,  and  in  print  request 
statements.  Node  sets  are  declared  using  the  following  statements. 

*NSET,  SETID  =  NID 
Ni,  N2,  ...  ,Nfc 

Where  SETID  =  NID  gives  the  set  identification  number,  and  Ni,  N2,  ...  ,Nfc  are  the  set 
of  node  numbers.  The  list  of  nodes  can  contain  any  number  per  line  delimited  by  either  commas 
or  blanks  and  continued  on  an  arbitrary  number  of  additional  lines.  An  alternate  statement  is 
available  for  nodes  defined  within  bounds. 

*NSET2 
NID,  Ni,  N2 
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where  NID  is  the  set  number,  and  Ni  and  N2  give  an  inclusive  range  of  node  numbers  to  be 
formed  into  a  set.  Only  nodes  within  the  given  range  that  are  actually  used  in  the  model  will  be 
stored.  This  line  may  be  repeated  to  define  multiple  sets. 

The  *NSET  and  NSET2  statements  can  be  used  to  exclusively  define  node  sets  or  used  in  conjunc¬ 
tion  with  the  optional  node  set  parameters  on  the  *NODE,  *NGEN,  or  *NFILL  statements. 
Depending  on  the  order  in  which  the  input  records  are  parsed,  any  statement  can  initiate  the 
creation  of  a  new  node  set  while  subsequent  statements  referencing  the  same  node  set  will  be  ac¬ 
cumulated  into  the  existing  set. 

Element  sets  are  defined  using  the  following  command. 

*ELSET,  SETID  =  EID 

El,  E2,  E3,  ...  ,Efc 

where  EID  is  the  assigned  element  set  number,  and  Ei  — >•  E^  is  the  sequence  of  element  numbers 
to  be  included  in  the  set.  The  number  of  element  ID’s  per  line  is  arbitrary,  and  any  number  of 
additional  lines  may  be  used  to  complete  the  set  definition. 

Element  sets  may  additionally  be  defined  using  the  following  command. 

♦ELSET2 
EID  El,  E2 

where  EID  is  the  assigned  element  set  number,  and  Ei  and  E2  give  a  range  of  element  numbers 
to  be  included  in  the  set.  An  arbitrary  number  of  element  sets  may  be  defined  using  this  statement. 

As  with  node  sets,  element  sets  may  be  defined  or  built  up  using  any  combination  of  ’•'ELE¬ 
MENT,  ’•'ELGEN,  ’•'ELSET,  and  ’•'ELSET2  statements. 

2.3.4  Boundary  Constraints  and  Load  Conditions 

Specific  point  constraints  on  displacement  degrees  of  freedom  may  be  applied  to  enforce  fixity  con¬ 
ditions  at  any  node.  Two  input  statements  are  available  to  apply  these  conditions. 

♦BOUNDARY 
NODE,  DOFi,  DOF2 

where  NODE  is  the  node  number.  DOFi  and  DOF2  specify  the  range  of  degrees  of  freedom 
to  be  fixed.  This  line  may  be  repeated  for  all  nodes  with  applied  displacements  constraints.  The 
three  degrees  of  freedom  at  each  node  are  denoted  by  1,2,3  and  correspond  to  the  global  (x,y,z) 
coordinate  or  the  u,  v,  and  w  translational  motions  with  respect  to  the  global  coordinate  system. 
This  is  shown  in  Figure  5.  The  second  statement  allows  a  node  set  to  be  used  to  define  fixity 
conditions. 


’•‘BOUND  ARY2 
NID,  DOFi,  DOF2 
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Figure  5.  Translational  degrees  of  freedom. 

where  NID  is  a  node  set  number,  and  DOFi  and  DOF2  define  the  range  of  degrees  of  free¬ 
dom  to  be  removed  from  each  node  in  the  set. 

External  loads  may  be  applied  as  either  concentrated  forces  at  nodes  or  as  pressure  distributions 
over  an  element  face.  For  concentrated  loads  the  input  entry  is  given  by 

*CLOAD 

NODE,  DOF,  MAG 

where  NODE  is  the  node  number  at  the  point  of  application,  DOF  is  the  direction  of  the  load 
along  the  (1,2,3)  degree  of  freedom  directions  which  correspond  to  the  global  (x,y,z)  coordinate 
system,  and  MAG  is  the  magnitude  of  the  load.  The  direction  of  the  load  is  specified  by  the  sign 
of  the  applied  force.  This  line  may  be  repeated  for  all  nodes  with  applied  concentrated  loads. 

Uniform  distributed  loading  is  applied  using  the  *PRESSURE  statement.  This  internally  in¬ 
tegrates  the  applied  pressures  and  forms  a  consistent  load  vector  which  is  included  with  any  other 
applied  point  loads.  The  syntax  is  given  by 

♦PRESSURE 
ElilD,  Fj,  P o 

where  ELID  is  the  element  number,  F,  is  the  element  face  number  as  shown  in  Figme  6,  and 
Po  is  the  applied  pressure.  The  sign  on  Po  dictates  whether  the  distributed  tractions  are  applied 
opposite  to  the  direction  of  the  outward  normal  to  the  face  yielding  compressive  loads  (-f-Po),  or 
in  the  same  sense  as  the  outer  normal  yielding  tensile  loading  (-Po)- 
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Figure  6.  Convention  for  element  face  numbering. 


2.3.5  Model  Generation 

It  is  expected  that  most  models  input  to  RESTRAN  will  be  generated  using  preprocessing  software 
such  as  PATRAN,  which  supports  templates  for  output  of  a  finite  element  model  in  various  for¬ 
mats.  The  bulk  of  an  input  file  typically  consists  of  node  and  element  definitions;  for  that  reason, 
the  adopted  form  of  RESTRAN  input  syntax  for  node,  element,  force,  and  displacement  constraint 
conditions  is  based  on  a  simplified  form  of  ABAQUS  input  formats.  Thus,  for  example,  output  from 
a  preprocessor  in  ABAQUS  format  can  be  cut  and  pasted  into  a  RESTRAN  input  file  with  only 
minor  modification,  and  only  those  statements  specific  to  RESTRAN  need  be  added  to  complete 
a  job  description.  Alternatively,  the  *INCLUDE  statement  can  be  used  to  temporarily  redirect 
input  to  an  alternate  file  containg  statements  to  be  used  by  RESTRAN. 

However,  in  the  absence  of  preprocessing  software  to  facilitate  the  creation  of  a  finite  element 
model,  an  input  option  is  available  to  create  a  simple  layered  plate  model  containing  a  central  cir¬ 
cular  or  elliptical  region  which  may  be  assigned  reduced  moduli  and/or  used  to  define  delamination 
planes.  The  use  of  this  rudimentary  feature  has  been  assumed  within  a  prepass  mode  after  which 
the  output  may  be  incorporated  into  a  complete  RESTRAN  input  file.  This  statement  is  given  as 
follows 


*MODEL  GENERATION 
NDIV,  NSURF,  A,  B,  NDi,  NDa,  ND3,  ND4 
SECT,  NQi,  NQ2,  NQs,  NQ4,  NI,  NORIENT 
Til,  Ti2,  Ti3 
T21,  T22>  T23 
T31,  T32,  T33 

1)  SCLi,  ZOi 

2)  MATIi,  MATOi,  NSIDi,  NOFj,  NEFj 
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where  the  various  input  parameters  are  identified  next  and  will  be  discussed  subsequently. 


NDIV 

Number  of  inner  divisions 

NSURF 

Number  of  surfaces 

A 

Length  of  rectangle 

B 

= 

Width  of  rectangle 

NDi 

= 

Number  of  divisions  in  outer  regions 

SFCT 

Sizing  parameter 

NQi 

= 

Flag  to  remove  elements  in  quadrants 

NI 

Flag  to  remove  elements  within  ellipse 

NORIENT 

= 

ID  of  orientation  vector 

Tij 

= 

Coordinate  transformation  matrix 

SCLi 

= 

Scale  factor  for  the  surface 

ZOi 

= 

Z-Coordinate  of  surface 

MATIi 

Layup  number  for  material  properties  in  inner  elliptical  region 

MATOi 

= 

Layup  number  for  material  properties  in  outer  elliptical  region 

NSIDi 

Node  set  number  for  nodes  within  elliptical  region 

NOFi 

= 

Starting  node  sequence  number 

NEFi 

Starting  element  sequence  number 

The  input  mixes  data  for  surfaces  and  layers  such  as  starting  node  number,  NOF,  and  start¬ 
ing  element  number,  NEF,  respectively.  The  number  of  surfaces  for  each  node  plane  is  equal  to 
the  number  of  layers  plus  one.  Element  data  is  assumed  to  correspond  to  the  layer  below  the 
current  surface.  Therefore,  input  lines  1  and  2  are  repeated  for  each  surface  starting  at  the  top  of 
the  laminate  and,  for  the  last  siuface,  arbitrary  values  for  layer  quantities  associated  with  elements, 
such  as  material  properties,  must  be  entered  to  avoid  null  values  on  the  input  line  but  are  not  used 
internally. 

The  model  is  formed  in  several  stages  for  each  layer.  The  orientation  is  initially  assumed  parallel 
with  the  global  (x,y,z)  coordinates  with  the  z-coordinate  defining  the  thickness  direction.  First,  an 
inner  square  region  defined  on  — 1  <  <  1  with  side  dimensions  equal  to  two  units  is  meshed 

containing  an  inner  circular  region  with  a  radius  fixed  at  0.65.  The  degree  of  discretization  is  de¬ 
termined  by  specifying  the  number  of  divisions,  NDIV,  as  shown  in  Figure  7.  Next,  a  larger  outer 
square  region  is  created  using  corner  vertices  located  at  (±SFCT,  ±SFCT),  yielding  a  square  with 
side  dimensions  2  x  SFCT.  As  shown  in  Figure  8,  the  corner  regions  are  discretized  according  to 
the  input  values  NDi,  ND2,  ND3,  ND4  which  specify  the  number  of  divisions  to  be  made  in  each 
region.  For  each  surface,  the  scale  factor  SCLj  may  be  used  to  shrink  or  expand  the  size  of  the 
delamination  through  the  thickness.  This  value  is  applied  to  the  node  coordinates  within  the  inner 
region  such  that  {x'^y')  =  SCLi(a:, y)  for  —1  <  (x, y)  <  1.  Acceptable  values  for  this  magnification 
factor  should  be  positive  and  less  than  or  equal  to  unity.  Figure  9  shows  an  enlarged  side  view  of 
a  four-layered  plate  with  a  varying  delamination  radius  created  by  using  SCL  =  1.0,  0.9,  0.8,  0.7 
and  0.6  in  the  five  different  surfaces.  The  relative  size  of  the  inner  square  to  the  outer  square  region 
can  be  varied  through  the  selection  of  SFCT  and  SCL,  and  is  given  by  the  ratio  SCL/SFCT.  A 
final  transformation  is  performed  wherein  the  coordinates  of  the  entire  node  field  is  multiplied  by 
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Figure  7.  Controlling  the  level  of  discretization  using  the  NDIV  parameter 


2SCLF 


H —  N1  — ^  —  N2 


Figure  8.  Discretization  of  outer  corner  regions. 


Surface  Depiction  of  Model 


X 

Figure  9.  Variable  delamination  size  through  thickness. 


the  factor  1/SFCT  and  then  scaled  according  to  the  stretch  factors  A  and  B  to  give  the  desired 
overall  plate  dimensions.  For  A  =  B,  the  inner  region  remains  circular,  while  for  A  ^  B,  different 
elliptical  regions  may  be  generated.  As  shown  in  Figure  10,  various  input  values  for  SFCT,  SCL, 
A,  and  B  can  yield  significantly  different  models. 

The  starting  numbers  for  nodes  in  each  surface  and  elements  in  each  layer  must  be  specified  such 
that  the  last  generated  node  in  a  surface  or  element  number  in  a  layer  does  not  exceed  the  starting 
number  assigned  to  the  node  or  element  in  the  next  surface  or  layer.  For  each  surface,  a  node  set 
number,  NSID,  is  input  to  identify  the  nodes  within  the  in  inner  circular  region.  These  nodes  are 
automatically  stored  as  a  node  set  and  can  be  used  to  define  an  inner  delamination  node  plane. 
For  each  layer,  material  property  identification  numbers  are  given  for  elements  within  and  outside 
the  elliptical  region  in  the  MATI  and  MATO  parameters,  respectively. 

If  symmetry  can  be  used  to  reduce  the  size  of  the  model,  quadrants  may  be  selected  for  removal 
using  the  flags  NQi,  NQ2,  NQ3,  NQ4  in  which  a  value  of  zero  causes  all  elements  and  nodes  in 
the  corresponding  quadrant  to  be  deleted.  The  designation  for  the  quadrants  and  models  exploiting 
various  symmetry  conditions  obtained  using  different  values  for  the  flags  NQi  is  shown  in  Figure 
11.  To  create  a  model  with  an  open  hole,  setting  the  modeling  flag  NI  equal  to  zero  will  cause  all 
elements  within  the  ellipse  to  be  removed. 


Finally,  the  coordinate  transformation  matrix  is  used  to  orient  the  model  in  space.  This  ma¬ 
trix  maps  the  coordinates  as 


"  1 

^11  ^12  ^13 

)■  = 

hi  h2  hs 

< 

y  > 

u- J 

^31  ^32  ^33 

i  -  J 

Figure  12  shows  an  example  with  a  transformation  matrix  given  by 


[T]  = 


0.0 

0.0 

1.0 

0.0 

1.0 

0.0 

1.0 

0.0 

0.0 

(4) 


together  with  the  modeling  option  flag  NI  =  0  which  creates  a  layered  model  of  a  vertical  plate 
with  an  open  elliptical  hole. 

This  routine  has  been  included  to  create  a  simple  model  and  was  designed  to  be  used  in  a 
*PREPASS  execution  mode  with  graphics  outputted  to  show  the  generated  model.  The  gen¬ 
erated  elements  and  nodes  are  output  to  the  file  rcstTun.dlg,  which  is  formatted  to  list  elements  and 
nodes  together  with  nodes  identified  along  the  various  boundaries.  This  format  is  shown  in  Figure 
13.  Upon  completion,  the  node  and  element  definitions  that  are  generated  may  be  subsequently 
included  into  a  new  input  file  with  all  loads  and  boundary  conditions  specified  in  order  to  perform 
atr  analysis. 
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Figure  11.  Selection  of  quadrants  to  model  symmetry  conditions. 


Figure  12.  Rotated  plate  in  space  with  an  open  elliptical  hole. 
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♦NODE 


1  O.OOOOOOOE+00  O.OOOOOOOE+00  0.3000000E+01 

2  O.OOOOOOOE+00  0 . lOOOOOOE+01  0 . 3000000E+01 

3  0.2000000E+01  0 . lOOOOOOE+01  0 . 3000000E+01 

4  0 .2000000E+01  O.OOOOOOOE+00  0 . 3000000E+01 

5  0.2000000E+01  -0 . lOOOOOOE+01  0 . 3000000E+01 

6  O.OOOOOOOE+00  -0 . lOOOOOOE+01  0 . 3000000E+01 

10  O.OOOOOOOE+00  0.2666667E+01  0 . 3000000E+01 


*NSET,  i 

WSID  = 

55 

1002 

1003 

1004 

1005 

1006 

1007 

1008 

1009 

1010 

1011 

1012 

1013 

1014 

1015 

1016 

1017 

1018 

1019 

1020 

1021 

1022 

1023 

1024 

*NSET,  1 

WSID  = 

60 

2002 

2003 

2004 

2005 

2006 

2007 

2008 

2009 

2010 

2011 

2012 

2013 

2014 

2015 

2016 

2017 

2018 

2019 

2020 

2021 

2022 

2023 

2024 

*MSET,  NSID  = 

65 

3002 

3003 

3004 

3005 

3006 

3007 

3008 

3009 

3010 

3011 

3012 

3013 

3014 

3015 

3016 

3017 

3018 

3019 

3020 

3021 

3022 

3023 

3024 

♦NSET,  NSID  = 

70 

4002 

4003 

4004 

4005 

4006 

4007 

4008 

4009 

4010 

4011 

4012 

4013 

4014 

4015 

4016 

4017 

4018 

4019 

4020 

4021 

4022 

4023 

4024 

♦ELEMENT,  LAYUP  =  1 


1 

1001 

1004 

1003 

1002 

1 

4 

3 

2 

2 

1006 

1005 

1004 

1001 

6 

5 

4 

1 

5 

1002 

1003 

1011 

1010 

2 

3 

11 

10 

6 

1003 

1013 

1012 

1011 

3 

13 

12 

11 

7 

1004 

1014 

1013 

1003 

4 

14 

13 

3 

8 

1005 

1015 

1014 

1004 

5 

15 

14 

4 

9 

1017 

1016 

1015 

1005 

17 

16 

15 

5 

10 

1018 

1017 

1005 

1006 

18 

17 

5 

6 

♦ELEMENT,  LAYUP  =  2 


65 

1058 

1059 

1075 

1074 

58 

59 

75 

74 

66 

1059 

1060 

1076 

1075 

59 

60 

76 

75 

67 

1060 

1078 

1077 

1076 

60 

78 

77 

76 

68 

1061 

1079 

1078 

1060 

61 

79 

78 

60 

69 

1062 

1080 

1079 

1061 

62 

80 

79 

61 

70 

1063 

1081 

1080 

1062 

63 

81 

80 

62 

71 

1064 

1082 

1081 

1063 

64 

82 

81 

63 

NODES 

CONTAINED 

ON  OUTER 

SIDE 

BOUNDARIES 

NODES 

ALONG  (0. 

.OOOE+00,0, 

,100E+02) 

<“> 

(0.200E+02,0 

.  lOOE+02) 

Figure  13.  Sample  format  of  file  restran.elg. 
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266 

267 

268 

269 

270 

271 

272 

273 

1266 

1267 

1268 

1269 

1270 

1271 

1272 

1273 

2266 

2267 

2268 

2269 

2270 

2271 

2272 

2273 

3266 

3267 

3268 

3269 

3270 

3271 

NODES  ALONG  (O.OOOE+00, 

-.lOOE+02)  <■ 

1 

V 

0 

0 

0 

0 

+ 

0 

0 

0 

.  lOOE+02) 

2 

6 

10 

18 

26 

34 

42 

50 

58 

66 

74 

86 

105 

120 

135 

150 

221 

236 

251 

266 

1001 

1002 

1006 

1010 

1018 

1026 

1034 

1042 

1050 

1058 

NODES  ALONG  (0. 

.200E+02, 

-.lOOE+02)  <• 

->  (0.200E+02,0 

.  lOOE+02) 

112 

127 

142 

157 

165 

173 

181 

189 

197 

205 

213 

228 

243 

258 

273 

1112 

1127 

1142 

1157 

1165 

1173 

1181 

1189 

1197 

1205 

1213 

1228 

1243 

1258 

1273 

2112 

2127 

2142 

2157 

2165 

2173 

2181 

2189 

2197 

2205 

Figure  13.  Sample  format  of  file  restran.elg  (continued). 


2.3.6  Material  Property  Input 

The  complete  specification  of  material  properties  encompasses  the  description  of  lamination  layup, 
elastic  moduli,  failure  criteria,  initial  material  damage,  and  property  degradation  laws.  The  input 
of  these  quantities  is  described  in  the  following  subsections. 


2.3.6. 1  Linear  Elastic  Material  Moduli 


The  assumption  of  a  layered  material  media  in  RESTRAN  requires  that  a  sequence  of  plies  be 
specified  even  in  the  case  of  a  homogeneous  material.  This  is  inputted  through  the  following  state¬ 
ment 


*LAYER,  LAYUP  =  M 
MATIDi,  THKi,  THETAj 

where  the  mandatory  parameter  LAYUP  =  M  designates  the  layer  identification  number  which 
is  matched  to  the  LAYUP  =  M  parameter  on  the  *ELEMENT  data  entry.  MATIDj  is 
the  identification  number  of  a  *MATERIAL  statement  which  gives  the  material  moduli  for  the 
i**  ply.  The  THKi  and  THETAj  entries  give  the  layer  thickness  and  ply  orientation  angle  with 
respect  to  the  global  or  local  coordinate  system,  respectively.  This  data  line  is  repeated  for  each 
ply  in  the  sublaminate. 

Each  layer  may  assume  different  orthotropic  material  properties  with  nine  independent  material 
moduli.  These  moduli  are  input  using  the  following  statement: 

♦MATERIAL,  MATID  =  Mj 
El,  E2,  E3,  Gi3,  G235  Gi2 
Vi3,  V23,  V12 

where  the  subscripts  (1,2,3)  refer  to  the  principle  fiber  coordinates  for  composite  materials  or 
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general  material  axes  for  homogeneous  orthotropic  materials. 


2. 3.6. 2  Nonlinear  Elastic  Material  Moduli 

The  support  of  nonlinear  material  behavior  in  RESTRAN  is  restricted  to  materials  which  exhibit  a 
nonlinear-elastic  stress-strain  relationship.  For  residual  strength  characterization,  nonlinear  inelas¬ 
tic  behavior  such  as  plastic  deformation  is  considered  a  failure  mechanism  and  may  be  accounted 
through  selected  failure  criteria  and  damage  laws.  Nonlinear  elastic  materials,  which  follow  the 
same  deformation  path  under  loading  and  unloading  as  shown  in  Figure  14,  are  indicated  through 
the  following  statement. 

♦NONLINEAR  MATERIAL 
ITLIM,  CVMTOL,  CVBTOL 
Pl»P2rP3»  •••  i'Pn 

where  ITLIM  is  the  maximum  number  of  iterations  permitted  to  obtain  a  converged  solution. 
CVMTOL  is  the  tolerance  of  the  difference  in  the  norm  of  the  residual  load  vector  between  iter¬ 
ations  below  which  convergence  of  material  failure  is  determined. 

llARill  <  CVMTOL  (5) 


CVBTOL  is  the  corresponding  prescribed  tolerance  for  predicting  convergence  of  buckling  eigen¬ 
values  such  that 


||Ai-Ai_i||  <  CVBTOL 


(6) 


The  *NONLINEAR  MATERIAL  statement  additionally  alerts  RESTRAN  to  treat  material 
properties  specified  on  the  *MATERIAL  entry  as  initial  values  and  to  subsequently  access  a 
user-defined  subroutine  denoted  USERNL  to  obtain  nonlinear  moduli  under  increasing  applied 
loads.  A  user-defined  subroutine  mode  of  input  was  selected  as  the  most  general  and  straightforward 
method  of  defining  nonlinear  elastic  properties.  As  opposed  to  tabular  entry  formats,  USERNL 
allows  all  material  moduli  to  be  varied  independently  provided  that  the  resulting  constitutive  ma¬ 
trix  remains  positive  definite.  Also,  specific  moduli  may  be  written  as  a  function  of  individual 
stress  or  strain  components,  or  based  on  combined  measures  of  unit  elastic  energy  or  strain  energy 
density.  RESTRAN  passes  interpolated  stress  and  strain  data  at  each  ply  transformed  into  the  local 
material  coordinate  system  of  the  current  element  and  requires  that  the  user  routine  return  scale 
factors  to  be  applied  to  the  initial  moduli  to  approximate  secant  material  properties.  The  input  of 
the  number  of  iterations  and  convergence  tolerance  parameters  is  optional,  and  default  values  of 
ITLIM  =  5,  CVMTOL  =  0.01,  and  CVBTOL  =  2  are  internally  assigned.  The  constants  P^v 
are  optional  input  parameters  that  are  passed  into  the  user-defined  subroutine  and  may  be  used  to 
define  secant  properties.  A  maximum  of  25  parameters  may  be  input  on  as  many  additional  lines 
as  needed. 
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The  subroutine  interface  written  in  FORTRAN  is  described  next.  The  ciurrent  implementation 
in  RESTRAN  assumes  a  single  nonlinear  material  that  can  be  described  by  a  maximum  of  25 
parameters  specified  in  the  input  file  (any  number  can  be  ’built-in’  to  the  user  written  subroutine). 
With  proper  consideration  of  the  argument  call  interface,  this  and  all  user-defined  subroutines  may 
be  developed  in  selected  languages,  compiled,  and  linked  with  the  main  RESTRAN  executable 
prior  to  job  execution.  Figure  15  shows  the  required  subroutine  format. 

SUBROUTINE  USERNL  (  STRESS, STRAIN, SXX,SYY,SZZ,TYZ,TZX,TXY,EXX,EYY, 

1  EZZ , GYZ , GZX , GXY , PARAM , THKN , ATHK , TLM , THETA , FEl , 

2  FE1,FE2,FE3,FG31,FG23,FG12,FV31,FV23,FV12,N0RD  ) 

C 

Q  *:^*3#s*s|t:*c*3(c3|c**:tc**5ft*:f:3|c**^:tc5fc******5tc*^****3|c*3|c*  ♦*♦♦*:»:** 

C  ** 

C  USER-DEFINED  ROUTINE  TO  COMPUTE  SECANT 

C  MODULI  IN  NONLINEAR  MATERIAL  ANALYSIS 

C  ♦♦  ** 

C 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

C 

DIMENSION  STRESS(27,6) ,STRAIN(27,6) ,PARAM(25) 

C 

***  USER  DEFINED  CODE  FOR  COMPUTING 

***  SECANT  MATERIAL  MODULI  *** 

C 

RETURN 

END 

Figure  15.  USERNL  subroutine  format. 
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This  routine  is  called  for  each  layer  in  each  element  for  the  calculation  of  material  properties  as  a 
function  of  stress/strain  state.  Stresses  and  strains  are  passed  to  the  user  subroutine  calculated  at 
Gauss  points  and  at  the  midpoint  of  each  element  layer  rotated  into  the  local  coordinate  system 
if  one  has  been  specified.  Gauss-point  tensor  quantities  are  stored  in  the  STRESS  and  STRAIN 
arrays  for  which  the  storage  order  depicted  in  Figure  16  is  the  same  for  both  stresses  and  strains. 
Interpolated  layer  stresses  are  provided  in  the  scalar  variables  SXX  TXY,  while  corresponding 
strain  values  are  stored  in  EXX  GXY.  The  total  laminate  and  layer  thicknesses  are  given  by  the 
TLM  and  ATHK  variables.  The  accumulated  thickness  to  the  current  ply  is  stored  in  the  THKN 
variable.  THETA  contains  the  orientation  angle  of  the  current  ply.  PARAM  is  an  array  containing 
input  parameters  to  define  the  material  nonlinearity.  The  user  routine  must  return  scalar  factors 
FEl  FV12  which  RESTRAN  uses  to  compute  secant  moduli  as  depicted  in  Figure  17  for  the 
current  load  level. 


where  Gi  = 

(^x  s'^xy  ) 

»^y  ’^z  ’Yyz’Yx2»Yxy) 


Figure  16.  Storage  format  for  stresses  and  strains. 


Figure  17.  Calculation  of  secant  moduli. 
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These  factors  are  applied  to  the  initial  moduli  as 


=  Fe\ 

jjk 

^31 

=  ^^31 

4lo 

—  Fe2 

rpk 

^2o 

^23 

=  Fv2Z 

4zo 

El 

=  Fe3 

El 

fjk 

^12 

=  Fv\2 

4o 

=  Fo23 

/^k 

^23o 

^31 

1 - 1 

CO 

II 

rik 

^Zlo 

r^k 

^12 

=  Fgi2 

rik 

^l2o 

This  user-defined  subroutine  must  be  compiled  and  linked  into  the  RESTRAN  executable  prior  to 
performing  residual  strength  analysis.  An  example  of  the  user- written  subroutine  USERNL  for 
returning  secant  moduli  is  presented  in  Appendix  A.  Finally,  it  must  be  noted  that  if  increased 
execution  time  in  converging  nonlinear  material  response  is  problematic  or  nonlinear  properties  are 
slight  or  inaccurately  characterized,  the  nonlinear  solution  algorithm  and  the  use  of  this  subroutine 
may  be  avoided  by  inputting  initial  moduli  as  effective  secant  approximations  based  on  measured 
ultimate  failure  strength  and  strain  data  as  depicted  in  Figure  18. 


Figure  18.  Linear  approximation  to  nonlinear  elastic  moduli. 

2.3.7  Failure  Mode  Selection 

For  residual  strength  prediction,  the  selection  of  the  METH  parameter  on  the  *SOLUTION  en¬ 
try  determines  whether  material  failure  modes  and/or  buckling  failure  modes  are  to  be  calculated 
in  the  algorithmic  assessment  of  progressive  failure. 

2. 3. 7.1  Material  Failure  Criteria 

RESTRAN  supports  a  robust  suite  of  material  failure  criteria.  The  specific  failure  criterion  to 
be  associated  with  an  input  material  labeled  Mi  is  specified  using 

♦FAILURE  CRITERION,  MATID  =  Mi 
FCTYPE 

Pi,  P2,  Ps)  P4>  Ps?  Pe 

P??  Psj  Pg?  Pro 
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where  PCTYPE  is  a  character  string  which  denotes  the  name  of  the  criterion,  and  Pi  through 
Pio  are  used  to  input  up  to  10  material  constants  required  by  the  criterion.  For  the  various  input 
parameters,  compressive  material  strengths  may  be  input  as  positive  or  negative,  the  sign  is  inter¬ 
nally  corrected.  The  criteria  designations  and  required  constants  are  listed  in  Table  2. 


Table  2:  Parameters  required  for  various  failure  criteria. 


Criterion 

Parameters 

MAX-STRESS 

Xt,  Xc,  Yt,  Yc,  Zt,  Zc,  R,  S,  T 

MAX-STRAIN 

Xc,  Yi,  Yc,  Tju  Rj  S,  T 

BELTRAMI 

^ult 

VON  MISES 

^ult 

HOFFMAN(l,2,3) 

Xcj  Yi,  Yc,  Zj,  Zc,  R,  S,  T 

HILL(1,2,3) 

X(,  Xc,  Yt,  Yc,  Zj,  Zc,  R,  S,  T 

TSAI-WU(1,2,3) 

Xt,  Xc,  Yt,  Yc,  Zt,  Zc,  R,  S,  T 

CHRISTENSEN 

X(,  Xc,  K,  a 

FENG 

Ai,  Aa,  As,  A4,  A5,  Ae 

HASHIN 

Xt,  Xc,  Yi,  Yc,  S 

USER-DEFINED 

NRCHK,  Pi  ^  Pio 

A  complete  description  of  the  various  failure  criteria,  including  the  required  material  constants, 
is  presented  in  the  RESTRAN  theoretical  manual  [3].  In  summary,  X^,  Xc,  Y^,  Yc,  Zf,  and  Zc  are 
the  normal  tensile  and  compressive  strengths  in  the  principle  1,  2,  and  3  directions,  respectively, 
and  jR,  5,  and  T  are  the  shear  strengths  defined  in  the  23,  13,  and  12  planes,  respectively.  For  the 
isotropic  Beltrami  and  Von  Mises  criteria,  is  the  combined  stress  mecisure  at  which  yielding  oc¬ 
curs  or  the  proportional  limit  reached.  The  constants  required  in  the  Christensen,  Feng,  and  Hashin 
criteria  are  explained  in  the  RESTRAN  Theoretical  Manual  [3].  The  variations  on  the  Hoffman, 
Hill,  and  Tsai-Wu  criteria  involve  specialization  of  these  criteria  to  account  for  two-dimensional 
or  three-dimensional  stress  states  and  whether  they  predict  a  single  unspecified  failure  mode  or 
dijfferentiate  between  fiber  and  matrix  modes.  The  variations  are  specified  by  the  appended  nu¬ 
merical  suffix  in  which  ’1’  designates  3-D  single  mode  failure,  ’2’  designates  2-D  plane  stress,  single 
mode  failure,  and  ’3’  specifies  3-D  mixed  mode  failure.  The  user-defined  option  is  discussed  in  the 
following  subsection. 


2.3. 7.2  User-Defined  Failure  Criteria 

Specifying  a  user-defined  failure  criteria  alerts  RESTRAN  to  access  a  user-supplied  subroutine 
designated  USERFC.  The  output  from  this  subroutine  is  an  unaltered  or  degraded  set  of  mate¬ 
rial  properties  for  the  current  ply  depending  on  whether  the  user-defined  criterion  has  predicted 
material  failure.  The  input  parameter  NRCHK  controls  error  checking  in  a  subroutine  driver  for 
USERFC.  If  NRCHK  is  equal  to  zero,  all  checks  will  be  bypassed.  A  nonzero  entry  will  cause  the 
following  tests  to  be  performed: 

1)  If  moduli  or  Poisson  ratios  have  increased,  a  warning  will  be  issued  and  execution  will  con- 
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tinue. 


2)  If  moduli  have  been  set  to  negative  values,  an  error  message  will  be  printed  and  execution 
will  continue  with  moduli  set  equal  to  zero. 

3)  If  scale  factors  are  returned  less  than  or  equal  to  zero,  an  error  message  will  be  printed  and 
an  absolute  value  of  the  scale  factor  will  be  used  in  subsequent  calculations. 

The  parameters  Pi  -)■  Pio  are  optional  and  may  be  used  to  pass  required  data  into  the  user- 
defined  subroutine.  The  format  of  USERFC  is  shown  in  Figure  19. 


SUBROUTINE  USERFC  (  SGPTS.EGPTS, STRESS, STRAIN, PTHK, THETA, IPLY, 

1  NELID , PARAM , El , E2 , E3 , G1 2 , G23 , G3 1 , V12 , V23 , 

2  V13,EM,EF,N0RD,NSTAT,NFSTAT,NACCS  ) 

C 

C 

C  **  USER-DEFINED  FAILURE  CRITERION  TO  COMPUTE  ** 

C  **  MATERIAL  FAILURE  UNDER  APPLIED  LOAD  ** 

C  **  ** 

C  ************* 

c 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

C 

DIMENSION  SGPTS(27,6) ,EGPTS(27,6) ,STRESS(6) ,STRAIN(6) ,PARAM(15) 
C 

***  SOURCE  CODE  FOR  USER-DEFINED  FAILURE  **=*= 

***  CRITERION  AND  DAMAGE  LAW  *** 

C 

RETURN 

END 


Figure  19.  USERFC  subroutine  format. 

The  input  and  output  data  passed  into  the  routine  through  the  argument  list  will  be  explained 
subsequently.  STRESS,  and  STRAIN  are  vectors  of  stress  and  strain  interpolated  to  the  cur¬ 
rent  ply  and  transformed  into  the  local  material  coordinate  system.  The  storage  convention  for 
the  stress/strain  vectors  is  given  by  [£^11,^22, £^33, £^23)  ■E'3i:-E’i2]-  SGPTS  and  EGPTS  are  arrays 
which  store  the  untransformed  stresses  and  strain  components  at  the  element  integration  or  Gauss 
points.  PTHK  and  THETA  are  ply  thickness  and  orientation  angle,  respectively.  IPLY  and  NELID 
are  ply  number  and  element  number,  respectively.  The  PARAM  array  contains  a  list  of  parameters 
from  the  *FAILURE  CRITERIA  entry.  Ei  — 1/31  are  material  moduli  and  Poisson  ratios.  EM,  and 
EF  are  scale  factors  for  matrix-  and  fiber-dominated  failure  modes.  NSTAT  is  the  ply  status  flag  in 
which  NSTAT  =  1  or  2  for  matrix-  or  fiber-dominated  failure  mode,  respectively,  and  NSTAT  =  3 
for  total  ply  failure.  In  addition,  the  option  exists  to  set  NSTAT  =  4,  which  causes  the  element  to 
be  excluded  from  further  degradation.  The  new  failure  status  is  returned  in  the  variable  NFSTAT. 
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NACCS  is  set  to  zero  if  no  failure  is  predicted  during  the  current  iteration  and  to  unity  if  failure  is 
predicted.  An  example  of  a  user-written  subroutine  for  determining  ply  failure  indices  is  presented 
in  Appendix  B. 

This  user-defined  subroutine  must  be  compiled  and  linked  into  the  RESTRAN  executable  prior 
to  performing  residual  strength  analysis. 


2.3.8  Material  Damage  Laws 

Independent  of  the  selected  failure  criterion,  several  types  of  material  damage  laws  may  be  invoked 
to  degrade  material  properties.  These  damage  laws  are  available  to  selectively  apply  varying  de¬ 
grees  of  material  degradation  to  simulate  progressive  ply  damage  modes.  For  complete  generality, 
a  user-written  subroutine  may  be  used. 


2.3.8. 1  Damage  Law  Selection 

A  specific  damage  law  is  selected  using  the  following  statement. 

♦DAMAGE  LAW,  MATID  =  Mi 
APPLIED-DAMAGE-LAW 
Pi5  P2?  P35  P45  P55  Pe 


where  MATID  =  M,  associates  the  damage  law  with  a  material  property  set  identification  number. 
The  following  parameter,  APPLIED-DAMAGE-LAW,  may  be  given  as  ’NULL’,  ’MIXED’, 
or  ’PARTIAL’  to  invoke  specific  damage  laws.  Information  regarding  failure  mode  determination 
by  the  selected  failure  criterion  is  stored  and  used  in  applying  the  selected  damage  law.  The  pa¬ 
rameters,  Pi,  P2,  P3,  P4j  Ps?  Per  may  be  used  in  the  computing  degraded  material  properties. 
Currently,  only  Pi  is  used  when  specifying  a  ’PARTIAL’  degradation  in  properties.  The  ’NULL’ 
option  causes  all  properties  of  the  current  ply  to  be  set  to  zero  as 


El  =  0.0 

Gzi  =  0.0 

1/12  =  0.0 

0 

0 

II 

0 

0 

11 

CO 

0 

b 

II 

CO 

(8) 

0 

0 

II 

» 

to 

II 

0 

b 

1^23  =  0-0 

When  selecting  a  single-mode  failure  criteria,  specifying 
be  degraded  according  to 

’PARTIAL’  will  cause 

properties  to  first 

El  =  aEi 

Gz\  =  olGzi 

Z/12  =  OIU12 

E2  =  OCE2 

C?23  =  OiG2Z 

Viz  =  OiViz 

(9) 

Ez  =  aEz 

Gi2  =  OiG\2 

^23  =  OiV2Z 

For  differentiating  between  fiber  and  matrix  failure  modes,  the  options  MIXED  or  PARTIAL 
may  be  selected.  Both  degrade  selected  material  properties  depending  on  predicted  failure  mode. 
For  a  partial  degradation,  a  reduction  factor,  a,  is  entered  through  the  Pi  parameter.  In  specifying 
a  mixed  mode  degradation,  the  factor  a  is  automatically  taken  as  zero.  The  degradation  laws  are 
given  by  the  following. 
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For  fiber  failure  mode 


E\  —  (xEi 

G31  =  aGzi 

Gi2  =  OiGi2 
vi2  =  avi2 

1/13  = 

For  matrix  dominated  failure 

E2  =  aE2 
E^  =  aE^ 
G25  —  <xG23 
^23  =  Q:1'23 

For  interlaminar  shearing  mode  failure 

G23  —  CtG23 

G31  =  aG3i 


(10) 


(11) 


(12) 


For  inplane  shearing  failure  mode 

Gi2  =  01G12  (13) 

After  initial  partial  ply  failure  is  accounted,  subsequent  failure  is  assumed  total. 

2. 3.8.2  User-Defined  Damage  Law 

A  user-defined  damage  law  is  applied  through  selecting  the  USER-DEFINED  failure  criterion 
on  the  *FAILURE  CRITERIA  entry.  The  USERFC  subroutine  detailed  in  Figure  19  is  then 
accessed  to  apply  the  specific  failure  criteria  and  return  reduced  material  properties  using  the 
user-defined  damage  law  encoded  into  this  subroutine. 

2.3.9  Initial  Damage  Description 

RESTRAN  has  been  specifically  designed  to  predict  residual  strength  of  composite  laminates  which 
contain  initial  damage  due  to  impact  events.  The  input  of  the  initial  damage  state  is  critical  in 
the  prediction  of  residual  strength  and  is,  perhaps,  the  most  difficult  data  to  obtain  from  a  user’s 
standpoint.  The  following  section  describes  the  input  requirements. 

2. 3.9.1  Initial  Material  Damage 

The  input  of  initial  damage  in  RESTRAN  is  performed  using  the  following  data  entry: 

♦INITIAL  DAMAGE 
ELSET,  DAMAGE-TYPE 
Pi)  P25  •••  )PlO 

where  ELSET  specifies  an  element  set  to  which  the  initial  damage  is  to  be  applied.  Damage 
is  assumed  to  apply  equally  to  all  plies  in  the  element.  Initial  property  degradation  may  be  applied 
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by  setting  the  the  DAMAGE-TYPE  parameter  equal  to  ’FIBER’,  ’MATRIX’,  ’TOTAL’, 
’USER-DEFINED’,  or  ’EXCLUDED’  to  specify  the  initial  effects  of  various  damage  modes. 
For  fiber  failure,  the  following  properties  are  degraded: 

£?i  =0 

C?31=0 

Gn  =  0  (14) 

Viz  =  0 

Vi2  =  0 

For  matrix  failure,  the  following  properties  are  degraded: 

^2=0 

Ez  =  0 

G23  =  0 

1^23  =  0 

For  total  failure,  all  material  properties  are  set  equal  to  zero  such  that 


El 

=  0.0 

G31 

=  0.0 

2^12 

=  0.0 

Ei 

=  0.0 

G23 

=  0.0 

=  0.0 

Ez 

=  0.0 

Gi2 

=  0.0 

1^23 

=  0.0 

Specifying  ’USER-DEFINED’  on  the  "^INITIAL  DAMAGE  entry,  the  input  parameters,  Pi 
Pio,  are  passed  into  the  user-defined  subroutine  USERID  for  initial  calculation  of  reduced  ma¬ 
terial  properties  to  simulate  initial  damage.  The  format  of  USERID  is  shown  in  Figure  20. 

SUBROUTINE  USERID  (PTHK, THETA, El, E2,E3,G13,G23,G12,V13,V23,V12, 

1  P1,P2,P3,P4,P5,P6,P7,P8,P9,NELID,IPLY,NSTAT  ) 

C 

Q  ****!)C  ♦*♦*♦****♦♦**♦**  **5*:******  St"****************** 


C  **  ** 
C  **  USER-DEFINED  SUBROUTINE  TO  ASSIGN  ** 
C  **  INITIAL  ELEMENT  MATERIAL  DAMAGE. 

C  **  ** 


Q  ♦  *:|c  ****♦*****:*!♦*  +  *****♦♦**=•'**=*=*’•'*  *******!|'!t'*>l'*!|‘!t:** 

C 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

C 

***  USER-DEFINED  INITIAL  MATERIAL  *** 

DAMAGE  SOURCE  CODE  *** 

C 

RETURN 

END 


Figure  20.  USERID  subroutine  format. 
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PTHK  and  THETA  are  the  current  ply  thickness  and  orientation  angle,  El  through  VI 2  are 
material  moduli,  PI  through  P9  are  input  parameters  provided  on  the  *MDAMAGE  data  entry, 
and  NELID  and  IPLY  are  the  current  element  and  ply  number.  NSTAT  is  a  ply  failure  status  flag 
which  is  set  to  T’  or  ’2’  for  matrix  or  fiber  (or  general)  failure,  respectively.  This  permits  additional 
material  failure  to  be  predicted.  If  total  failure  is  predicted,  NSTAT  is  set  equal  to  ’3’.  This  is 
usually  assigned  when  the  material  moduli  have  been  reduced  to  zero.  If  NSTAT  is  set  equal  to 
’4,’  then  the  current  damage  state  is  maintained  and  further  degradation  is  excluded.  This  may  be 
enforced  initially  through  a  damage- type  specification  discussed  next.  An  example  of  a  user- written 
subroutine  to  assign  initial  material  damage  is  contained  in  Appendix  C. 

The  specification  of  ’EXCLUDED’  is  provided  to  indicate  that  the  current  element  set  is  free 
of  initial  damage  and  is  excluded  from  accumulating  damage  during  subsequent  residual  strength 
analysis.  This  feature  permits,  in  perhaps  the  easiest  manner,  an  arbitrary  set  of  elements  to  be 
declared  failure  free  or  indestructible.  This  option  may  be  used  to  prevent  elements  from  failing  at 
which  external  loads  are  applied  -  which  will  normally  terminate  the  analysis  -  or  to  avoid  failing 
some  internal  structure  which  is  deemed  not  to  influence  residual  strength  predictions.  The  use  of 
this  method  for  protecting  elements  or  by  the  parameter  LAPFAIL  is  advocated  over  defining  a 
set  of  elements  with  high  material  moduli  or  strength  measures  because  these  approaches  could,  in 
uncommon  situations,  cause  these  exaggerated  element  properties  to  be  used  in  predicting  ultimate 
failure,  thus  yielding  incorrect  measures  of  residual  strength. 

2.3.9. 2  Delamination  Description 

A  delamination  is  defined  by  a  node  set  which  identifies  the  node  points  lying  within  the  delam¬ 
ination  plane.  RESTRAN  automatically  generates  a  corresponding  set  of  coincident  nodes  and 
redefines  element  connectivity  to  sever  connection  between  elements  lying  above  and  below  this 
plane  defined  by  the  normal  to  the  delamination  surface.  The  input  statement  used  to  define  de¬ 
lamination  planes  is  given  by 

*DELAMINATION 
NSETi,  NSETs,  NSET3,  ...  ,NSET^ 

where  NSET^  lists  node  set  identification  numbers  defined  in  the  *NODE  SET  entry.  Any 
number  of  node  sets  may  be  listed  per  line,  and  this  line  may  be  repeated  as  many  times  as  needed 
to  define  the  delaminations  present  in  the  model. 

RESTRAN  currently  calculates  buckling  loads  based  on  a  linear  stability  analysis.  To  approxi¬ 
mate  structural  behavior  in  the  post-buckled  regime,  modifications  to  the  material  properties  of 
those  elements  participating  in  the  buckling  mode  must  be  specified.  A  full  description  of  allowable 
post-buckled  element  response  is  discussed  in  the  RESTRAN  Theoretical  Manual  [3].  The  following 
statement  is  mandatory  for  residual  strength  analysis  which  includes  potential  sublaminate  buck¬ 
ling  failure.  The  entry  statement  is  given  by 

♦POSTBUCKLED  MATERIAL 
CLFAC,  PBFAC,  NITER 

where  CLFAC  specifies  whether  a  constant  load  at  buckling  failure  is  to  be  approximately  main¬ 
tained  in  the  buckled  elements  during  the  additional  progressive  failure  analysis.  This  parameter 
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is  a  real  number  ranging  from  0.0  to  1.0.  For  CLFAC  =  0.0,  no  internal  scaling  is  performed  to 
approximate  the  constant  load  state,  and  for  values  greater  than  0.0  up  to  1.0,  varying  degrees  of 
maintaining  the  internal  element  loads  present  at  failure  are  enforced.  The  parameter  PBFAC  is 
used  to  reduce  element  material  properties  to  account  for  buckling  failure.  This  factor  is  used  to 
reduce  the  constitutive  properties  as 

[Cy]  =  PBFAC[Cy]  (16) 

Both  CLFAC  and  PBMFAC  may  be  specified  to  account  for  assumed  post-buckled  sublaminate 
behavior.  For  CLFAC  >  0,  the  global  stiffness  matrix  becomes  a  nonlinear  function  of  loads  at 
which  prior  sublaminates  have  undergone  buckling  failure,  and  an  iterative  procedure  is  invoked 
to  converge  subsequent  eigen-analyses.  The  parameter  NITER  is  used  to  specify  the  maximum 
number  of  iterations  to  be  performed  in  calculating  further  sublaminate  buckling  loads. 


2.3.10  Coordinate  Transformation 

An  arbitrary  coordinate  system  {x,  y,  z)  may  be  specified  for  each  element  in  RESTRAN  to  define  a 
local  reference  system.  This  local  coordinate  system  may  be  used  to  define  ply  material  properties 
and  account  for  any  laminate  orientation  with  respect  to  the  element  (x’,y’,z’)  system  which  is 
aligned  with  the  global  coordinate  system.  On  output,  element  or  layer  stresses  and  strains  may 
be  output  in  this  local  system.  These  Cartesian  coordinate  systems  are  depicted  in  Figure  21. 


Figure  21.  Global  (x,y,z),  element  (x’,y’,z’)  and  local  {x,y,z)  coordinate  systems. 
A  local  coordinate  system  is  defined  using  the  following  input  statement 

♦ORIENTATION 

CID,  Vxi,  Vx2,  Vx3,  Vyi,  Vy2,  Vyg 
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where  CID  is  the  coordinate  system  number,  and  Vxi  -)■  Vya  are  the  components  of  two  vectors 
defined  in  the  global  coordinate  system  that  define  a  plane.  If  these  two  vectors  are  not  orthogonal, 
the  Vy  vector  is  made  orthogonal  to  Vx  and  the  third  coordinate  Vz  is  determined  from  the 
vector  cross  product  Vz  =  VxxVy.  The  coordinate  system  is  associated  with  selected  elements 
using  the  ORIENTATION  parameter  on  the  *ELEMENT  data  entry.  The  generation  of  a  local 
coordinate  system  is  depicted  in  Figure  22. 


Figure  22.  Local  coordinate  system  generation. 

To  set  a  tolerance  for  deformed  geometry  check,  the  following  statement  is  required: 

♦DEFORMED  GEOMETRY 
TOL 

where  TOL  is  the  angular  tolerance  in  degrees  such  that  warnings  wiU  be  issued  if  the  calcu¬ 
lated  internal  element  angles,  6,  are  outside  the  bounds 

TOL  <  0  <  180  -  TOL  (17) 


The  default  value  for  the  tolerance  is  45°.  This  will  only  alert  the  user  to  possibly  incorrect  element 
definitions  or  excessive  element  distortions  but  will  not  stop  execution. 

Nodes  may  be  equivalenced  through  the  following  command 

♦EQUIVALENCE 
NSETl,  NSET2 

where  NSETl  and  NSET2  are  node  set  identifiers.  Any  number  of  node  set  pairs  may  be  entered 
on  separate  lines  and  must  be  of  equal  size.  This  command  specifies  that  each  node  in  sequence  in 
the  first  node  set  is  to  be  replaced  by  the  associated  node  in  the  second  set.  This  operation  may 
be  used  to  join  separate  portions  of  a  finite  element  model. 
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Elements  may  be  excluded  from  failure  using  the  following  general  statement 

♦EXCLUDE  ELEMENT 
ELSETl,  ELSET2,  ...  ,  ELSETn 

where  ELSETi  axe  element  set  identifiers.  Any  number  of  element  sets  may  be  entered  on  each 
line.  The  total  number  of  input  sets  is  currently  limited  to  100.  This  command  specifies  that  all 
elements  specified  in  the  various  sets  be  excluded  from  failure.  This  feature  may  be  used  for  models 
where  only  failure  in  a  portion  of  the  model  is  of  interest,  or  regions  of  the  overall  model  may  be 
coarsely  modeled  and  local  failure  in  these  regions  is  not  to  be  accounted  for  in  the  prediction  of 
residual  strength. 

To  terminate  parsing  of  the  input  file,  the  following  statement  is  required: 

♦ENDDATA 

Any  statements  following  *ENDDATA  in  the  input  file  are  ignored. 

2.4  Execution  Modification 

Various  modifications  to  the  execution  of  solution  procedmres  may  be  specified  through  the  use  of 
the  *PARAMETER  statement. 

2.4.1  User-Defined  Buckling  Mode  Interpretation 

RESTRAN  has  built-in  algorithms  for  processing  and  interpreting  sublaminate  buckling  modes. 
These  algorithms,  however,  assume  that  delaminations  form  a  planar  surface.  As  described  in  the 
RESTRAN  users  manual  [3],  this  surface  need  not  be  continuous,  but  a  delamination  which  cannot 
be  described  by  a  sequence  of  plane  sections  cannot  be  processed.  To  overcome  this  limitation,  a 
user-written  subroutine  may  be  used  to  analyze  the  buckling  mode  of  complicated  delamination 
configurations.  Specialized  interpretation  of  simultaneous  buckling  of  multiple  delaminations  or 
stepped  fracture  surfaces,  such  as  depicted  in  Figure  23,  may  then  be  performed. 


Figure  23.  Complex  delamination  states. 
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To  invoke  the  access  of  a  user-written  subroutine  to  analyze  buckling  modes,  the  following  en¬ 
try  is  included  anywhere  in  the  input  stream. 

♦PARAMETER  USRBKL 

This  causes  RESTRAN  to  bypass  the  existing  interpretive  algorithms  and  access  two  user-defined 
routines.  These  subroutine  are  designated  USERBl  and  USERB2.  The  routine  USERBl  is 
accessed  to  analyze  the  current  buckling  mode  shape  and  determine  whether  the  buckling  mode  is 
acceptable  and  should  be  used  as  a  possible  instability  failure  mode.  The  format  of  this  routine  is 
given  in  Figure  24. 


C 

C 

C 

C 

C 

C 

C 

C 

C 


C 

C 

C 


SUBROUTINE  USERBl  (  EMODE,NELEM,NDSET,NXSET,NDID,EIGV,NUMNOD, 

1  NBPRC ,LCNSL , ITERB , NUMEXC ,NACC , NSURF ,NINC , 

2  MDX0,MDX1,MAXB,MDX3  ) 


**  USER-DEFINED  SUBROUTINE  TO  EVALUATE  ** 

**  BUCKLING  FAILURE  MODE  ** 


*****************************  J(c******:)t*******!(t**** 

IMPLICIT  DOUBLE  PRECISION (A-H,0-Z) 

DIMENSION  EMODE(MDXO,MAXB) ,NELEM(MDX3,12) ,NDSET(MDX1) 
DIMENSION  NXSET(MDXl) 

SOURCE  CODE  FOR  EVALUATING  BUCKLING  FAILURE  *** 

RETURN 

END 


Figure  24.  USERBl  subroutine  format. 


This  routine  is  passed  the  mode  shape  contained  in  the  emode  array  stored  as 


NODE 

POSITIVE  NORMAL 

NODE 

NEGATIVE  NORMAL 

ID 

U 

V 

w 

ID 

U 

V 

W 

Ni 

Ui 

Wi 

Ml 

Ui 

Vi 

Wi 

N2 

U2 

V2 

W2 

Ma 

U2 

V2 

W2 

Ns 

Us 

Vs 

Ws 

Ms 

Us 

Vs 

Ws 

Nn 

Un 

Vn 

Wn 

Mn 

Un 

Vn 

Wn 
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where  the  displacements  are  stored  in  each  record  as  coincident  node  pairs,  Ni  Sc  Mi,  corresponding 
to  the  positive  and  negative  normal  to  the  delamination  surface.  NDSET  is  a  linear  array  contain¬ 
ing  node  numbers  defining  the  current  delamination,  NXSET  is  an  array  to  store  NUMEXC  nodes 
in  the  delamination  set  at  which  coincident  nodes  are  to  be  removed  to  satisfy  contact  constraints. 
NDID  is  the  node  set  ID,  EIGV  is  the  eigenvalue  of  the  buckling  mode,  NODNUM  is  the  total 
number  of  nodes  in  the  model,  NBPRC  and  ITERB  store  the  buckling  mode  interpretation  proce¬ 
dure  fiag  and  number  of  iterations,  and  LCNSL  is  the  current  iteration  count  for  the  delamination. 
NACC  is  set  equal  to  zero  if  the  mode  shape  is  unacceptable  and  equal  to  one  if  this  buckling 
mode  is  to  be  considered  in  determining  the  next  failure  mode.  NSURF  is  a  parameter  that  will  be 
passed  into  the  user-written  subroutine  USERB2  which  will  assess  element  failures  in  the  buckled 
sublaminate.  The  format  of  the  USERB2  user-defined  subroutine  is  given  in  Figure  25. 


SUBROUTINE  USERB2  (  NELEM,NODSET,NSETI,NDELM,NDID,NELFL, 

1  ELFAC , NFAIL , NUMELF , NSURF , MAXI , MAX2 , 

2  MAX3,MAX4  ) 

C 

C  ***********:)c****************************!(c*****J|c:|c* 


C  **  ** 

C  **  USER  DEFINED  SUBROUTINE  TO  DETERMINE 

C  **  ELEMENT  FAILURE  DUE  TO  BUCKLING  ** 

C  **  ** 


C 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

C 

DIMENSION  NELEM(MAX3,12) ,N0DSET(MAX2) 
DIMENSION  NELFL(MAXl) , ELFAC (MAXI) ,NDELM(MAX4) 
C 

USER  DEFINED  CODE  TO  SET  ELEMENT  *** 

***  FAILURE  DUE  TO  BUCKLING 
C 

RETURN 

END 


Figure  25.  USERB2  subroutine  format. 

The  array  NELEM  contains  the  element  connectivity  data.  The  format  of  this  array  is 


NELID 

N1 

N2 

N3 

N4 

N5 

N6 

N7 

N8 

in  which  NELID  is  the  element  number,  and  Nl->  N8  are  the  node  numbers.  These  node  numbers 
have  been  renumbered  internally.  To  obtain  the  input  node  numbers,  the  conversion  table,  NSETI, 
can  be  used.  The  original  input  node  number  for  the  node  and  the  element  is  obtained  as 
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MPDE  =  MSETI(NELEM(I,J+1)) 

The  nodset  array  contains  the  node  ID’s  defining  the  current  delamination.  This  array  has  the 
form: 


NODSET(NSID,  Nl,  N2,  ...  ,Nn  ) 


where  NSID  is  the  node  set  ID,  and  Ni  axe  the  node  numbers.  The  routine  assigns  element  failure 
such  that  NELFM(I)  stores  the  element  ID,  and  the  property  reduction  factor  for  this  element 
is  stored  in  ELFAC(I).  NDID  is  the  delamination  ID  number.  NUMELF  is  the  total  number  of 
elements  that  have  been  determined  to  fail,  which  is  set  by  the  user  subroutine. 

Examples  of  user-written  subroutines  to  interpret  buckling  modes  and  assign  element  failures  are 
presented  in  Appendices  D  and  E. 


2.4.2  Element  Failure  Exclusion 
♦PARAMETER  LAPFAIL 

This  parameter  causes  the  elements  to  which  external  loads  are  applied  to  not  be  analyzed  for 
failure.  If  an  element  is  completely  failed  to  which  external  forces  are  applied,  the  definition  of 
residual  strength  used  in  RESTRAN  involving  a  scalar  multiplier  to  the  initial  load  vector  loses 
definition  and  the  analysis  is  terminated.  Although  these  particular  elements  may  be  identified 
and  assigned  moduli  and/or  strength  ultimates  to  preclude  failure,  in  certain  cases  where  elements 
in  this  set  are  involved  in  the  final  catastrophic  failure  event,  RESTRAN  will  use  the  material 
constants  associated  with  these  elements  and  predict  an  erroneously  large  strength  measure  before 
terminating.  Thus,  the  use  of  this  parameter  is  generally  recommended  but  is  maintained  as  an 
option. 


2.4.3  Job  Execution  Status 

In  running  an  iterative  solution  sequence  which  may  involve  lengthy  execution  times  for  large  mod¬ 
els,  the  following  entry  can  be  included  to  request  an  intermediate  record  of  where  execution  is 
proceeding 

♦PARAMETER  STATUS 

The  additional  output  containing  execution  description  and  timing  information  is  written  to  the 
file  restran.log. 
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2.5  Output  Control 

Depending  on  the  operating  system  used  to  run  RESTRAN,  partitions  set  up  for  scratch  files  may 
be  inadequate  to  hold  the  temporary  files  generated  during  solution.  A  parameter  may  be  input  to 
direct  the  directory  in  which  scratch  files  are  to  be  stored.  This  parameter  is  given  by 

^PARAMETER  DIRECTORY  =  FPATH 

where  FPATH  may  be  set  equal  to  CURRENT  or  SCRATCH.  If  the  CURRENT  option  is 
used,  all  temporary  files  are  created  in  the  current  directory  in  the  fort.*  format  and  must  be 
manually  deleted.  The  default  is  SCRATCH. 

Nodal  displacement  output  is  specified  using  the  following  statement: 

♦NODE  PRINT,  NODSET  =  NDS 
NSTRING 

where  NSTRING  is  a  character  string  containing  the  characters  D  and/or  M  in  any  order  to 
print  displacements  and/or  buckling  mode  shapes,  respectively,  after  convergence  of  the  i*^  load 
increment.  The  optional  parameter  NDS  may  be  used  to  specify  a  node  set  number  to  limit  output. 

Output  of  element  stress  and  strain  data  is  requested  using  the  statement 

♦ELEMENT  PRINT,  OUTPUT=LOC,  SYSTEM=CRD,  ELSET=NES 
ESTRING 

where  ESTRING  is  a  character  string  containing  the  characters  S  and/or  E  in  any  order  to 
print  stresses  and/or  strains,  respectively.  The  optional  character  parameter  LOC  may  be  set  to 
’GAUSS’  (default)  or  ’LAYER’  for  output  at  element  Gauss  points  or  at  the  midpoint  of  ply  layers, 
respectively.  The  character  parameter  CRD  may  be  set  to  ’LOCAL’  for  output  in  the  local  element 
coordinate  system  or  ’GLOBAL’  (default)  for  output  in  the  global  coordinate  system.  The  integer 
parameter  NES  may  be  used  to  specify  an  element  set  to  limit  output. 

In  an  iterative  residual  strength  analysis,  ply  failure  statistics  can  be  output  after  convergence 
at  each  analysis  cycle  through  the  following  statement: 

♦PLY  FAILURE  PRINT,  NELSET=NPS 

where  ply  failures  in  each  element  are  output  showing  the  percentages  of  failures  due  to  ma¬ 
trix,  fiber,  general,  and  buckling  failure.  The  optional  parameter  NPS  may  be  used  to  specify  an 
element  set  number  to  limit  output.  This  statement  may  also  be  used  in  a  linear  static  analysis 
to  cause  the  computation  of  element  failure  indices  after  stress  recovery  is  completed.  Output  is 
limited  to  the  element  ply  exhibiting  the  maximum  failure  index,  the  predicted  failure  mode,  and 
the  associated  factor  to  scale  the  applied  loads  to  cause  failure  of  the  ply.  At  the  end  of  the  list, 
the  minimum  scale  factor  is  identified  and  may  be  considered  a  residual  strength  measure  based  on 
an  overall  first-ply  failure  criterion. 


34 


2.6  Graphics 

The  graphics  option  in  RESTRAN  has  been  developed  with  two  primary  aims.  First,  because  of 
the  specialized  nature  of  RESTRAN  in  predicting  progressive  failure,  the  graphical  information 
of  interest  is  assumed  to  concentrate  on  the  element  level  failure  depictions  together  with  overall 
deformation  and/or  local  buckling  mode  shapes.  This  consideration  limited  the  need  to  provide 
detailed  graphical  information  on  local  stress/strain/deformation  contours  as  is  available  using 
commercial  pre-  and  postprocessing  programs  such  as  PATRAN,  which  are  designed  for  general 
structural  analysis.  Secondly,  a  commercial  software  package  in  widespread  use  with  graphical 
capabilities  and  available  for  use  on  computer  platforms  down  to  the  PC  level  was  also  of  interest. 
This  would  free  the  essential  operations  of  the  RESTRAN  analysis  program  from  being  linked  to 
built-in,  platform-dependent  graphics  routines  which  would  require  several  versions  of  the  program 
to  be  released.  Therefore,  graphics  output  can  be  selected  in  two  formats:  MATHEMATICA  and 
TECPLOT.  Both  are  common  commercial  codes  which  offer  a  robust  capability  in  processing  prim¬ 
itive  graphical  objects  such  as  polygonal  surfaces  and  wireframe  meshes.  The  generated  output  file 
can  be  directly  input  into  the  selected  external  code  for  viewing.  The  MATHEMATICA  format 
offers  more  parameters  for  viewing  the  input  structural  model  and  progressive  failure  through  a 
color  coding  and  gray-scale  scheme.  This  input  file  contains  complete  information  on  the  structural 
geometry  and  can  be  edited  to  change  any  viewing  parameters  such  as  those  effecting  view  ori¬ 
entation,  lighting,  and  magnification  factor  for  deformation/mode  shape  depiction.  The  graphics 
option  is  invoked  by  including  the  following  statement  in  the  input  file: 

^GRAPHICS,  FORMAT  =  FMT 
NDFRM,  NWTRN,  GMFAC 

where  the  parameter,  FMT,  is  set  equal  to  ’MATHEMATICA’,  ’TECPLOT’,  or  ’USER- 
DEFINED’  to  specify  the  required  format  of  the  output  file.  Additional  optional  parameters  are 
given  by 


NDFRM  =  1  For  only  undeformed  geometry 

=  2  For  both  deformed  and  undeformed  geometry 
NWTRN  =  1  For  a  wireframe  depiction  in  PREPASS 
=  2  For  surface  depiction  in  PREPASS 
GMFAC  =  Magnification  factor  for  deformed  plots 


In  Prepass,  RESTRAN  will  generate  a  solid  or  wireframe  depiction  of  the  model,  and  output  is 
written  to  a  file  designated  restran.prp.  In  any  other  solution  option,  RESTRAN  will  generate 
deformation/mode  shape  graphics  for  linear  static  or  instability  analysis.  For  residual  strength 
analysis,  a  visual  depiction  of  damage  state  and  deformation  plot  will  be  generated  at  the  end  of 
each  converged  progressive  failure  cycle.  This  information  is  stored  in  a  file  designated  restran.grf. 
These  files  may  then  be  directly  input  into  the  selected  graphics  program  to  generate  the  graph¬ 
ical  display.  An  example  of  a  generated  graphics  file  in  MATHEMATICA  format  is  presented  in 
Appendix  F. 

The  ’’USER-DEFINED”  option  causes  RESTRAN  to  call  the  user- written  subroutine  USRGRF. 
The  format  of  this  subroutine  is  shown  in  Figure  26. 
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SUBROUTINE  USRGRF  (NELEM, COORDS, U,NSETI,NELNUM,NUMNOD, 

1  NWTRN , NDFRM , GMFAC , SCALE , NPATH , NTG , 

2  NFLAG,NDEV,MDX1,MDX2,MDX3  ) 

C 

Q  ****:  +  *,|c*****:+;*:)c*  +  **  +  +  **  +  *:(c******!t:*********  +  * 

C  **  ** 

C  **  USER-DEFINED  GRAPHICS  OUTPUT  ROUTINE  ** 

C  **  ** 

Q  *********************  *********** 

C 

IMPLICIT  DOUBLE  PRECISION (A-H,0-Z) 

C 

DIMENSION  NELEM (MDX3, 12) , COORDS (MDX2, 4) ,U(MDX1) 
DIMENSION  NSETI(MDXl) 

C 

C 

C  NODE  CONVENTION 

C 


c 

Y,V 

c 

2,W 

/ 

c 

1 

/ 

SURFACE 

c 

1  8 — 

- 7 

FACI 

NODES 

c 

1  /I 

/I 

FI 

(1,2, 6, 5) 

c 

1  /  1 

/  I 

F2 

(2, 3, 7, 6) 

c 

1/  1 

/  1 

F3 

(3, 4, 8, 7) 

c 

5 - 

— 6  1 

F4 

(1,5, 8, 4) 

c 

1  1 

1  1 

F5 

(1,4, 3, 2) 

c 

1  4 — 

—  1 - 3 

F6 

(5, 6, 7, 8) 

c 

1  / 

1  / 

c 

1  / 

1  / 

c 

1/ 

1/ 

c 

1 - 

- — 2 - 

—  X,U 

C 

C 

C  «<  PROGRAM  STATEMENTS  TO  OUTPUT  GRAPHICAL  »> 

C  «<  INFORMATION  IN  A  USER-DEFINED  FORMAT  »> 

C 

RETURN 

END 

Figure  26.  Format  of  user-defined  subroutine  USRGRF. 

In  USRGRF,  the  storage  format  of  element  connectivity,  coordinate  data,  and  node  number  con¬ 
version  contained  in  the  NELEM,  COORDS,  and  NSETI  arrays  have  been  previously  defined  in 
Section  2.4.1.  Other  quantities  passed  into  the  user-defined  routine  are  the  U  array  which  contains 
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the  displacement  or  mode  shape  vector  and  a  sequence  of  scalar  parameters.  NELNUM  and  NOD- 
NUM  contain  the  number  of  elements  and  nodes,  respectively.  SCALE  and  NTG  store  the  current 
load  multiplier  and  global  analysis  cycle  number.  NPATH  stores  the  solution  procedure  currently 
being  processed.  This  is  given  by 


NPATH  =  i 


1 

2 

4 

5 


MaterialFailureCalculation 
BucklingFailureCalculation 
Linear  Static  Analysis 
Linear  Buckling  Analysis 


(18) 


NFLAG  =  1  or  2,  depending  on  whether  the  solution  is  in  a  PREPASS  or  ANALYSIS  mode,  respec¬ 
tively.  The  remaining  parameters,  MDXl,  MDX2,  and  MDX3  are  size  parameters  set  internally 
by  RESTRAN  and  must  not  be  altered.  An  example  of  a  USRGRF  subroutine  is  presented  in 
Appendix  G. 


2.7  Miscellaneous 

To  include  a  heading  or  title  to  the  program  output,  the  following  input  record  is  used: 


♦HEADING 
JOB  TITLE 


where  the  job  title  may  be  up  to  80  characters  in  length. 

For  convenience,  the  input  file  may  be  included  with  the  program  output  by  including  the  fol¬ 
lowing  statement 

♦ECHO 


This  will  echo  the  input  file  verbatim  to  the  output  stream. 

Comments  may  be  inserted  between  statements  anywhere  in  the  input  file  using  a  ’♦*’  format 
in  fields  1  and  2.  For  example 

♦♦  This  line  is  a  comment 

will  cause  the  input  file  interpreter  to  skip  this  line.  Comments  may  also  be  inserted  within 
the  range  of  input  statements  that  have  a  fixed  number  of  entries  per  record  and  can  have  an  arbi¬ 
trary  number  of  records  or  lines.  These  include  the  ♦NODE,  ♦ELEMENT,  and  ♦BOUNDARY 
statements. 

Model  data  created  in  programs  such  as  PATRAN  in  ABAQUS  format  may  be  directly  included 
by  using  the  following  statement 

♦INCLUDE,  FILE=’filename.inp’ 

This  will  cause  input  to  be  read  from  the  specified  file.  The  filename  must  be  enclosed  in  quotes. 
The  included  file  may  contain  any  input  statements  that  are  equivalent  to  ABAQUS  format  such 
as  ♦NODE,  ♦ELEMENT,  ♦BOUNDARY,  or  ♦CLOAD.  All  optional  parameter  inputs  should 
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be  in  accordance  with  RESTRAN  input  formats  and,  in  addition,  the  *ENDDATA  or  *END  STEP 
line  must  be  present  in  the  included  file  to  return  lO  processing  to  the  standard  RESTRAN  input 
file  restran.inp. 

3  RESTRAN  Error  Checking 

An  extensive  set  of  internal  testing  of  initial  input  and  run-time  operations  has  been  incorporated 
into  RESTRAN  to  alert  the  user  to  inconsistencies  and  errors  encountered  in  the  model  generation 
and  execution  of  the  various  RESTRAN  analysis  options.  The  range  of  error  testing  is  described 
in  the  following  subsections. 

3.1  Tests  on  Initial  Input 

The  essential  tests  on  the  basic  input  involve  checking  whether  the  RESTRAN  size  parameters  are 
adequate  to  process  the  current  model  and  verify  that  a  valid  model  has  been  input.  Testing  of 
model  validity  includes  checking  node  definition,  element  connectivity,  assigned  material  properties, 
support  conditions,  and  failure  criteria.  These  tests  are  designed  to  ensure  that  a  viable  model  can 
be  created  for  performing  subsequent  analyses.  Fatal  errors  are  determined  as  those  that  preclude 
the  subsequent  operation  of  an  analysis  algorithm  and  will  cause  program  execution  to  cease. 
Warnings  are  issued  if  inconsistencies  are  detected,  such  as  unused  nodes,  but  execution  is  allowed 
to  proceed. 

3.2  Run-Time  Error  Determination 

A  valid,  consistent  structural  model  may  still  be  ill-posed  and  produce  meaningless  results.  During 
execution,  numerous  tests  axe  made  to  assess  elastic  and  differential  stiffness  properties,  eigenmode 
determination  and  interpretation,  material  degradation,  and  processing  performed  by  user-defined 
subroutines.  Although  any  set  of  error  checking  is  bound  to  be  incomplete,  the  internal  testing 
in  RESTRAN  will  alert  the  user  to  the  most  simple  and  to  many  subtle  errors  in  the  analysis  of 
residual  strength. 

4  Input  File  Format 

The  RESTRAN  input  file  is  designated  restran.inp.  The  file  format  has  been  implemented  with 
a  degree  of  flexibility  in  the  parsing  of  input  such  that  all  statements  are  case  insensitive,  and 
arbitrary  space  or  comma  delimitators  can  be  used  in  listing  statement  options  and  data  items. 
The  input  stream  in  RESTRAN  is  parsed  sequentially  and  may  generally  be  input  in  any  order.  The 
only  exceptions  are  the  node  and  element  generation  statements  *NGEN,  *NFILL,  *ELGEN, 
and  *MODEL  GENERATION,  which  are  invoked  during  the  input  parsing.  The  routines  that 
are  invoked  require  node  or  element  data  to  have  been  read  prior  to  performing  the  requested  node 
or  element  generation.  Termination  of  input  is  specified  by  the  *ENDDATA  entry. 

5  Computer  Implementation  of  RESTRAN 

RESTRAN  is  written  in  FORTRAN  77  with  the  aim  of  providing  maximum  portability  across 
different  computer  platforms.  There  is  no  overriding  driver  or  script  to  automate  various  support 
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or  implementation  functions,  such  as  setting  array  size  parameters  and  processing  user-defined  sub¬ 
routines.  The  installation  of  RESTRAN  requires  that  a  master  memory  size  parameter  be  set  in 
the  preface  to  the  MAIN  section  of  the  RESTRAN  code.  This  parameter  is  denoted  MAXRAM 
and  represents  the  maximum  number  of  double  precision  words  that  can  be  held  in  core  memory. 
This  segment  of  code  is  shown  in  Figure  27. 

C  ®@(§ 

C  @@@  RESTRAN,  REsidual  STRength  ANalysis,  IS  A  FINITE  @@@ 

C  @@@  ELEMENT  BASED  PROGRAM  TO  PREDICT  THE  RESIDUAL  @®@ 

C  @®@  STRENGTH  OF  IMPACT  DAMAGED  COMPOSITE  LAMINATES  ®®® 

C  ®@@  ®®® 

C  ®®®®®®®@®®®@®®®®®®®®®®®®®®®®@@@®®®@®@®®®®@@®@®®®®®@®®®®®®®@ 

C 

IMPLICIT  DOUBLE  PRECISI0N(A-H,0-Z) 

C 

C  ®®@@®®®®®®®®®®®®®®®®®®@®@@®®®®®®@®@@®@®©®@@@@ 


c  @@®  @®® 

C  ®®®  RESTRAN  MEMORY  PARMETERS  ®@® 
C  ®®®  ®®® 
C  ®@®@@@@©®®@@®®®®@@@®@@@®®®®®®®®®@@@®@@®@®®@@® 
C 

C  «<  SIZE  PARAMETERS  USED  IN  PROGRAM:  »> 
C  «<  »> 
C  «<  MAXRAM:  MAXIMUM  AVAILABLE  RAM  STORAGE  IN  »> 
C  «<  DOUBLE  PRECISION  WORDS  »> 
C 

PARAMETER  (  MAXRAM  =  16500000  ) 


Figure  27.  MAXRAM  parameter  in  RESTRAN. 

During  program  execution,  judicious  segments  of  RAM  memory  are  internally  allocated  to  perform 
the  preface  operations  involving  model  input  and  processing.  Minimum  array  sizes  are  automati¬ 
cally  determined,  and  memory  is  reallocated  to  minimize  both  internal  storage  requirements  and 
the  size  of  external  disk  files  to  streamline  the  execution  of  a  RESTRAN  analysis. 

User-defined  subroutines  must  be  created  using  a  standardized  subroutine  parameter  interface. 
For  convenience,  sample  routines  are  provided  in  separate  appendices.  These  routines  need  not 
be  written  in  FORTRAN,  but  if  other  languages  are  used,  care  must  be  given  to  differences  in 
conventions  such  as  how  arrays  are  ordered  as  in  C.  These  routines  must  then  be  compiled  and 
linked  with  the  RESTRAN  object  code  to  form  an  executable  module. 

6  Computational  Efficiency 

The  computational  efficiency  of  RESTRAN  is  primarily  dictated  by  the  matrix  operations  associ¬ 
ated  with  decomposition  and  eigenvalue  extraction.  As  discussed  in  more  detail  in  the  RESTRAN 
Theoretical  Manual  [3],  the  finite  element  basis  naturally  leads  to  the  generation  of  large,  sparse. 
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banded  matrices  in  representing  the  complete  model.  Internal  algorithms  have  been  created  to  pro¬ 
cess  matrix  storage  modes  in  various  formats.  These  include  full  matrix  storage  for  small  problems, 
half-bandwidth  storage  format  (default),  and  out-of-core  storage  for  large  problems  that  exceed  the 
internal  RAM  memory  capacity  and  require  most  data  storage  to  be  resident  in  external  disk  files. 
The  basic  determinator  of  execution  speed  is  the  selection  of  internal  data  storage  format  required 
to  hold  sufficient  data  in  high-speed  in-core  memory.  Bandwidth  storage  assumes  enough  RAM 
memory  to  store  a  complete  global  matrix  at  any  one  time  with  intermittent  I/O  operations  to 
swap  one  global  matrix  quantity  for  another.  The  out-of-core  solution  mode  is  slowest  due  to  the 
high  I/O  overhead  of  transferring  data  between  core  and  external  files  which  must  be  performed 
constantly  during  program  execution. 

Additional  considerations  pertain  to  the  direct  assembly  procedure  currently  used,  which  is  more 
computationally  and  I/O  costly  than  frontal  solution  methods,  and  to  the  condensation  of  internal 
equations  involved  in  isolating  individual  delaminations  and  imposition  of  applied  displacement 
boundary  conditions.  These  constraints  are  applied  using  row  and  column  shifting,  which  is  more 
intensive  than  simple  equation  isolation  using  unit  diagonal-zeroed  row/column  procedures  that  are 
straightforward  in  linear  equation  solution  but  introduce  difficulties  in  performing  eigenanalyses. 

Thus,  the  execution  time  of  performing  an  analysis  using  RESTRAN  is  significantly  higher  than 
that  obtained  in  performing  similar  solution  procedures  in  the  widely  used  commercial  codes  such 
as  NASTRAN,  ABAQUS,  and  ANSYS,  and  may  be  addressed  in  future  enhancements.  The  spe¬ 
cialized  analysis  of  residual  strength  accounting  for  combined  progressive  material  and  structural 
failure  performed  by  the  RESTRAN  code  is,  however,  simply  not  available  in  a  single  commercial 
program. 

7  Numerical  Examples 

The  following  subsections  present  sample  problems  illustrating  the  use  of  the  RESTRAN  code. 
Each  problem  uses  a  model  of  a  simple  laminated  composite  casing  in  order  to  contrast  different 
finite  element  discretizations,  applied  loading,  existing  damage,  and  analysis  procedures.  Timing 
summaries  are  provided  to  demonstrate  the  performance  of  RESTRAN  in  a  Linux  environment  on 
a  Carrera  workstation  with  a  530  MHz  Alpha  cpu. 

7.1  Linear  Static  Analysis 

A  static  analysis  of  the  finite  element  model  shown  in  Figure  25  is  performed.  An  arbitrary  layup 
based  on  a  nominal  ply  thickness  of  0.0052in  was  selected  for  the  various  model  sections.  The  top 
and  bottom  plates  were  composed  of  (906/i454  |  ±45/028  !)« laminates  where  ’|’  designates  element 
layer  interfaces.  The  vertical  panel  laminate  was  modeled  as  (Og/ ±454/905/ ±455  I  Oe/ ±455/04  l)^. 
Ply  properties  were  selected  as 


El 

=  2.48E6 

E2 

=  2.48E6 

Ez 

=  0.71E6 

G23 

=  3.0E5 

Gi3 

=  3.0E5 

Gi2 

=  9.9E5 

=  0.28 

1^13 

=  0.28 

=  0.243 
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For  first  ply  failure,  the  maximum  stress  criterion  was  selected  with  the  following  material  strengths 


Xt  =  21.3E3 

Yt 

=  21.3E3 

Zi  =  21.2E3 

Xc  =  47.8E3 

Yc 

=  47.8E3 

Zc  =  21.0E3 

R  =  1.65E3 

S 

=  1.65E3 

T  7.0E3 

The  model  input,  job  log,  and  abbreviated  output 

are  shown  subsequently,  together  with  the 

suiting  deformed  plot. 


Figure  28.  Model  used  for  linear  static  analysis. 


7.1.1  Input  Data  File  for  a  Linear  Static  Problem 

^HEADING 

MODEL  OF  A  COMPOSITE  CASING 
♦SOLUTION,  METH  =  LSA 
6,15 
♦♦ECHO 
♦♦PREPASS 
♦NODE  PRINT 
D,M 

♦PLY  FAILURE  PRINT 
♦MEMORY  ALLOCATION 
BAND 

♦GRAPHICS,  format  =  mathematica 
2,  2,  1.0 
♦  ♦ 

♦♦  NODE  DEFINITIONS 
♦  ♦ 

♦node 

1  0.0  0.0  20.0 

5  0.0  10.0  20.0 
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51 

o 

o 

0.0 

20.0 

55 

10.0 

10.0 

20.0 

151 

10.0 

0.0 

0.0 

155 

10.0 

10,0 

0.0 

201 

0.0 

0.0 

0.0 

205 

0.0 

10.0 

0.0 

1001 

0.0 

0.0 

19.875 

1005 

0.0 

10.0 

19.875 

1051 

9.875 

0.0 

19.875 

1055 

9.875 

10.0 

19.875 

1151 

9.875 

0.0 

0.125 

1155 

9.875 

10.0 

0.125 

1201 

0.0 

0.0 

0.125 

1205 

0.0 

10.0 

0.125 

2001 

0.0 

0.0 

19.75 

2005 

0.0 

10.0 

19.75 

2051 

9.75 

0.0 

19.75 

2055 

9.75 

10.0 

19.75 

2151 

9.75 

0,0 

0.25 

2155 

9.75 

10.0 

0.25 

2201 

0.0 

0.0 

0.25 

2205 

0.0 

10.0 

0.25 

3001 

0.0 

0.0 

19.625 

3005 

0.0 

10.0 

19.625 

3051 

9.625 

0.0 

19.625 

3055 

9.625 

10.0 

19.625 

3151 

9.625 

0.0 

0,375 

3155 

9.625 

10.0 

0.375 

3201 

0.0 

0.0 

0.375 

3205 

0.0 

10.0 

0.375 

4001 

0.0 

0,0 

19.5 

4005 

0.0 

10.0 

19.5 

4051 

9.5 

0.0 

19,5 

4055 

9.5 

10.0 

19.5 

4151 

9.5 

0.0 

0.5 

4155 

9.5 

10.0 

0.5 

4201 

0.0 

0.0 

0.5 

4205 

0.0 

10.0 

0.5 

>K*  GENERATE  NODE  SETS 


*MGEN,  NID=10 
1,  5,  1 
*NGEN,  NID=20 
51,  55,  1 
*NGEN,  NID=30 
151,  155,  1 
*WGEN,  NID=40 
201,  205,  1 
*NGEN,  NID=50 
1001,  1005,  1 
*NGEN,  NID=60 
1051,  1055,  1 
*NGEM,  NID=70 
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1151,  1155,  1 
*NGEN,  NID=80 
1201,  1205,  1 
*NGEN,  NID=90 
2001,  2005,  1 
*NGEN,  NID=100 
2051,  2055,  1 
*NGEN,  NID=110 
2151,  2155,  1 
♦NGEN,  NID=120 
2201,  2205,  1 
*NGEN,  NID=130 
3001,  3005,  1 
*NGEN,  NID=140 
3051,  3055,  1 
*NGEN,  NID=150 
3151,  3155,  1 
*NGEN,  NID=160 
3201,  3205,  1 
*NGEN,  NID=170 
4001,  4005,  1 
*NGEN,  NID=180 
4051,  4055,  1 
♦NGEN,  NID=190 
4151,  4155,  1 
♦NGEN,  NID=200 
4201,  4205,  1 

**  generate  node  planes 
*nfill 

10,  20,  9,  5 

20,  30,  19,  5 

30,  40,  9,  5 

50,  60,  9,  5 

60,  70,  19,  5 

70,  80,  9,  5 

90,  100,  9,  5 

100,  110,  19,  5 
110,  120,  9,  5 

130,  140,  9,  5 

140,  150,  19,  5 
150,  160,  9,  5 

170,  180,  9,  5 

180,  190,  19,  5 
190,  200,  9,  5 

** 

**  DEFINE  GENERATOR  ELEMENTS 

*+  top  horizontal  component 
** 

♦ELEMENT,  LAYUP  =  1 

1,  1001,  1006,  1007,  1002,  1, 

♦ELEMENT,  LAYUP  =  2 


100,  2001,  2006,  2007,  2002,  1001,  1006,  1007,  1002 
♦ELEMENT,  LAYUP  =  3 

200,  3001,  3006,  3007,  3002,  2001,  2006,  2007,  2002 
♦ELEMENT,  LAYUP  =  4 

300,  4001,  4006,  4007,  4002,  3001,  3006,  3007,  3002 
♦♦ 

♦*  bottom  horizontal  component 
♦♦ 

♦ELEMENT,  LAYUP  =  1 

400,  201,  196,  197,  202,  1201,  1196,  1197,  1202 

♦ELEMENT,  LAYUP  =  2 

500,  1201,  1196,  1197,  1202,  2201,  2196,  2197,  2202 
♦ELEMENT,  LAYUP  =  3 

600,  2201,  2196,  2197,  2202,  3201,  3196,  3197,  3202 
♦ELEMENT,  LAYUP  =  4 

700,  3201,  3196,  3197,  3202,  4201,  4196,  4197,  4202 
♦♦ 

♦♦  vertical  component 
♦♦ 

♦ELEMENT,  LAYUP  =  5,  orientation  =  1 
800,  1151,  151,  152,  1152,  1146,  146,  147,  1147 

♦ELEMENT,  LAYUP  =  6,  orientation  =  1 
900,  2151,  1151,  1152,  2152,  2146,  1146,  1147,  2147 
♦ELEMENT,  LAYUP  =  7,  orientation  =  1 
1000,  3151,  2151,  2152,  3152,  3146,  2146,  2147,  3147 
♦ELEMENT,  LAYUP  =  8,  orientation  =  1 
1100,  4151,  3151,  3152,  4152,  4146,  3146,  3147,  4147 
** 

**  element  generation 
** 

♦ELGEN,  ELID  =  10 

1,10,5,1,4,1,10, , , 

100.10.5.1.4.1.10.. , 

200.10.5.1.4.1.10.. , 

300.10.5.1.4.1.10.. , 

+ELGEN,  ELID  =  20 

800,0,0,0,4,1,1,20,-5,4 
900,0,0,0,4,1,1,20,-5,4 
1000,0,0,0,4,1,1,20,-5,4 
1100,0,0,0,4,1,1,20,-5,4 
*ELGEN,  ELID  =  30 

400.10, -5,5,4,1,4,,, 

500.10, -5,5,4,1,4, , , 

600.10, -5,5,4,1,4, , , 

700.10, -5,5,4,1,4,,, 

♦♦ 

♦♦  DEFINE  SET  OF  ELEMENTS  AT  UPPER  AND  LOWER  PLATE  JUNCTIONS 


+ELSET,  ESID  =  500 


10 

20 

30 

40 

109 

119 

129 

139 

209 

219 

229 

239 

309 

319 

329 

339 

445 

449 

453 

457 

545 

549 

553 

557 

645 

649 

653 

657 

745 

749 

753 

757 

800 

801 

802 

803 

900 

901 

902 

903 

1000 

1001 

1002 

1003 

1100 

1101 

1102 

1103 

876 

877 

878 

879 

976 

977 

978 

979 

1076 

1077 

1078 

1079 
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1176  1177  1178  1179 

**  EXCLUDE  ELEMENT  SET  500  FROM  FIRST  PLY  FAILURE  PREDICTION 
** 

♦EXCLUDE  ELEMENTS 
500 
** 

*♦  RELAX  TOLERANCE  IN  DEFORMED  GEOMETRY  CHECKS 
** 

♦DEFORMED  GEOMETRY 

20.0 

♦♦  LAMINATE  DESCRIPTION 


♦♦ 


♦LAYER,  LAYUP  =  1 

1,  0.0416, 

90.0 

1,  0.0208, 

45.0 

1,  0.0208, 

-45.0 

♦LAYER,  LAYUP  =  2 

1,  0.0208, 

45.0 

1,  0.0208, 

-45.0 

1,  0.1456, 

0.0 

♦LAYER,  LAYUP  =  3 

1,  0.1456, 

0.0 

1,  0.0208, 

-45.0 

1,  0.0208, 

45.0 

♦LAYER,  LAYUP  =  4 

1,  0.0208, 

-45.0 

1,  0.0208, 

45.0 

1,  0.0416, 

90.0 

♦* 

♦LAYER,  LAYUP  =  5 

1,  0.0312, 

0.0 

1,  0.0208, 

45.0 

1,  0.0208, 

-45.0 

1,  0.0260, 

90.0 

1,  0.0260, 

45.0 

1,  0.0260, 

-45.0 

♦LAYER,  LAYUP  =  6 

1,  0.0312, 

0.0 

1,  0.0260, 

45.0 

1,  0.0260, 

-45.0 

1,  0.0208, 

0,0 

♦LAYER,  LAYUP  =  7 

1,  0.0208, 

0.0 

1,  0.0260, 

-45.0 

1,  0.0260, 

45.0 

1,  0.0312, 

0.0 

♦LAYER,  LAYUP  =  8 

1,  0.0260, 

-45.0 

1,  0.0260, 

45.0 

1,  0.0260, 

90.0 

1,  0.0208, 

-45.0 

1,  0.0208, 

45.0 

45 


1,  0.0312,  0.0 

**  MATERIAL  DEFINITION 
** 

★MATERIAL,  MATID  =  1 

2.48E6,2.48E6,0.71E6,0.3E6,0.3E6,0.99E6 
0.28,0.28,0.243 
★FAILURE  CRITERIA,  FCID  =  1 
MAX-STRESS 

21.2E3,  47.8E3,  21.2E3,  47.8E3,  21.2E3,  21.0E3 
1.65E3,  1.65E3,  7.0E3 
★DAMAGE  LAW,  DLID  =  1 
NULL 
** 

★★  ESTABLISH  LOCAL  COORDINATE  SYSTEM  FOR  VERTICAL  SECTION 
*★ 

★ORIENTATION 

1,  0.0,  1.0,  0.0,  0.0,  0.0,  1.0 

★★ 

★BOUND ARY2 
10,  1,  1 
50,  1,  1 
90,  1,  1 
130,  1,  1 
170,  1,  1 
★BOUND ARY2 
40,  1,  3 
80,  1,  3 
120,  1,  3 
160,  1,  3 
200,  1,  3 

★CLOAD 

46,  3,  -50.0 

47,  3,  -100.0 

48,  3,  -100.0 

49,  3,  -100.0 

50,  3,  -50.0 

51,  3,  -50.0 

52,  3,  -100.0 

53,  3,  -100.0 

54,  3,  -100.0 

55,  3,  -50.0 

101,  1,  -100.0 
102,  1,  -200.0 

103,  1,  -200.0 

104,  1,  -200.0 

105,  1,  -100.0 

★ENDDATA 


46 


7.1.2  Timing  Summary  of  Job  Execution 


@@  @@ 

@@  U.S.  ARMY  RESEARCH  LABORATORY  @@ 

m  @0 

m  RESTRAN  m 

@@  m 

@©  RESIDUAL  STRENGTH  ANALYSIS  OF  IMPACT  m 
m  DAMAGED  COMPOSITE  LAMINATES  @@ 

m  m 

m  VERSION  1.0  m 

m  i9@ 


MODEL  SIZE  PARAMETERS 


NUMBER  OF  ELEMENTS  =  640 
NUMBER  OF  NODES  =  1025 
DEGREES  OF  FREEDOM  =  3075 
SYSTEM  BANDWIDTH  =  183 


TOTAL  RAM  MEMORY  AVAILABLE  =  14500000  WORDS 

TOTAL  RAM  MEMORY  REQUIRED  =  784080  WORDS 

:|e  3|c  :)jc  ^  :fc  :4c  9|e  3|c  :(c  9)c  3te  3|c ^  :(c  :4e  :(c  :(c  :|c  3tc 

TIMING  SUMMARY  OF  PROGRAM  EXECUTION  ** 

*:4t 

:4c  :(c  :4e  4e  *  ie  3|c  :i|c  :|c  :4c  34c  :(e  :4e  :4c  :<(  :f;  :4c  :f;  :4e  :fe  %  :f;  :tc  :4c  :4(:^  :(c  :4c  :4c  :tc  :4e  :itc  :4c  :4c  :|e  :4c  :4e  :4c  :tc  :fc  % 


RESTRAN  PREFACE  OPERATIONS:  5.3739  SECONDS 


@@@  @@@ 
<9@@  BEGINNING  LINEAR  STATIC  ANALYSIS 

@@@ 


ELASTIC  STIFFNESS  ASSEMBLY  TIME: 
BOUNDARY  CONSTRAINT  PROCESSING  TIME: 
LINEAR  SOLUTION  TIME: 

STRESS  RECOVERY  TIME: 

GRAPHICAL  OUTPUT  GENERATION  TIME: 


13.8543  SECONDS 
0.8140  SECONDS 
2.1745  SECONDS 
13.9383  SECONDS 
2.6908  SECONDS 
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TOTAL  JOB  EXECUTION  TIME:  43.2183  SECONDS 


7.1.3  Linear  Static  Solution  Output  File 


(§@  @@ 

@@  U.S.  ARMY  RESEARCH  LABORATORY  @@ 

@@  @@ 

@@  RESTRAN  @@ 

@@  @@ 

@0  RESIDUAL  STRENGTH  ANALYSIS  OF  IMPACT  @@ 
@@  DAMAGED  COMPOSITE  LAMINATES  @@ 

@@  @@ 

@@  VERSION  1.0  ®@ 

@@  @@ 


***  WARNING:  MATERIAL  ID  1  HAS  A  MIXED  MODE  FAILURE  CRITERIA 

ASSOCIATED  WITH  A  SINGLE  MODE  DAMAGE  LAW  (NULL) . 

NULL  ACCEPTED 

###################################### 

###  ### 

###  MODEL  DEFINITION### 

###  ### 

###################################### 

JOB  HEADING: 


MODEL  OF  A  COMPOSITE  CASING 

SOLUTION  CONTROL: 

SOLUTION  ALGORITHM: 

LINEAR  STATIC  ANALYSIS 
MEMORY  ALLOCATION: 

BANDWIDTH  STORAGE  MODE 

*  BANDWIDTH  MINIMIZATION  WILL  BE  PERFORMED 
PROGRAM  OPTIONS: 
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*  EXTENDED  JOB  EXECUTION  STATISTICS  WILL  BE  PRINTED 


*  THE  FOLLOWING  ELEMENT  SETS  ARE  EXCLUDED  FROM  FAILURE  PREDICTION 
500 


REQUESTED  PRINT  OPTIONS: 


*  PRINT  DISPLACEMENTS  AFTER  CONVERGENCE  AT  THE  iTH  LOAD  INCREMENT 

*  PRINT  BUCKLING  MODE  SHAPE  CORRESPONDING  TO  THE  LOWEST  BUCKLING 
LOAD  MULTIPLIER  AT  THE  iTH  ITERATION 

*  PRINT  ELEMENT  PLY  FAILURE  STATISTICS 

GRAPHICS: 

@  GRAPHICAL  OUTPUT  IN  MATHEMATICA  FORMAT  WILL  BE  GENERATED 

@  BOTH  UNDEFORMED  AND  DEFORMED  GEOMETRY  IS  REQUESTED 
DEFORMATION/MODE  MAGNIFICATION  FACTOR  =  O.lOOOE+01 

@  RESTRAN  GRAPHICS  CONTAINED  IN  FILE:  restran.grf 

NODE  DEFINITIONS: 


DDE  ID 

X 

Y 

Z 

1 

O.OOOE+00 

O.OOOE+00 

0.200E+02 

2 

O.OOOE+00 

0.250E+01 

0.200E+02 

3 

O.OOOE+00 

0.500E+01 

0.200E+02 

4 

O.OOOE+00 

0.750E+01 

0.200E+02 

5 

O.OOOE+00 

O.lOOE+02 

0.200E+02 

4201 

O.OOOE+00 

O.OOOE+00 

0.500E+00 

4202 

O.OOOE+00 

0.250E+01 

0.500E+00 

4203 

O.OOOE+00 

0.500E+01 

0.500E+00 

4204 

O.OOOE+00 

0.750E+01 

0.500E+00 

4205 

O.OOOE+00 

O.lOOE+02 

0.500E+00 

S  P 

E  C  I  F  I  E 

D  NODE 

SETS: 

SET  ID  10: 
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1 


2 


3 


4 


5 


SET  ID  20 : 


51 

52 

53 

54 

55 

SET  ID 

190: 

4151 

4152 

4153 

4154 

4155 

SET  ID 

200: 

4201 

4202 

4203 

4204 

4205 

E  L  E 

M  E  N  T 

D  E  F  I 

N  I  T  I 

0  N  S: 

ELEMENT 

N1 

N2 

N3 

N4 

N5 

N6 

N7 

N8 

1 

1001 

1006 

1007 

1002 

1 

6 

7 

2 

2 

1006 

1011 

1012 

1007 

6 

11 

12 

7 

3 

1011 

1016 

1017 

1012 

11 

16 

17 

12 

4 

1016 

1021 

1022 

1017 

16 

21 

22 

17 

5 

1021 

1026 

1027 

1022 

21 

26 

27 

22 

1175 

4064 

3064 

3065 

4065 

4059 

3059 

3060 

4060 

1176 

4056 

3056 

3057 

4057 

4051 

3051 

3052 

4052 

1177 

4057 

3057 

3058 

4058 

4052 

3052 

3053 

4053 

1178 

4058 

3058 

3059 

4059 

4053 

3053 

3054 

4054 

1179 

4059 

3059 

3060 

4060 

4054 

3054 

3055 

4055 

E  L  E  M  E  N 

T  S  E 

T  S: 

ELEMENT  SET  ID 

10: 

2 

3 

4 

5 

ELEMENT  SET  ID 

20: 

801 

802 

803 

804 

ELEMENT  SET  ID 

30: 

405 

410 

415 

420 

... 
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ELEMENT  SET  ID  500: 


10  20  30  40  . . . 

MATERIAL  PROPERTY  DATA: 


MATERIAL  ID  1 

El  =  0.248E+07  G23  =  0.300E+06  MU23  =  0.28000 

E2  =  0.248E+07  G13  =  0.300E+06  MU13  =  0.28000 

E3  =  0.710E+06  G12  =  0.990E+06  MU12  =  0.24300 

ASSOCIATED  FAILURE  LAW: 

MAXIMUM  STRESS  CRITERION 

XT  =  0.212E+05  XC  =  0.478E+05  YT  =  0.212E+05 
YC  =  0.478E+05  ZT  =  0.212E+05  ZC  =  0.210E+05 
R  =  0.165E+04  S  =  0.165E+04  T  =  0.700E+04 

ASSOCIATED  DAMAGE  LAW: 

NULL  POST-FAILURE  PROPERTIES 


MATERIAL  LAYER  DEFINITIONS 


LAYUP  ID  =  1 

NUMBER  OF  PLIES  =  3 

MATERIAL  ID  PLY  THICKNESS  ORIENTATION 

1  0.04160  90.00000 

1  0 . 02080  45 . 00000 

1  0.02080  -45.00000 

LAYUP  ID  =  2 

NUMBER  OF  PLIES  =  3 

MATERIAL  ID  PLY  THICKNESS  ORIENTATION 

1  0.02080  45.00000 

1  0.02080  -45.00000 

1  0 . 14560  0 . 00000 

LAYUP  ID  =  3 

NUMBER  OF  PLIES  =  3 

MATERIAL  ID  PLY  THICKNESS  ORIENTATION 

1  0.14560  0.00000 

1  0.02080  -45.00000 
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1 


0.02080 


45.00000 


LAYUP  ID  =  4 

NUMBER  OF  PLIES  =  3 

MATERIAL  ID  PLY  THICKNESS  ORIENTATION 


1 

0.02080 

-45.00000 

1 

0.02080 

45.00000 

1 

0.04160 

90.00000 

LAYUP  ID  = 

5 

NUMBER  OF  PLIES  =  6 

MATERIAL  ID 

PLY  THICKNESS 

ORIENTATION 

1 

0.03120 

0.00000 

1 

0.02080 

45.00000 

1 

0.02080 

-45.00000 

1 

0.02600 

90.00000 

1 

0 . 02600 

45.00000 

1 

0.02600 

-45.00000 

LAYUP  ID  = 

6 

NUMBER  OF  PLIES  =  4 

MATERIAL  ID 

PLY  THICKNESS 

ORIENTATION 

1 

0.03120 

0.00000 

1 

0.02600 

45.00000 

1 

0.02600 

-45.00000 

1 

0.02080 

0.00000 

LAYUP  ID  = 

7 

NUMBER  OF  PLIES  =  4 

MATERIAL  ID 

PLY  THICKNESS 

ORIENTATION 

1 

0.02080 

0.00000 

1 

0.02600 

-45.00000 

1 

0.02600 

45.00000 

1 

0.03120 

0.00000 

LAYUP  ID  = 

8 

NUMBER  OF  PLIES  =  6 

MATERIAL  ID 

PLY  THICKNESS 

ORIENTATION 

1 

0.02600 

-45.00000 

1 

0.02600 

45.00000 

1 

0,02600 

90.00000 

1 

0.02080 

-45.00000 

1 

0.02080 

45.00000 

1 

0.03120 

0.00000 

MATERIAL  COORDINATE  SYSTEMS 


ID  Vll  V12  V13  V21  V22 

1  O.OOOE+00  O.lOOE+01  O.OOOE+00  O.OOOE+00  O.OOOE+00 


BOUNDARY  CONDITIONS: 


SINGLE  POINT  CONSTRAINTS: 

NODE  DOFl  TO  D0F2 

1  11 

2  11 

3  1  1 

4  11 


4202 

1 

3 

4203 

1 

3 

4204 

1 

3 

4205 

1 

3 

APPLIED  CONCENTRATED  FORCES: 

NODE 

DIRECTION 

MAGNITUDE 

46 

3 

-.500E+02 

47 

3 

-.lOOE+03 

48 

3 

- . lOOE+03 

49 

3 

-.lOOE+03 

50 

3 

-.500E+02 

51 

3 

-.500E+02 

52 

3 

-.lOOE+03 

53 

3 

- . lOOE+03 

54 

3 

- . lOOE+03 

55 

3 

-.500E+02 

101 

1 

-.lOOE+03 

102 

1 

-.200E+03 

103 

1 

-.200E+03 

104 

1 

- . 200E+03 

105 

1 

- . lOOE+03 

FORCE/MOMENT  RESULTANTS  AT  ORIGIN: 
IRI  =  0.113E+04  |M|  =  0.567E+04 


Rx  =  -.800E+03  Ry  =  O.OOOE+00  Rz  =  -.800E+03 


V23 

O.lOOE+01 
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Mx  =  0.400E+04  My  =  0.400E+03  Mz  =  -.400E+04 
Xo  =  0.950E+01  Yo  =  0.500E+01  Zo  =  O.lOOE+02 


***END  OF  MODEL 
MODEL  SIZE 


DEFINITION*** 

PARAMETERS 


NUMBER  OF  ELEMENTS  =  640 
NUMBER  OF  NODES  =  1025 
DEGREES  OF  FREEDOM  =  3075 
SYSTEM  BANDWIDTH  =  183 


############################################### 
##  ## 

##  LINEAR  STATIC  ANALYSIS## 
##  ## 

############################################### 


DISPLACEMENT  FIELD 
(SCALE  FACTOR  =  O.lOOE+01) 


NODE 

U 

V 

w 

1 

O.OOOOOE+00 

0.27644E-03 

-.24278E+00 

2 

O.OOOOOE+00 

0.60853E-04 

-.24592E+00 

3 

O.OOOOOE+00 

" . 12759E-03 

-.24693E+00 

4 

O.OOOOOE+00 

-.31570E-03 

- . 24594E+00 

5 

O.OOOOOE+00 

-.53044E-03 

-.24281E+00 

4201 

0 . OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 

4202 

O.OOOOOE+00 

O.OOOOOE+00 

0 . OOOOOE+00 

4203 

O.OOOOOE+00 

O.OOOOOE+00 

0 . OOOOOE+00 

4204 

0 . OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 

4205 

0 . OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 

MAXIMUM  ELEMENT  FAILURE  INDICES 


ELEMENT 

PLY 

FAILURE  INDEX 

MODE 

MULTIPLIER 

1 

3 

0.614E-01 

GENERAL 

0.163E+02 

2 

3 

0,613E-01 

GENERAL 

0.163E+02 

3 

3 

0.610E-01 

GENERAL 

0.164E+02 

4 

3 

0.610E-01 

GENERAL 

0.164E+02 

5 

3 

0.603E-01 

GENERAL 

0.166E+02 

1174 

5 

0.266E-01 

GENERAL 

0.377E+02 

1175 

5 

0.322E-01 

GENERAL 

0.311E+02 

1176 

1 

0.691E-01 

GENERAL 

0 . 145E+02 
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1177 

1 

0.633E-01 

GENERAL 

0.158E+02 

1178 

1 

0.633E-01 

GENERAL 

0.158E+02 

1179 

1 

0.692E-01 

GENERAL 

0 . 144E+02 

HIGHEST  GLOBAL  FAILURE  INDEX  OCCURRED  IN  ELEMENT  902 
IN  PLY  NUMBER  1  PREDICTING  FIRST  PLY  FAILURE  AT  A 
FACTOR  TO  APPLIED  LOADS  EQUAL  TO  0.174E+01 


@®  ®@ 

@®  STATIC  ANALYSIS  COMPLETED  ®@ 
@®  ®® 


@®®@@®®®@@®®@@@@®@®@®@@@®®®®@®®®®®®®®®®@®®@®®®®®®@®®®®®® 


According  to  a  first  ply  failure  criterion,  failure  is  predicted  at  a  load  given  by 

{P}  =  1.74{Papplied}  (19) 

A  graphical  depiction  of  the  deformed  shape  is  shown  in  Figure  29. 


Figure  29.  Deformation  under  applied  static  loads. 


7.2  Linear  Buckling  Analysis 

A  similar  model  of  a  composite  casing  as  shown  in  Section  7.1  is  analyzed  for  buckling.  The  same 
laminate  layup  and  material  properties  were  used.  The  model,  shown  in  Figure  30,  was  constructed 
to  model  circular  delamination  at  each  of  the  three  interior  sublaminate  interfaces  in  the  vertical 
plate.  The  location  of  these  delaminations  are  indicated  by  ||  in  the  layup  description  given  by 
(Oe/ ±454/905/ ±455  II  06/±455/04  ||  O4/T455/O6  ||  ^465/^05/ The  instability  of  these 
delaminations  was  analyzed  simultaneously.  For  comparison,  the  critical  buckling  load  of  the  entire 
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component  without  the  presence  of  delaminations  was  computed  by  removing  the  *DELAMINA 
TION  statement  from  the  input  file.  The  model  input,  job  log,  and  abbreviated  output  files  are 
shown  subsequently,  together  with  the  resulting  mode  shape  plots. 


Figure  30.  Model  used  for  linear  buckling  analysis. 


7.2.1  Input  Data  File  for  a  Linear  Buckling  Problem 

^HEADING 

STABILITY  ANALYSIS  OF  A  DELAMINATED  COMPOSITE  CASING 
^SOLUTION,  METH  =  LBA 
6,15 
**ECH0 
**PREPASS 
*N0DE  PRINT 
D,M 

*MEM0RY  ALLOCATION 
BAND 

^GRAPHICS,  format  =  mathematica 
2,  2,  1.0 
** 

**  NODE  DEFINITIONS 
*N0DE 

**  TOP  HORIZONTAL  COMPONENT,  LEFT 


1 

0.0 

-0.1000005E+02 

10.0 

2 

0.0 

-0.9009584E+01 

10.0 

3 

0.0 

-0.80191i9Et01 

10.0 
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TOP 

** 

HORIZONTAL  COMPONENT, 

RIGHT 

151 

10.0 

-0.1000005E+02 

10.0 

152 

10.0 

-0.9009584E+01 

10.0 

153 

10.0 

-0.8019119E+01 

10.0 

**  BOTTOM  HORIZONTAL  COMPONENT,  LEFT 


** 

1501 

0.0 

-0.1000005E+02 

-lO.O 

1502 

0.0 

-0.9009584E+01 

-10.0 

1503 

0.0 

-0.8019119E+01 

-10.0 

**  BOTTOM  HORIZONTAL  COMPONENT,  RIGHT 


1651 

10.0 

-0.1000005E+02 

-10.0 

1652 

10.0 

-0.9009584E+01 

-10.0 

1653 

10.0 

-0.8019119E+01 

-10.0 

**  DEFINE  NODE  SETS  FOR  FILL  OPERATION 


*NSET, 

NSID 

=  10 

1, 

2, 

3, 

4, 

5, 

6, 

7, 

8, 

9, 

10 

11, 

12, 

13, 

14, 

15 

♦NSET, 

NSID 

=  20 

151, 

152, 

153, 

154, 

155, 

156, 

157, 

158, 

159, 

160 

161, 

162, 

163, 

164, 

165 

♦NSET, 

NSID 

=  30 

301, 

302, 

303, 

304, 

305, 

306, 

307, 

308, 

309, 

310 

311, 

312, 

313, 

314, 

315 

*iiset , 

ns  id 

=  620 

7098, 

7113, 

7128, 

7143, 

7158, 

7166, 

7174, 

7182,  7190, 

7198 

7206, 

7214, 

7229, 

7244, 

7259 

★nset , 

nsid 

=  630 

8098, 

8113, 

8128, 

8143, 

8158, 

8166, 

8174, 

8182,  8190, 

8198 

8206, 

8214, 

8229, 

8244, 

8259 

*nset , 

nsid 

=  640 

9098, 

9113, 

9128, 

9143, 

9158, 

9166, 

9174, 

9182,  9190, 

9198 

9206, 

9214, 

9229, 

9244, 

9259 
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**  EQUIVALENCE  NODES  TO  JOIN  SEPARATE  MODEL  DOMAINS 

♦EQUIVALENCE 
500,  20 

510,  40 

520,  60 

530,  80 

540,  100 
600,  120 
610,  140 
620,  160 
630,  180 
640,  200 
** 

♦*  GENERATE  NODE  PLANES 


♦NFILL 


10, 

20, 

9, 

15 

30, 

40, 

9, 

15 

50, 

60, 

9, 

15 

70, 

80, 

9. 

15 

90, 

100, 

9, 

15 

110, 

120, 

9, 

15 

130, 

140, 

9, 

15 

150, 

160, 

9, 

15 

170, 

180, 

9, 

15 

190, 

200, 

9, 

15 

♦♦  NODES,  ELEMENTS,  AND  NODE  SETS  FOR  DELAMINATED  VERTICAL 
♦♦  COMPONENT  OBTAINED  FROM  THE  FOLLOWING  ♦MODEL  GENERATION 
♦♦  STATEMENT: 


♦MODEL  GENERATION 

4,5,10.0,10.0,5,5,5,5 

1.5, 1,1, 1,1, 1,1 

♦♦ 

0.0,  0.0,  1.0 

0.0,  1.0,  0.0 

** 

1.0,  0.0,  0.0 

b 

o 

o 

** 

1,6,250,5000,5000 

** 

1.0,9.875 

** 

2,6,255,6000,6000 

** 

1.0,9.750 

** 

3,6,260,7000,7000 

** 

1.0,9.625 

** 

4,6,265,8000,8000 

1.0,9.50 

** 

5,6,270,9000,9000 

** 

♦NODE 

5001 

0 . lOOOOOOE+02  0 . OOOOOOOE+00  0 . OOOOOOOE+00 

5002 

0 . 1000000E+02  0 . 1333340E+01  0 . OOOOOOOE+00 

5003 

0 . lOOOOOOE+02  0 . 1333340E+01  0 . 1333340E+01 

5004 

0 . lOOOOOOE+02  0 . OOOOOOOE+00  0 . 1333340E+01 
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5005 

O.lOOOOOOE+02 

-0.1333340E+01 

0.1333340E+01 

9268 

0.9500000E+01 

0 . 1000005E+02 

0. 489451 lE+01 

9269 

0.9500000E+01 

0.1000005E+02 

0.6038188E+01 

9270 

0.9500000E+01 

0.1000005E+02 

0.7028653E+01 

9271 

0.9500000E+01 

0.1000005E+02 

0.8019119E+01 

9272 

0.9500000E+01 

0.1000005E+02 

0.9009584E+01 

9273 

0.9500000E+01 

0.1000005E+02 

0.1000005E+02 

*♦  NODE  SETS  TO  DEFINE  DELAMINATION  PLANES 


♦NSET,  NSID  =  255 


6001 

6002 

6003 

6004 

6005 

6006 

6007 

6008 

6009 

6010 

6011 

6012 

6013 

6014 

6015 

6016 

6017 

6018 

6019 

6020 

6021 

6022 

6023 

6024 

6025 

6026 

6027 

6028 

6029 

6030 

6031 

6032 

6033 

6034 

6035 

6036 

6037 

6038 

6039 

6040 

6041 

*NSET,  NSID  = 

260 

7001 

7002 

7003 

7004 

7005 

7006 

7007 

7008 

7009 

7010 

7011 

7012 

7013 

7014 

7015 

7016 

7017 

7018 

7019 

7020 

7021 

7022 

7023 

7024 

7025 

7026 

7027 

7028 

7029 

7030 

7031 

7032 

7033 

7034 

7035 

7036 

7037 

7038 

7039 

7040 

7041 

♦NSET,  NSID  = 

265 

8001 

8002 

8003 

8004 

8005 

8006 

8007 

8008 

8009 

8010 

8011 

8012 

8013 

8014 

8015 

8016 

8017 

8018 

8019 

8020 

8021 

8022 

8023 

8024 

8025 

8026 

8027 

8028 

8029 

8030 

8031 

8032 

8033 

8034 

8035 

8036 

8037 

8038 

8039 

8040 

8041 

**  DEFINE  GENERATOR  ELEMENTS 


**  top  horizontal  component 


♦ELEMENT,  LAYUP  = 
1,  301,  316, 

1 

317, 

302, 

1, 

16, 

17, 

2 

♦ELEMENT,  LAYUP  = 
300,  601,  616, 

2 

617, 

602, 

301, 

316, 

317, 

302 

♦ELEMENT,  LAYUP  = 
600,  901,  916, 

3 

917, 

902, 

601, 

616, 

617, 

602 

♦ELEMENT,  LAYUP  = 
900,  1201,  1216, 

4 

1217, 

1202, 

901, 

916, 

917, 

902 

** 

**  bottom  horizontal  component 

** 

♦ELEMENT,  LAYUP  =  1 

1200,  1801,  1816,  1817,  1802, 

1501, 

1516, 

1517, 

1502 

♦ELEMENT,  LAYUP  = 
1500,  2101,  2116, 

2 

2117, 

2102, 

1801, 

1816, 

1817, 

1802 
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★ELEMENT,  LAYUP  =  3 

1800,  2401,  2416,  2417,  2402,  2101,  2116,  2117,  2102 
★ELEMENT,  LAYUP  =  4 

2100,  2701,  2716,  2717,  2702,  2401,  2416,  2417,  2402 
*★ 

★★  ELEMENT  GENERATION 


★ELGEN,  ELID  =10 

1,  10,15,1,  14,1,10,,, 

300,  10,15,1,  14,1,10,,, 

600,  10,15,1,  14,1,10,,, 

900,  10,15,1,  14,1,10,,, 

1200,  10,15,1,  14,1,10,,, 

1500,  10,15,1,  14,1,10,,, 

1800,  10,15,1,  14,1,10,,, 

2100,  10,15,1,  14,1,10,,, 

★★ 

★★  ELEMENTS  IN  VERTICAL  COMPONENT  FROM  PRIOR  *MODEL  GENERATION  RUN 
★★ 

★ELEMENT,  LAYUP  =  5,  ORIENTATION  =  1 


5001 

6001 

6004 

6003 

6002 

5001 

5004 

5003 

5002 

5002 

6006 

6005 

6004 

6001 

5006 

5005 

5004 

5001 

5003 

6007 

6006 

6001 

6008 

5007 

5006 

5001 

5008 

5004 

6008 

6001 

6002 

6009 

5008 

5001 

5002 

5009 

5005 

6002 

6003 

6011 

6010 

5002 

5003 

5011 

5010 

8239 

9252 

9253 

9268 

9267 

8252 

8253 

8268 

8267 

8240 

9253 

9254 

9269 

9268 

8253 

8254 

8269 

8268 

8241 

9254 

9255 

9270 

9269 

8254 

8255 

8270 

8269 

8242 

9255 

9256 

9271 

9270 

8255 

8256 

8271 

8270 

8243 

9256 

9257 

9272 

9271 

8256 

8257 

8272 

8271 

8244 

9257 

9258 

9273 

9272 

8257 

8258 

8273 

8272 

★★ 

★★  RELAX  TOLERANCE  OF  DEFORMED  GEOMETRY  CHECKS 

★DEFORMED  GEOMETRY 
20.0 


★★  LAMINATE  DESCRIPTION 
*★ 

★LAYER,  LAYUP  =  1 
1,  0.0416,  90.0 

1,  0.0208,  45.0 

1,  0.0208,  “45.0 
★LAYER,  LAYUP  =  2 
1,  0.0208,  45.0 

1,  0.0208,  -45.0 
1,  0.1456,  0.0 

★LAYER,  LAYUP  =  3 
1,  0.1456,  0.0 

1,  0.0208,  “45.0 
1,  0.0208,  45.0 
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★LAYER,  LAYUP  =  4 


1, 

0.0208, 

-45.0 

1. 

0.0208, 

45.0 

1, 

0.0416, 

90.0 

** 

★LAYER,  LAYUP  =  5 

1. 

0.0312, 

0.0 

1, 

0.0208, 

45.0 

1, 

0.0208, 

-45.0 

1, 

0.0260, 

90.0 

1, 

0.0260, 

45.0 

1, 

0.0260, 

-45.0 

♦LAYER,  LAYUP  =  6 

1, 

0.0312, 

0.0 

1. 

0.0260, 

45.0 

1, 

0.0260, 

-45.0 

1, 

0.0208, 

0.0 

♦LAYER,  LAYUP  =  7 

1, 

0.0208, 

0.0 

1, 

0.0260, 

-45.0 

1, 

0.0260, 

45.0 

1, 

0.0312, 

0.0 

★LAYER,  LAYUP  =  8 

1, 

0.0260, 

-45.0 

1. 

0.0260, 

45.0 

1, 

0.0260, 

90.0 

1, 

0.0208, 

-45.0 

1, 

0.0208, 

45.0 

1, 

0.0312, 

0.0 

★★  MATERIAL  DEFINITION 
★★ 

★MATERIAL,  MATID  =  1 

2.48E6,2.48E6,0.71E6,0.3E6,0.3E6,0.99E6 
0.28,0.28,0.243 
★FAILURE  CRITERIA,  FCID  =  1 
MAX-STRESS 

21.2E3,  47.8E3,  21.2E3,  47.8E3,  21.2E3,  21.0E3 
1.65E3,  1.65E3,  7.0E3 
★DAMAGE  LAW,  DLID  =  1 
NULL 

★★  ESTABLISH  LOCAL  COORDINATE  SYSTEM 
★ORIENTATION 

1,  0.0,  1.0,  0.0,  0.0,  0.0,  1.0 

★B0UNDARY2 
10,  1,  1 
30,  1,  1 
50,  1,  1 
70,  1,  1 
90,  1,  1 
★B0UNDARY2 
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110,  1,  3 
130,  1,  3 
150,  1,  3 
170,  1,  3 
190,  1,  3 


+CL0AD 

136,  3, 

-50.0 

137,  3, 

-100.0 

138,  3, 

-100.0 

139,  3, 

-100.0 

140,  3, 

-100.0 

141,  3, 

-100.0 

142,  3, 

-100.0 

143,  3, 

-100.0 

144,  3, 

-100.0 

145,  3, 

-100.0 

146,  3, 

-100.0 

147,  3, 

-100.0 

148,  3, 

-100.0 

149,  3, 

-100.0 

150,  3, 

-50.0 

151,  3, 

-50.0 

152,  3, 

-100.0 

153,  3, 

-100.0 

154,  3, 

-100.0 

155,  3, 

-100.0 

156,  3, 

-100.0 

157,  3, 

-100.0 

158,  3, 

-100.0 

159,  3, 

“100.0 

160,  3, 

-100.0 

161,  3, 

-100.0 

162,  3, 

-100.0 

163,  3, 

-100.0 

164,  3, 

-100.0 

165,  3, 

-50.0 

♦ENDDTA 

7.2.2  Solution  Timing  Summary  for  Linear  Buckling  Analysis 


@@  @@ 

®@  U.S.  ARMY  RESEARCH  LABORATORY  @@ 

@@  @@ 

@@  RESTRAN  m 

@®  @@ 

@@  RESIDUAL  STRENGTH  ANALYSIS  OF  IMPACT  @@ 
m  DAMAGED  COMPOSITE  LAMINATES  ®@ 

m  @@ 

@@  VERSION  1.0  ®@ 

@@ 
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MODEL  SIZE  PARAMETERS 


NUMBER  OF  ELEMENTS  =  2096 
NUMBER  OF  NODES  =  2865 
DEGREES  OF  FREEDOM  =  8595 
SYSTEM  BANDWIDTH  =  1005 


TOTAL  RAM  MEMORY  AVAILABLE  =  14500000  WORDS 
TOTAL  RAM  MEMORY  REQUIRED  =  9232360  WORDS 


H!3tcitt:llf^:i):^st:^iitiH:itiHc*st:*********************'***** 

♦*  TIMING  SUMMARY  OF  PROGRAM  EXECUTION  ** 

:|c  3jc  ***  3fc  *  3*c  *  :ie  :te  ♦  3te  :|e  :<c  Jfc  ****  :*e  3(£  :tc  sjc  ***♦♦♦**  +  **  3tc 


RESTRAN  PREFACE  OPERATIONS:  24.5313  SECONDS 


@@(9  @@@ 

BEGINNING  LINEAR  BUCKLING  ANALYSIS  @@@ 
@@@  @@@ 


ELASTIC  STIFFNESS  ASSEMBLY  TIME; 
BOUNDARY  CONSTRAINT  PROCESSING  TIME: 
LINEAR  SOLUTION  TIME: 

STRESS  RECOVERY  TIME: 

DIFFERENTIAL  STIFFNESS  ASSEMBLY  TIME: 
BOUNDARY  CONSTRAINT  PROCESSING  TIME: 
BOUNDARY  CONSTRAINT  PROCESSING  TIME: 
EIGENVALUE  EXTRACTION  TIME: 

GRAPHICAL  OUTPUT  GENERATION  TIME: 


42.2373  SECONDS 
11.7441  SECONDS 

2  MINUTES  33.4424  SECONDS 
37.7783  SECONDS 
14.0586  SECONDS 
11.5098  SECONDS 
11.9873  SECONDS 

1  HOURS  53  MINUTES  45.6260  SECONDS 
12.9336  SECONDS 


TOTAL  JOB  EXECUTION  TIME:  1  HOURS  58  MINUTES  44.2803  SECONDS 


7*2,3  Linear  Buckling  Solution  Output  File 
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@@  U.S.  ARMY  RESEARCH  LABORATORY  @@ 
@0  @@ 
@@  RESTRAN  @@ 
@0  @@ 
@@  RESIDUAL  STRENGTH  ANALYSIS  OF  IMPACT  @@ 
@0  DAMAGED  COMPOSITE  LAMINATES  @@ 
@0  @@ 
00  VERSION  1.0  @@ 
@0  @@ 


00@@00@@000@@@00@000@000000000@@@@@@@@@@@®@®@@ 

***  MESSAGE:  ELEMENT  NODE  ORDER  IS  BEING  CONVERTED  TO  RESTRAN  FORMAT 

###################################### 

###  ### 

###  MODEL  DEFINITION### 

###  ### 

###################################### 

JOB  HEADING: 


MODEL  OF  A  COMPOSITE  CASING 


SOLUTION  CONTROL: 


SOLUTION  ALGORITHM: 
LINEAR  BUCKLING  ANALYSIS 


MEMORY  ALLOCATION: 
BANDWIDTH  STORAGE  MODE 


*  BANDWIDTH  MINIMIZATION  WILL  BE  PERFORMED 
PROGRAM  OPTIONS: 


♦  EXTENDED  JOB  EXECUTION  STATISTICS  WILL  BE  PRINTED 


REQUESTED  PRINT  OPTIONS: 


*  PRINT  DISPLACEMENTS  AFTER  CONVERGENCE  AT  THE  iTH  LOAD  INCREMENT 


*  PRINT  BUCKLING  MODE  SHAPE  CORRESPONDING  TO  THE  LOWEST  BUCKLING 
LOAD  MULTIPLIER  AT  THE  iTH  ITERATION 


GRAPHICS: 


@  GRAPHICAL  OUTPUT  IN  MATHEMATICA  FORMAT  WILL  BE  GENERATED 

0  BOTH  UNDEFORMED  AND  DEFORMED  GEOMETRY  IS  REQUESTED 
DEFORMATION/MODE  MAGNIFICATION  FACTOR  =  0 . lOOOE+01 

e  RESTRAN  GRAPHICS  CONTAINED  IN  FILE:  restran.grf 

NODE  DEFINITIONS: 


NODE  ID 

X 

Y 

Z 

1 

0,000E+00 

-.lOOE+02 

O.lOOE+02 

2 

O.OOOE+00 

-.901E+01 

O.lOOE+02 

3 

O.OOOE+00 

-.802E+01 

O.lOOE+02 

4 

O.OOOE+00 

-.703E+01 

O.lOOE+02 

5 

O.OOOE+00 

-.604E+01 

O.lOOE+02 

9268 

0.950E+01 

0 . lOOE+02 

0.489E+01 

9269 

0.950E+01 

0 . lOOE+02 

0.604E+01 

9270 

0.950E+01 

0 . lOOE+02 

0.703E+01 

9271 

0.950E+01 

O.lOOE+02 

0.802E+01 

9272 

0.950E+01 

O.lOOE+02 

0.901E+01 

S  P 

E  C  I  F  I  E 

D  NODE 

SETS: 

SET  ID  10: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

SET  ID 

265: 

8001 

8011 

8002 

8012 

8003 

8013 

8004 

8014 

8005 

8015 

8006 

8016 

8007 

8017 

8008 

8018 

8009 

8019 

8010 

8020 
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8021 

8031 

8041 

8022 

8032 

8023  8024 

8033  8034 

8025  8026  8027 

8035  8036  8037 

8028 

8038 

8029 

8039 

8030 

8040 

ELEMENT 

DEFINITIONS: 

ELEMENT 

N1 

N2 

N3 

N4 

N5 

N6 

N7 

N8 

1 

301 

316 

317 

302 

1 

16 

17 

2 

2 

316 

331 

332 

317 

16 

31 

32 

17 

3 

331 

346 

347 

332 

31 

46 

47 

32 

4 

346 

361 

362 

347 

46 

61 

62 

47 

5 

361 

376 

377 

362 

61 

76 

77 

62 

8240 

9253 

9268 

9269 

9254 

8253 

8268 

8269 

8254 

8241 

9254 

9269 

9270 

9255 

8254 

8269 

8270 

8255 

8242 

9255 

9270 

9271 

9256 

8255 

8270 

8271 

8256 

8243 

9256 

9271 

9272 

9257 

8256 

8271 

8272 

8257 

8244 

9257 

9272 

1365 

1364 

8257 

8272 

1065 

1064 

E  L 

EMENT  SET 

ELEMENT  SET  ID 

S: 

10: 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

2218 

2219 

2220 

2221 

2222 

2223 

2224 

2225 

2226 

2227 

2228 

2238 

2229 

2239 

2230 

2231 

2232 

2233 

2234 

2235 

2236 

2237 

M 

A  T  E  R  I 

A  L 

P 

ROPER 

T  Y  D 

A  T  A: 

MATERIAL  ID  1 

El 

= 

0 . 248E+07 

G23 

0.300E+06 

MU23  = 

0.28000 

E2 

= 

0.248E+07 

G13 

=: 

0.300E+06 

MU13  = 

0.28000 

E3 

= 

0.710E+06 

G12 

= 

0.990E+06 

MU12  = 

0.24300 

ASSOCIATED  FAILURE  LAW: 

MAXIMUM  STRESS  CRITERION 

XT  =  0.212E+05  XC  =  0.478E+05  YT  =  0.212E+05 
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YC  =  0.478E+05  ZT  =  0.212E+05  ZC  =  0.210E+05 
R  =  0.165E+04  S  =  0.165E+04  T  =  0.700E+04 

ASSOCIATED  DAMAGE  LAW: 

NULL  POST-FAILURE  PROPERTIES 


MATERIAL  ID  2 

El  =  O.lOOE+06  G23  =  0.385E+05  MU23  =  0.30000 
E2  =  O.lOOE+06  G13  =  0.385E+05  MU13  =  0.30000 
E3  =  O.lOOE+06  G12  =  0.385E+05  MU12  =  0.30000 

ASSOCIATED  FAILURE  LAW: 

MAXIMUM  STRESS  CRITERION 

XT  =  0.212E+05  XC  =  0.478E+05  YT  =  0.212E+05 
YC  =  0.478E+05  ZT  =  0.212E+05  ZC  =  0.210E+05 
R  =  0.165E+04  S  =  0.165E+04  T  =  0.700E+04 

ASSOCIATED  DAMAGE  LAW: 

NULL  POST-FAILURE  PROPERTIES 


MATERIAL  LAYER  DEFINITIONS 


LAYUP  ID  =  1 

NUMBER  OF  PLIES  =  3 


MATERIAL  ID 

1 

1 

1 


PLY  THICKNESS 

0.04160 
0 . 02080 
0.02080 


ORIENTATION 

90.00000 

45.00000 

-45.00000 


LAYUP  ID  =  88 

NUMBER  OF  PLIES  =  6 


MATERIAL  ID 

2 

2 

2 

2 

2 

2 


PLY  THICKNESS 

0.02600 

0.02600 

0.02600 

0.02080 

0.02080 

0.03120 


ORIENTATION 

-45.00000 

45.00000 

90.00000 

-45.00000 

45.00000 

0.00000 
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MATERIAL  COORDINATE  SYSTEMS 


ID  Vll  V12  V13  V21  V22 

1  O.OOOE+00  O.lOOE+01  O.OOOE+00  O.OOOE+00  O.OOOE+00 


BOUNDARY  CONDITIONS: 


SINGLE  POINT  CONSTRAINTS: 
NODE  DOFl  TO  D0F2 

1  13 

2  1  3 

3  13 

4  13 

5  1  3 


2711  1  3 

2712  1  3 

2713  1  3 

2714  1  3 

2715  1  3 


APPLIED  CONCENTRATED  FORCES: 


NODE 

DIRECTION 

MAGNITUDE 

1651 

2 

0.500E+01 

5099 

2 

O.lOOE+02 

5100 

2 

O.lOOE+02 

5101 

2 

O.lOOE+02 

9269 

2 

“ . lOOE+02 

9270 

2 

-.lOOE+02 

9271 

2 

-.lOOE+02 

9272 

2 

- . lOOE+02 

1365 

2 

■“.500E+01 

FORCE/MOMENT  RESULTANTS  AT  ORIGIN: 

|R|  =  O.OOOE+00  |M|  =  O.OOOE+00 
Rx  =  O.OOOE+00  Ry  =  O.OOOE+00  Rz  =  O.OOOE+OO 


V23 
O.lOOE+01 
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Mx  =  O.OOOE+00  My  =  O.OOOE+00  Mz  =  O.OOOE+00 
Xo  =  O.OOOE+00  Yo  =  O.OOOE+00  Zo  =  O.OOOE+00 


♦+*END  OF  MODEL  DEFINITION*** 


MODEL  SIZE  PARAMETERS 

NUMBER  OF  ELEMENTS  =  2096 
NUMBER  OF  NODES  =  2865 
DEGREES  OF  FREEDOM  =  8595 
SYSTEM  BANDWIDTH  =  1005 


################################################### 


## 


## 


##  LINEAR  BUCKLING  ANALYSIS## 
##  ## 
################################################### 

CONVERGED  EIGENVALUES: 

NUMBER  LAMBDA 


1  0.1782068E+03 

2  0.1763134E+03 

3  0.8268226E+02 

4  0.7179661E+02 


BUCKLING  MODE  SHAPE 
(SCALE  =  0.718E+02) 


NODE 

U 

V 

W 

1 

O.OOOOOE+00 

O.OOOOOE+00 

0 . OOOOOE+00 

2 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 

3 

O.OOOOOE+00 

O.OOOOOE+00 

0 . OOOOOE+00 

4 

O.OOOOOE+00 

O.OOOOOE+00 

0 . OOOOOE+00 

5 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 

9268 

0.67520E+00 

0.41156E-01 

-.27246E-01 

9269 

0.52133E+00 

0.32320E-01 

-.31375E-01 

9270 

0.37198E+00 

0.22749E-01 

-.33309E-01 

9271 

0.21898E+00 

0.12126E-01 

-.31172E-01 

9272 

0.83349E-01 

0.21711E-02 

-.25241E-01 
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A  graphical  depiction  of  the  deformed  shape  is  shown  in  Figure  31. 


Buckling  Mode  Buckling  Mo$le 


GLOBAL  LOCAL 

A  =  71.7966  A  =  29.8007 


Figure  31.  Global  and  local  buckling  mode  shapes. 

The  presence  of  the  delaminations  at  the  specified  ply  interfaces  is  thus  shown  to  reduce  the  critical 
buckling  load  by  66%. 

7.3  Residual  Strength  Analysis 

To  illustrate  the  execution  of  a  residual  strength  analysis,  a  model  of  an  elastically  supported  lam¬ 
inated  composite  face  plate  is  analyzed.  Due  to  the  lack  of  data  of  sufficient  resolution  regarding 
the  spatial  distribution  of  specific  internal  damage  modes  and  subsequent  experimental  determi¬ 
nation  of  residual  strength,  a  specific  state  of  internal  damage  is  assumed  for  the  following  example. 

The  geometry  and  applied  loading  is  shown  in  Figure  32.  Material  properties  were  selected  as 
S2-Glass/3501  Epoxy  tape  with  a  nominal  ply  thickness  of  0.0052  in.  The  material  elastic  proper¬ 
ties  are  given  by 


El 

=  7.150E6 

E2 

=  2.13E6 

Ez 

=  2.13E6 

G23 

=  0.71E6 

Gi3 

=  0.98E6 

Gi2 

=  0.98E6 

l'2Z 

=  0.499 

=  0.306 

=  0.296 

with  strengths  given  by 

Xren  =  2.43E5  Xcomp  =  1-77E5  Yren  =  7.0E3 

Ycomp  —  3.06E4  Zxen  —  7.03E3  Zcomp  —  3.5E4 

R  =  1.7E4  S  =  1.57E4  T  =  1.57E4 
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Figure  32.  Geometry  and  loading  of  an  elastically  supported  composite  plate. 


The  face  plate  is  composed  of  96  plies  and  is  assumed  to  have  been  subjected  to  a  normal  impact 
leading  to  an  inner  circular  region  where  all  properties  are  effectively  zero  and  the  region  is  modeled 
as  an  open  hole.  The  laminate  was  divided  into  three  layers,  where  ’|’  indicates  a  layer  interface. 
The  layup  is  given  by  [±45/0g  |  ±455/906±455/06±452/08/t452/06T455/906/t455  |  06/t45].  A 
residual  strength  analysis  is  performed  which  is  selected  to  evaluate  material  failure  with  nonlinear 
mciterial  properties.  Secant  moduli  were  computed  for  the  inplane  Young’s  moduli  as 


e[  =^i(1.0-0.75 


), 


and  for  shear  moduli  as 


G[j  =  Gijil.0-0.5 


Tii 

Pi 


where  a  and  /3  are  taken  as  material  strengths  given  by 


(20) 

(21) 


ai  =  2.43E5  02  =  7.0E3  03  =  7.03E3 

ft  =  1.7E4  ft  =  1.57E4  ft  =  1.57E4 

The  maximum  reduction  in  Young’s  moduli  was  limited  to  25%  of  the  initial  values  and  50%  for 
the  shear  moduli.  For  this  example,  the  maximum  stress  failure  criteria  is  used  with  a  nonspecific 
damage  law  that  assigns  zero  modulus  to  plies  experiencing  failure  without  regard  to  failure  mode. 
Abbreviated  input  and  output  files  are  shown,  followed  by  a  graphical  presentation  of  the  progres¬ 
sion  of  failure. 
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7.3.1  Input  Data  File  for  a  Residual  Strength  Problem 

♦HEADING 

MODEL  OF  AN  ELASTICALLY  SUPPORTED  COMPOSITE  PLATE  WITH  AN  OPEN  HOLE 
** 

**  [45/-45/0_6 I (45/-45) _5/90_6/ (45/-45) _5/0_6/ (45/-45) _2/0_8/ 

♦♦  (-45/45) _2/0_6/ (-45/45) _5/90_6/ (-45/45) _5 I 0_6/-45/45] 

♦♦  LAYUP  TAPE  THICKNESS  =  0.0052” 

** 

^SOLUTION,  METH  =  MTL 
6,15,1.05 
**ECHO 
**PREPASS 
**NODE  PRINT 
**  D  ,M 

**PLY  FAILURE  PRINT 
^NONLINEAR  MATERIAL 
10,  0.01,  2.0 
3 

2.43E5,  1.77E5 
7.0E3,  3.06E4 
7.0E3,  3.5E4 
1.7E4,  1.57E4,  1.57E4 
^parameter  lapfail 
^parameter  directory  =  current 
^MEMORY  ALLOCATION 
BAND 

^GRAPHICS,  format  =  mathematica 


2 

** 

,  2,  1.0 

NODE  1 

DEFINITIONS 

** 

*N0DE 

** 

**  elastic 

supports 

** 

1  5.0 

“■10.0 

10.0 

2  5.8 

-10.0 

10.0 

3  6.3 

-10.0 

10.0 

4  6.8 

-10.0 

10.0 

5  7.3 

-10.0 

10.0 

4708 

00 

-8.6396 

-5.9189 

4709 

8.9875 

-8.6396 

-5.9189 

4710 

9.1875 

-8.6396 

-5.9189 

4711 

9.3375 

-8.6396 

-5.9189 

4712 

9.4375 

-8.6396 

-5.9189 

4713 

9.5 

-8.6396 

-5.9189 

** 

** 

NODES,  ELEMENTS,  AND  NODE  SETS  FOR  DELAMINATED  VERTICAL 

** 

COMPONENT 

OBTAINED 

FROM  THE  FOLLOWING  *MODEL  GENERATION 
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STATEMENT: 

** 

♦MODEL  GENERATION 

** 

8,4,10.0,10.0,4,4,4,4 

** 

2. 0,1, 1,1, 0,1 

** 

0.0,  0.0,  1.0 

** 

0.0,  1.0,  0.0 

** 

1.0,  0.0,  0.0 

1.0,10.0 

** 

1,5,250,5000,5000 

0.95,9.875 

** 

2,6,255,7000,7000 

** 

0.85,9.625 

3,7,265,11000,11000 

** 

0.8,9.50 

** 

4,8,270,13000,13000 

★NODE 

5114 

0 . lOOOOOOE+02  0.3250000E+01 

0 . OOOOOOOE+00 

5115 

0 . lOOOOOOE+02  0.3198816E+01 

0.5745243E+00 

5116 

0 . lOOOOOOE+02  0.3040097E+01 

0.1149049E+01 

5117 

0 . lOOOOOOE+02  0 . 2755322E+01 

0 . 1723573E+01 

5118 

0 . lOOOOOOE+02  0.2298097E+01 

0.2298097E+01 

13517 

0.9500000E+01 

O.lOOOOOOE+02 

0.4558612E+01 

13518 

0.9500000E+01 

O.lOOOOOOE+02 

0.5918959E+01 

13519 

0.9500000E+01 

O.lOOOOOOE+02 

0.7279306E+01 

13520 

0.9500000E+01 

O.lOOOOOOE+02 

0.8639653E+01 

13521 

0.9500000E+01 

O.lOOOOOOE+02 

O.lOOOOOOE+02 

** 

**  DEFINE  GENERATOR  ELEMENTS 
♦ELEMENT,  LAYUP  =  4,  ELSET  =  300 


★  ★  top 

front 

elastic  support 

1, 

401, 

402,  902,  901,  1,  2,  102, 

101 

101, 

501, 

502,  802,  801,  401,  402,  902, 

901 

201, 

601, 

602,  702,  701,  501,  502,  802, 

801 

301, 

901, 

902,  1002,  1001,  101,  102,  202, 

201 

401, 

1001, 

1002,  1102,  1101,  201,  202,  302, 

301 

★*  top 

back  elastic  support 

501, 

1601, 

1602,  1702,  1701,  1201,  1202,  1302, 

1301 

601, 

1701, 

1702,  1802,  1801,  1301,  1302,  1402, 

1401 

701, 

1801, 

1802,  1902,  1901,  1401,  1402,  1502, 

1501 

801, 

2001, 

2002,  2102,  2101,  1801,  1802,  1902, 

1901 

901, 

2201, 

2202,  2302,  2301,  2001,  2002,  2102, 

2101 

**  bottom  front  elastic  support 
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1001, 

3601, 

3602, 

3702,  3701 

1101, 

4001, 

4002, 

4102,  4101 

1201, 

4401, 

4402, 

4502,  4501 

1301, 

3701, 

3702, 

3802,  3801 

1401, 

3801, 

3802, 

3902,  3901 

**  bottom  back  elastic  support 


1501, 

2401, 

2402, 

2502, 

2501, 

1601, 

2501, 

2502, 

2602, 

2601, 

1701, 

2601, 

2602, 

2702, 

2701, 

1801, 

3001, 

3002, 

3102, 

3101, 

1901, 

3201, 

3202, 

3302, 

3301, 

4001,  4002, 
4401,  4402, 
4601,  4602, 
4101,  4102, 
4201,  4202, 


2801,  2802, 
2901,  2902, 
3001,  3002, 
3201,  3202, 
3401,  3402, 


4102,  4101 
4502,  4501 
4702,  4701 
4202,  4201 
4302,  4301 


2902,  2901 
3002,  3001 
3102,  3101 
3302,  3301 
3502,  3501 


5|C* 


ELEMENT  GENERATION  FOR  ELASTIC  SUPPORTS 


♦ELGEN,  ELID  =  300 
1 ,  12 ,1,1,,,,,, 
101,  12,1,1,,,,,, 


1801,  12,1,1, 

y  >  9  f  9 

1901,  12,1,1, 

9  9  9  9  9 

**  CREATE  NODE  SETS  FOR  EQUIVALENCING 

** 

*NSET,  NSID  = 

10 

13,  113, 

213,  313, 

413, 

913, 

1013, 

1113 

513,  813, 

613,  713 

*NSET,  NSID  = 

20 

13300,  13319, 

13338,  13357, 

13299, 

13318, 

13337, 

13356 

13298,  13317, 

13297,  13316 

*NSET,  NSID  = 

70 

3413,  3513, 

3213,  3313, 

2813, 

2913, 

3013, 

3113 

2413,  2513, 

2613,  2713 

*NSET,  NSID  = 

80 

13487,  13506, 

13486,  13505, 

13447, 

13466, 

13485 , 

13504 

13446,  13465, 

13484,  13503 

**  EQUIVALENCE  NODE  SETS 

^EQUIVALENCE 
10,  20 
30,  40 
50,  60 
70,  80 
** 

**  ELEMENTS  IN  VERTICAL  COMPONENT  FROM  PRIOR  +MODEL  GENERATION  RUN 
♦ELEMENT,  LAYUP  =  1  ORIENTATION  =  1 
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5129 

7114 

7115 

5130 

7115 

7116 

5131 

7116 

7117 

5482 

7499 

7500 

5483 

7500 

7501 

5484 

7501 

7502 

♦ELEMENT, 

LAYUP  = 

7129 

11114 

11115 

7130 

11115 

11116 

7131 

11116 

11117 

7482 

11499 

11500 

7483 

11500 

11501 

7484 

11501 

11502 

♦ELEMENT, 

LAYUP  = 

11129 

13114 

13115 

11130 

13115 

13116 

11131 

13116 

13117 

11482 

13499 

13500 

11483 

13500 

13501 

11484 

13501 

13502 

7147 

7146 

5114 

7148 

7147 

5115 

7149 

7148 

5116 

7519 

7518 

5499 

7520 

7519 

5500 

7521 

7520 

5501 

2  ORIENTATION 

= 

11147 

11146 

7114 

11148 

11147 

7115 

11149 

11148 

7116 

11519 

11518 

7499 

11520 

11519 

7500 

11521 

11520 

7501 

3  ORIENTATION  = 

13147  13146  11114 

13148  13147  11115 

13149  13148  11116 


13519  13518  11499 

13520  13519  11500 

13521  13520  11501 


5115 

5147 

5146 

5116 

5148 

5147 

5117 

5149 

5148 

5500 

5519 

5518 

5501 

5520 

5519 

5502 

5521 

5520 

1 

7115 

7147 

7146 

7116 

7148 

7147 

7117 

7149 

7148 

7500 

7519 

7518 

7501 

7520 

7519 

7502 

7521 

7520 

1 

11115 

11147 

11146 

11116 

11148 

11147 

11117 

11149 

11148 

11500 

11519 

11518 

11501 

11520 

11519 

11502 

11521 

11520 

**  EXCLUDE  ELEMENTS  IN  ELASTIC  SUPPORTS 


^EXCLUDE  ELEMENT 
300 

**  RELAX  TOLERANCE  OF  DEFORMED  GEOMETRY  CHECKS 

♦DEFORMED  GEOMETRY 
20.0 
♦♦ 

♦♦  LAMINATE  DESCRIPTION 
♦♦ 

♦LAYER,  LAYUP  =  1 
1,  0.0052,  45.0 

1,  0.0052,  -45.0 
1,  0.0312,  0.0 

♦LAYER,  LAYUP  =  2 
1,  0.0260,  45.0 

1,  0.0260,  -45.0 
1,  0.0312,  90.0 

1,  0.0260,  45.0 

1,  0.0260,  -45.0 
1,  0.0312,  0.0 

1,  0.0104,  -45.0 
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1. 

0.0104, 

45.0 

1, 

0,0416, 

o 

o 

1. 

0.0104, 

-45.0 

1. 

0.0104, 

45.0 

1, 

0.0312, 

o 

o 

1, 

0.0260, 

-45.0 

1, 

0.0260, 

45.0 

1, 

0.0312, 

90.0 

1, 

0.0260, 

o 

LO 

1 

1, 

0.0260, 

45.0 

♦LAYER,  LAYUP=3 

1, 

0.0312, 

0.0 

1, 

0.0052, 

-  45.0 

1, 

0.0052, 

45.0 

*LAYER,  LAYUP  =  4 
2,  1.0,  0.0 
** 

**  MATERIAL  DEFINITIONS 

**  COMPOSITE  PLY  PROPERTIES  (S2-GLASS/3501  EPOXY) 
♦MATERIAL,  MATID  =  1 

7,150E6,  2.13E6,  2.13E6,  0.98E6,  0.71E6,  0.98E6 
0.306,  0.499,  0.296 
♦FAILURE  CRITERIA,  FCID  =  1 
MAX-STRESS 

2.43E5,  1.77E5,  7.0E3,  3.06E4,  7.0E3,  3.5E4 
1.7E4,  1.57E4,  1.57E4 
♦DAMAGE  LAW,  DLID  =  1 
NULL 
*♦ 

♦♦  MEOTALLIC  PROPERTIES 
♦MATERIAL,  MATID  =  2 

1.0E8,  1.0E8,  1.0E8,  0.3846E8,  0.3846E8,  0.3846E8 
0.3,  0.3,  0.3 

♦FAILURE  CRITERIA,  FCID  =  2 
MAX-STRESS 

5.3E4,  2.2E5,  5.3E4,  2.2E5,  5.3E4,  2.2E6 
4.8E4,  4.8E4,  4.8E4 
♦DAMAGE  LAW,  DLID  =  2 
NULL 

♦♦  ESTABLISH  LOCAL  COORDINATE  SYSTEM 
♦ORIENTATION 

1,  0.0,  1.0,  0.0,  0.0,  0.0,  1.0 

♦♦  NODE  SET  DEFINITION  FOR  BOUNDARY  CONSTRAINT  INPUT 
** 

♦NSET,  NSID  =  100 

1,  101,  201,  301,  401,  501,  601,  701,  801 

901,  1001,  1101 
♦NSET,  NSID  =110 
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1201, 

1301, 

1401, 

1501, 

1601,  1701, 

1801, 

1901, 

2001 

2101, 

2201, 

.2301 

♦NSET, 

NSID 

=  120 

2401, 

2501, 

2601, 

2701, 

2801,  2901, 

3001, 

3101, 

3201 

3301, 

3401, 

3501 

♦NSET, 

NSID 

=  130 

3601, 

3701, 

3801, 

3901, 

4001,  4101, 

4201, 

4301, 

4401 

4501, 

4601, 

4701 

♦BOUND ARY2 
100,  1,  3 
110,  1,  3 
120,  1,  3 
130,  1,  3 

♦*  biaxial  y-z  plane  loading 
♦CLOAD 

5300,  3,  - . lOOOOOOE+02 
5319,  3,  -.lOOOOOOE+02 
5338,  3,  -.lOOOOOOE+02 


13519,  2,  -.2000000E+02 

13520,  2,  - . 2000000E+02 

13521,  2,  -.2000000E+02 

♦ENDDATA 


7.3.2  Output  Data  File  for  a  Residual  Strength  Problem 


@0  @@ 

@®  U.S.  ARMY  RESEARCH  LABORATORY  ®® 

@@  @@ 

@®  RESTRAN  @® 

®@  @@ 

@@  RESIDUAL  STRENGTH  ANALYSIS  OF  IMPACT  @@ 
@@  DAMAGED  COMPOSITE  LAMINATES  @@ 

@@  ®® 

@@  VERSION  1.0  ®@ 

@@  ®@ 


@®@@@@@@@@@@@@@@@@@@@@@@@@@®®@@®@®@@®@@@®®®®®® 

***  MESSAGE;  ELEMENT  NODE  ORDER  IS  BEING  CONVERTED  TO  RESTRAN  FORMAT 

***  WARNING:  MATERIAL  ID  1  HAS  A  MIXED  MODE  FAILURE  CRITERIA 

ASSOCIATED  WITH  A  SINGLE  MODE  DAMAGE  LAW  (NULL) . 

NULL  ACCEPTED 

***  WARNING:  MATERIAL  ID  2  HAS  A  MIXED  MODE  FAILURE  CRITERIA 

ASSOCIATED  WITH  A  SINGLE  MODE  DAMAGE  LAW  (NULL) . 

NULL  ACCEPTED 
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###################################### 

###  ### 

###  MODEL  DEFINITION### 

###  ### 

###################################### 

JOB  HEADING: 

MODEL  OF  AN  ELASTICALLY  SUPPORTED  COMPOSITE  PLATE 

SOLUTION  CONTROL: 

SOLUTION  ALGORITHM: 

RESIDUAL  STRENGTH  ANALYSIS:  MATERIAL  FAILURE  ONLY 

ITERATION  PARAMETERS: 

MAXIMUM  NUMBER  OF  GLOBAL  ANALYSIS  CYCLES  =  6 

NUMBER  OF  ITERATIONS  TO  CONVERGE  MATERIAL  FAILURE  =  15 

MEMORY  ALLOCATION: 

BANDWIDTH  STORAGE  MODE 

*  BANDWIDTH  MINIMIZATION  WILL  BE  PERFORMED 
MATERIAL  BEHAVIOR: 

NONLINEAR  ELASTIC 

LOCAL  NONLINEAR  SOLUTION  PARAMETERS: 

MAXIMUM  NUMBER  OF  ITERATIONS  =  10 

MATERIAL  CONVERGENCE  TOLERANCE  =  O.lOOE-01 
BUCKLING  CONVERGENCE  TOLERANCE  =  0.200E+01 

PROGRAM  OPTIONS: 

*  ELEMENTS  ASSOCIATED  WITH  LOAD  APPLICATION  POINTS 
ARE  EXCLUDED  FROM  EXHIBITING  FAILURE 

*  EXTENDED  JOB  EXECUTION  STATISTICS  WILL  BE  PRINTED 

*  THE  FOLLOWING  ELEMENT  SETS  ARE  EXCLUDED  FROM  FAILURE  PREDICTION 
300 
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REQUESTED  PRINT  OPTIONS: 

=K  PRINT  ELEMENT  PLY  FAILURE  STATISTICS 

GRAPHICS: 


@  GRAPHICAL  OUTPUT  IN  MATHEMATICA  FORMAT  WILL  BE  GENERATED 

@  BOTH  UNDEFORMED  AND  DEFORMED  GEOMETRY  IS  REQUESTED 
DEFORMATION/MODE  MAGNIFICATION  FACTOR  =  O.lOOOE+01 

@  RESTRAN  GRAPHICS  CONTAINED  IN  FILE:  restran.grf 


NODE  DEFINITIONS: 


ODE  ID 

X 

Y 

Z 

1 

0.500E+01 

- . lOOE+02 

O.lOOE+02 

2 

0.580E+01 

- . lOOE+02 

O.lOOE+02 

3 

0.630E+01 

- . lOOE+02 

O.lOOE+02 

4 

0.680E+01 

-.lOOE+02 

0 . lOOE+02 

5 

0.730E+01 

-.lOOE+02 

O.lOOE+02 

13517 

0.950E+01 

O.lOOE+02 

0.456E+01 

13518 

0.950E+01 

O.lOOE+02 

0.592E+01 

13519 

0.950E+01 

O.lOOE+02 

0.728E+01 

13520 

0.950E+01 

O.lOOE+02 

0.864E+01 

13521 

0.950E+01 

O.lOOE+02 

O.lOOE+02 

S 

P  E  C  I  F  I  E 

D  NODE 

SETS: 

SET  ID 

10: 

13 

113 

213 

313 

413 

913 

1013 

1113 

513 

813 

613 

713 

SET  ID 

20: 

13300 

13319 

13338 

13357 

13299 

13318 

13337 

13356 

13298 

13317 

13297 

13316 

• 

E  L  E 

M  E  N  T 

D  E  F  I 

N  I  T  I 

0  N  S: 

ELEMENT 

N1 

N2 

N3 

N4 

N5 

N6 

N7 

N8 
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1 

401 

402 

902 

901 

1 

2 

102 

101 

2 

402 

403 

903 

902 

2 

3 

103 

102 

3 

403 

404 

904 

903 

3 

4 

104 

103 

4 

404 

405 

905 

904 

4 

5 

105 

104 

5 

405 

406 

906 

905 

5 

6 

106 

105 

11480 

13497 

13516 

13517 

13498 

11497 

11516 

11517 

11498 

11481 

13498 

13517 

13518 

13499 

11498 

11517 

11518 

11499 

11482 

13499 

13518 

13519 

13500 

11499 

11518 

11519 

11500 

11483 

13500 

13519 

13520 

13501 

11500 

11519 

11520 

11501 

11484 

13501 

13520 

13521 

13502 

11501 

11520 

11521 

11502 

E  L 

E  M  E  N 

T  S  E 

T  S: 

ELEMENT  SET  ID  300: 

1 

101  201 

301  ...  1909  1910  1911 

1912 

M  A  T  E  R  I 

AL  PROPERTY  DATA: 

MATERIAL  ID  1 

El 

E2 

E3 

=  0.715E+07 

=  0.213E+07 

=  0.213E+07 

G23  =  0.980E+06  MU23  =  0.30600 

G13  =  0.710E+06  MU13  =  0.49900 

G12  =  0.980E+06  MU12  =  0.29600 

ASSOCIATED  FAILURE  LAW: 
HOFFMAN  CRITERION  (3-D  MIXED  MODE) 
ASSOCIATED  DAMAGE  LAW: 


NULL  POST-FAILURE  PROPERTIES 

MATERIAL  ID  2 

El  =  O.lOOE+09  G23  =  0.385E+08  MU23  =  0.30000 
E2  =  O.lOOE+09  G13  =  0.385E+08  MU13  =  0.30000 
E3  =  O.lOOE+09  G12  =  0.385E+08  MU12  =  0.30000 

ASSOCIATED  FAILURE  LAW: 

MAXIMUM  STRESS  CRITERION 

XT  =  0.530E+05  XC  =  0.220E+06  YT  =  0.530E+05 
YC  =  0.220E+06  ZT  =  0.530E+05  ZC  =  0.220E+07 
R  =  0.480E+05  S  =  0.480E+05  T  =  0.480E+05 

ASSOCIATED  DAMAGE  LAW: 
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NULL  POST-FAILURE  PROPERTIES 


I  A  T  E  R  I  A 

L  LAYER  D 

E  F  I  N  I  T  I 

LAYUP  ID  = 

1 

NUMBER  OF  PLIES  =  3 

MATERIAL  ID 

PLY  THICKNESS 

ORIENTATION 

1 

0.00520 

45.00000 

1 

0.00520 

-45.00000 

1 

0.03120 

0.00000 

LAYUP  ID  = 

2 

NUMBER  OF  PLIES  =  17 

MATERIAL  ID 

PLY  THICKNESS 

ORIENTATION 

1 

0.02600 

45.00000 

1 

0.02600 

-45.00000 

1 

0.03120 

90.00000 

1 

0.02600 

45.00000 

1 

0.02600 

-45.00000 

1 

0.03120 

0.00000 

1 

0.01040 

-45.00000 

1 

0.01040 

45.00000 

1 

0.04160 

0.00000 

1 

0.01040 

-45.00000 

1 

0.01040 

45.00000 

1 

0.03120 

0.00000 

1 

0.02600 

-45.00000 

1 

0.02600 

45.00000 

1 

0.03120 

90.00000 

1 

0.02600 

0.00000 

1 

0.02600 

45.00000 

LAYUP  ID  = 

3 

NUMBER  OF  PLIES  =  3 

MATERIAL  ID 

PLY  THICKNESS 

ORIENTATION 

1 

0.03120 

0.00000 

1 

0.00520 

0.00000 

1 

0.00520 

45.00000 

LAYUP  ID  = 

4 

NUMBER  OF  PLIES  =  1 

MATERIAL  ID 

PLY  THICKNESS 

ORIENTATION 

2 

1.00000 

0.00000 
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MATERIAL  COORDINATE  SYSTEMS 


ID  Vll  V12  V13  V21  V22 

1  O.OOOE+00  O.lOOE+01  O.OOOE+00  O.OOOE+00  O.OOOE+00 


BOUNDARY  CONDITIONS: 


SINGLE  POINT  CONSTRAINTS: 

NODE  DOFl  TO  D0F2 

1  11 

101  1  1 

201  1  1 

301  1  1 


4401  1  3 

4501  1  3 

4601  1  3 

4701  1  3 


APPLIED  CONCENTRATED  FORCES: 


NODE  DIRECTION  MAGNITUDE 

5300  3  - . lOOE+02 
5319  3  -.lOOE+02 
5338  3  - . lOOE+02 


13519  2  -.200E+02 

13520  2  -.200E+02 

13521  2  -.200E+02 

FORCE/MOMENT  RESULTANTS  AT  ORIGIN: 


|R|  =  O.OOOE+00  |M|  =  O.OOOE+00 


Rx  =  O.OOOE+00  Ry  =  O.OOOE+00  Rz  =  O.OOOE+00 
Mx  =  O.OOOE+00  My  =  O.OOOE+00  Mz  =  O.OOOE+00 
Xo  =  O.OOOE+00  Yo  =  O.OOOE+00  Zo  =  O.OOOE+00 


***END  OF  MODEL  DEFINITION 
MODEL  SIZE  PARAMETERS 


V23 

O.lOOE+01 
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NUMBER  OF  ELEMENTS  =  1308 

NUMBER  OF  NODES  =  2208 

DEGREES  OF  FREEDOM  =  6624 

SYSTEM  BANDWIDTH  =  2574 

############################################### 
##  ## 

##  BEGIN  FAILURE  ANALYSIS## 
##  ## 


############################################### 


@@ 

@@  PRIMARY  ANALYSIS  CYCLE.  PASS  NUMBER 
@@ 


1 


(S® 

@@ 

@@ 


@@@@@@@@@@@@@@@@@@®@@@@@@®@@@®®®®@®®®®®®®®®®®®®®®®®®®®® 


#################################### 

###  ### 

###  NONLINEAR  SOLUTION  PROCEDURE  ### 

###  ### 

#################################### 

«  GLOBAL  ANALYSIS  CYCLE:  1  ITERATION:  1  » 

*  MATERIAL  FAILURE  LOAD  CONVERGENCE  FACTOR:  0.177E+03 

«  GLOBAL  ANALYSIS  CYCLE:  1  ITERATION:  2  » 

*  MATERIAL  FAILURE  LOAD  CONVERGENCE  FACTOR:  0.196E+00 

«  GLOBAL  ANALYSIS  CYCLE:  1  ITERATION:  3  » 

*  MATERIAL  FAILURE  LOAD  CONVERGENCE  FACTOR:  0.267E-01 

«  GLOBAL  ANALYSIS  CYCLE:  1  ITERATION:  4  » 

*  MATERIAL  FAILURE  LOAD  CONVERGENCE  FACTOR:  0.212E-01 

«  GLOBAL  ANALYSIS  CYCLE:  1  ITERATION:  5  » 

*  MATERIAL  FAILURE  LOAD  CONVERGENCE  FACTOR:  0.657E-02 

***  NONLINEAR  MATERIAL  LOADS  HAVE  CONVERGED  *** 
IN  5  ITERATIONS.  *** 


ALGORITHMIC  PATH  FOR  MINIMUM  LOAD  INCREMENT  ** 
♦♦  TO  NEXT  FAILURE:  MATERIAL  DEGRADATION  ** 
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** 

** 

MATERIAL 

BUCKLING 

** 

**  SCALE: 

0.151E+03 

EXCLUDED 

** 

** 

************5|£****5|C5t:5k***3fC+:3tf  *******  *******  ********* 


********************************* 


*  * 

*  MATERIAL  FAILURE  ANALYSIS  * 

*  * 

*  ITERATION  NO.  SCALE  FACTOR  * 

*  -  - ^ 

*  1  0.151E+03  * 

*  * 


********************************* 
*  NUMBER  OF  ELEMENTS  DEGRADED  =  3 


********************************* 


*  ♦ 

*  MATERIAL  FAILURE  ANALYSIS  ★ 

♦  ♦ 

*  ITERATION  NO.  SCALE  FACTOR  * 

*  :jc 

*  2  0.151E+03  * 

*  * 


********************************* 
*  NUMBER  OF  ELEMENTS  DEGRADED  =  14 

******************************  **34e 


*  ♦ 

*  MATERIAL  FAILURE  ANALYSIS  ♦ 

*  * 

*  ITERATION  NO.  SCALE  FACTOR  * 

*  -  - * 

*  3  0.151E+03  ♦ 

*  * 


********************************* 
*  NUMBER  OF  ELEMENTS  DEGRADED  =  7 

********************************* 


*  * 

*  MATERIAL  FAILURE  ANALYSIS  * 

*  * 

*  ITERATION  NO.  SCALE  FACTOR  * 

*  -  - * 

*  4  0.151E+03  ★ 

*  * 


********************************* 


*  NUMBER  OF  ELEMENTS  DEGRADED  =  6 
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jjt:^:*********  **************  ******** 


MATERIAL  FAILURE  ANALYSIS 
ITERATION  NO.  SCALE  FACTOR 
5  0.151E+03 


**s|t:(e3»t:f::Jeatc:4c*:f::<«****  *******  ********** 


*  NUMBER  OF  ELEMENTS  DEGRADED  =  10 

3jCJtC3tC5|c:<C*  *****************  ********** 


*  * 

*  MATERIAL  FAILURE  ANALYSIS  ♦ 

*  * 

*  ITERATION  NO.  SCALE  FACTOR  * 

:,c  -  - * 

*  6  0.151E+03  * 

*  * 


********************************* 


*  NUMBER  OF  ELEMENTS  DEGRADED  = 


20 


********************************* 


MATERIAL  FAILURE  ANALYSIS 


* 
* 
* 

ITERATION  NO.  SCALE  FACTOR  ♦ 

-  - * 

7  0.151E+03  ♦ 

* 


********************************* 


*  NUMBER  OF  ELEMENTS  DEGRADED  =  54 

********************************* 


*  * 

*  MATERIAL  FAILURE  ANALYSIS  * 

*  * 

*  ITERATION  NO .  SCALE  FACTOR  * 

^  -  - * 

*  8  0.151E+03  * 

*  * 


********************************* 


*  NUMBER  OF  ELEMENTS  DEGRADED  =  140 

********************************* 
*  * 

*  MATERIAL  FAILURE  ANALYSIS  * 

*  * 

*  ITERATION  NO.  SCALE  FACTOR  * 

*  -  - * 
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*  9  0.151E+03  * 

*  * 

*  NUMBER  OF  ELEMENTS  DEGRADED  =  188 

*  ELEMENT  FAILURE  HAS  ALTERED  MODEL  STABILITY  SUCH 
THAT  RIGID  BODY  MODES  HAVE  BEEN  DETECTED. 

TOTAL  FAILURE  IS  ASSUMED. 

ELEMENT  PLY  FAILURE  STATUS  AT  CYCLE  NO  1 

ELEMENT  PLY  FAILURE  MODES: 

ID  ’/.FIBER  ’/.MATRIX  ’/.BUCKLING  ’/.TOTAL 

5129  0  0  0  100 

5130  00  0  100 

5131  0  0  0  33 

5132  0  0  0  33 


11461 

11462 

11463 

11464 

############################################### 


##  ## 

##  ANALYSIS  PREDICTS  CATASTROPHIC  FAILURE  OF  ## 

##  THE  MODEL  AT  AN  ULTIMATE  LOAD  GIVEN  BY:  ## 

##  ## 

##  P(ULT)  =  (0.15140E+03)  *  P(INITIAL)  ## 

##  ## 


############################################### 


000  100 

000  100 

000  100 

0  0  0  33 


7.3.3  Graphical  Output  Using  MATHEMATICA 

The  initial  graphics  generated  by  MATHEMATICA  show  the  color  codes  used  to  display  various 
element  failure  modes.  In  the  specific  demonstration  problem  analyzed,  complete  element  failure 
is  indicated  by  a  wireframe  depiction;  because  no  specific  failure  modes  were  predicted  using  the 
maximum  stress  criterion,  material  failure  is  shown  using  the  general  material  failure  color  coding 
scheme. 
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Equilibrium  between  external  and  internal  forces  was  converged  in  six  iterations.  With  the  material- 
only  failure  option  set,  in  each  element  all  plies  were  analyzed  using  the  maximum  stress  criterion 
to  determine  the  scale  factor  to  by  applied  to  the  external  loads  to  cause  failure  in  a  specific  ply. 
The  lowest  factor  was  selected,  multiplied  by  the  acceleration  factor  (which  in  this  sample  run  was 
selected  as  1.05  or  5%  above  the  minimum  ply  failure  load),  and  then  applied  to  the  external  loads. 
During  the  first  analysis  cycle,  the  initial  failure  of  three  plies  was  sufficient  to  cause  a  sequence 
of  load  redistribution  and  subsequent  additional  failures  that  led  to  a  cascade  process  resulting  in 
ultimate  total  failure  as  listed  in  Table  3. 
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Table  3.  Cascade  failure  in  residual  strength  analysis. 


Iteration 

Element  Failures 

1 

3 

2 

14 

3 

7 

4 

6 

5 

10 

6 

20 

7 

54 

8 

140 

9 

188 

This  was  automatically  determined  by  the  presence  of  fragmentation  which  caused  the  global  stiff¬ 
ness  matrix  to  become  singular,  which  terminated  execution.  The  final  failure  state  is  shown  below. 
Mediating  element  failures  are  shown  to  have  progressed  to  the  outer  boundary.  Failure  in  the  line 
of  elements  comprising  the  outer  boundary  have  been  precluded  using  the  LAPFAIL  parameter, 
which  prevents  elements  to  which  applied  external  loads  have  been  applied  from  exhibiting  failure. 


ANALYSIS  CYCLE  1 

Material  Failure  at  Scale  =  0 . 151E+03 


Y 

Failure  State 


.V 


The  final  residual  strength  prediction  yields  the  ultimate  load  carrying  capability  as  a  multiple 
of  the  initial  biaxial  loads  equal  to  151.4.  Because  a  5%  acceleration  was  used,  the  resolution  of 
the  residual  strength  measure  R  falls  within  the  bounds  144.19  <  R<  151.4. 


8  Conclusion 

The  RESTRAN  computer  program  is  a  finite  element  based  design  tool  for  the  analysis  of  resid¬ 
ual  strength  in  composite  structures  with  arbitrary  three-dimensional  geometry,  loading,  support 
conditions,  and  material  properties.  The  unique  feature  of  RESTRAN  is  the  combined  algorithmic 
accounting  of  both  material  and  structural  failure  modes.  Material  failure  modes  axe  predicted  us¬ 
ing  a  robust  suite  of  failure  criteria  and  damage  laws.  Structural  failure  due  to  sublaminate  buckling 
of  multiple  delaminations  is  performed,  allowing  for  an  arbitrary  number  of  assumed  delaminations. 
User-defined  subroutine  interfaces  are  provided  to  permit  user  modification  and  enhancement  of 
basic  analysis  procedures.  A  progressive  failure  analysis  is  performed  until  ultimate  structural 
failure  is  simulated,  yielding  an  estimate  of  the  residual  strength.  Thus,  the  RESTRAN  program 
constitutes  a  robust  design/research  tool  which  can  provide  an  accurate  assessment  of  remaining 
strength  in  composite  structures  containing  impact  damage. 
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Appendix  A: 

Sample  User-Defined  Subroutine  to  Compute  Nonlinear  Material  Properties 
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oo  oooo  oooooooo  ooooooooooo  o  oooo 


c 

c 

c 

c 


SUBROUTINE  USERNL  (  SIRES S,STRA IN, SXX,SYY,SZZ,TYZ,TZX,TXY,EXX, 

1  EYY , EZZ , G YZ , GZX , GXY , THKN , ATHK , TLM , THETA , 

2  PNLRM , FE 1 , FE2 , FES , FG3 1 , FG23 , FG 12 , FV3 1 , 

3  FV23,FV12,N0RD  ) 

**  ** 

**  USER-DEFINED  ROUTINE  TO  COMPUTE  SECANT  ** 

**  MODULI  IN  NONLINEAR  MATERIAL  ANALYSIS  ** 

** 

ie:«e:)c*3t(*if:****************:^*N<:*:4c***3^***************3(c* 


IMPLICIT  DOUBLE  PRECISION (A-H,0-Z) 

DIMENSION  STRESS (27, 6) , STRAIN (27, 6) , PNLRM (25) 

THREE  DIFFERENT  USER-DEFINED  NONLINEARITIES  ARE  CONTAINED  @@@ 
@@@  IN  THIS  SUBROUTINE.  THE  FIRST  SPECIFIES  A  PIECE-WISE  LINEAR  @@@ 
@@@  LINEAR  STRESS-STRAIN  RELATIONSHIP,  THE  SECOND  DEFINES  @@@ 
m@  UNEQUAL  MODULI  IN  TENSION  AND  COMPRESSION.  THE  THIRD  CASE  @@@ 
YIELDS  MODULI  THAT  ARE  A  CONTINUOUS  FUNCTION  OF  STRESSES  @@@ 
@@@  UP  TO  SET  LEVEL.  THESE  CASES  ARE  SELECTED  THROUGH  SETTING  @@@ 
THE  FIRST  ELEMENT  OF  THE  PNLRM  ARRAY  WHICH  CONTAINS,  IN 
@©@  SEQUENCE,  THE  PARAMETERS,  Pn,  SPECIFIED  IN  THE  ^NONLINEAR  @©© 
MATERIAL  INPUT  STATEMENT. 

IF  (  INT(PNLRM(1))  .EQ.  1  )  THEN 

@@@  USER-DEFINED  PIECE-WISE  LINEAR  STRESS-STRAIN  LAW  ©©© 

@@@ 

PNLRM (2, 3, 4)  :  CONTAINS  YIELD  STRESS  VALUES  @©@ 

©@©  PNLRM (5, 6, 7)  :  CONTAINS  MODULI  REDUCTION  FACTOR  @@@ 

©@©  COMPUTE  VON  MISES  STRESS  ©©© 

VMSS  ^  SQRT((  (SXX-SYY)**2  +  (SYY-SZZ)**2  +  (SZZ-SXX)**2)/2 .0) 

@©©  DETERMINE  WHICH  LINEAR  SECTION  OF  THE  @@@ 

STRESS/STRAIN  CURVE  WE  ARE  IN  ©@© 

IF  (  VMSS  .LT.  PNLRM (1)  )  THEN 


C 

C 

C 

C 

C 

C 

C 


©©@  INITIAL  MODULI 
FAC  =  1 

ELSE  IF  (  VMSS  .GE.  PNLRM(l)  .AND.  VMSS  .LT.  PNLRM(2)  )  THEN 
©Q©  SECOND  LINEAR  SEGMENT  @©© 

FAC  =  PNLRM (4) 

ELSE  IF  (  VMSS  .GE.  PNLRM (2)  .AND.  VMSS  .LT.  PNLRM (3)  )  THEN 
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C  @(§@  THIRD  LINEAR  SEGMENT  @@@ 

C 

FAC  =  PNLRM(5) 

C 

ELSE  IF  (  VMSS  .GE.  PNLRM(3)  )  THEN 
C 

C  @®@  FINAL  LINEAR  SEGMENT  @@@ 

C 

FAC  =  PNLRM(5) 

C 

END  IF 
C 

C  @@@  ASSIGN  MODULI  REDUCTION  FACTORS 
C 

FEl  =  FAC 
FE2  =  FAC 
FES  =  FAC 
FG31  =  FAC 
FG23  =  FAC 
FG12  =  FAC 
C 

ELSE  IF  (  INT(PNRLM(1))  .EQ.  2  )  THEN 
C 


C  @@@  UNEQUAL  TENSION  AND  COMPRESSION  MODULI  @@@ 
C  @00  @0® 
C  @00  PNRLM(2-7)  CONTAIN  MULTIPLIERS  TO  YIELD  TENSION  MODULI  @00 
C  00©  VALUES  FOR  Exx,Eyy,Ezz,Gzx,Gyz,Gxy  @00 
C  @00  PNRLM(8-15)  CONTAIN  MULTIPLIERS  TO  YIELD  COMPRESSION  000 


C  000  MODULI  VALUES  FOR  Exx.Eyy ,Ezz,Gzx,Gyz,Gxy  000 

C 

IF  (  SXX  .GT.  O.ODO  )  THEN 
FEl  =  PNRLM(2) 

ELSE 

FEl  =  PNRLM(8) 

END  IF 

IF  (  SYY  .GT.  O.ODO  )  THEN 
FE2  =  PNRLM(3) 

ELSE 

FE2  =  PNRLMO) 

END  IF 

IF  (  SZZ  .GT.  O.ODO  )  THEN 
FES  =  PNRLM(4) 

ELSE 

FES  =  PNRLM(IO) 

END  IF 

IF  (  TYZ  .GT.  O.ODO  )  THEN 
FG31  =  PNRLM(5) 

ELSE 

FG31  =  PNRLM(ll) 

END  IF 

IF  (  TZX  .GT.  O.ODO  )  THEN 
FG23  =  PNRLM(6) 

ELSE 

FG32  =  PNRLM(12) 
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o  o 


c 

c 

c 

c 

c 

c 

c 

c 

c 


END  IF 

IF  (  TXY  .GT.  O.ODO  )  THEN 
FG12  =  PNRLM(7) 

ELSE 

FG12  =  PNRLM(13) 

END  IF 


ELSE  IF  (  INT(PNLRM(1))  .EQ.  3  )  THEN 


@@@  USER-DEFINED  NONLINEAR  STRESS-STRAIN  LAW  @®@ 


@@@  PNLRM(2-11) 


CONTAINS  ULTIMATE  STRENGTHS  @0® 
SXX-TEN,  SXX-COMP,  ®®@ 

SYY-TEN,  SYY-COMP,  ®@® 

SZZ-TEN,  SZZ-COMP,  @®@ 

R,  S,  T 


@@®  REDUCTION  LAW:  FACi  =  1.0  -  BETA*(Si/STRENGTHi)  @@® 

C 

BETAl  =0.75 
BETA2  =0.50 
IF  (  SXX  .GE.  0.0  )  THEN 

FEl  =  1.0  -  BETA1*ABS(SXX/PNLRM(2)) 

ELSE 

FEl  =  1.0  -  BETA1*ABS(SXX/PNLRM(3)) 

END  IF 
C 

IF  (  SYY  .GE.  0.0  )  THEN 

FE2  =  1.0  -  BETA1*ABS(SYY/PNLRM(4)) 

ELSE 

FE2  =  1.0  -  BETA1*ABS(SYY/PNLRM(5)) 

END  IF 
C 

IF  (  SZZ  .GE.  0.0  )  THEN 

FES  =  1.0  -  BETA1*ABS(SZZ/PNLRM(6)) 

ELSE 

FES  =  1.0  -  BETA1*ABS(SZZ/PNLRM(7)) 

END  IF 
C 

FG31  =  1.0  -  BETA2*ABS(TYZ/PNLRM(8)) 

FG23  =  1.0  -  BETA2*ABS(TZX/PNLRM(9)) 

FG12  =  1.0  -  BETA2*ABS(TXY/PNLRM(10)) 

C 

C  ®®@  LIMIT  MODULI  REDUCTION  TO  FACTOR  OBTAINED  TO  25*/.  @®® 
C  e®@  FOR  YOUNG’S  MODULI  AND  50*/.  FOR  SHEAR  MODULI  ®®® 

C 

IF  (  FEl  .LT.  0.75  )  FEl  =  0.75 

IF  (  FE2  .LT.  0.75  )  FE2  =  0.75 

IF  (  FES  .LT.  0.75  )  FES  =  0.75 

IF  (  FG31  .LT.  0.50  )  FG31  =  0.50 

IF  (  FG23  .LT.  0.50  )  FG23  =  0.50 

IF  (  FG12  .LT.  0.50  )  FG12  =  0.50 

C 

END  IF 
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END  EXAMPLE 


RETURN 

END 


Appendix  B: 

Sample  User-Defined  Subroutine  to  Compute  Material  Failure 
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o  o 


C 

C 

C 

C 

C 

C 

C 

C 

C 


SUBROUTINE  USERFC  (  SOFTS, EGPTS , STRESS, STRAIN, PTHK, THETA, IPLY, 

1  NELID , PARAM , El , E2 , E3 , 012 , G23 , 031 , V12 , V23 , 

2  V 1 3 , EM , EF , NORD , NSTAT , NFSTAT , NACCS  ) 

************************  *********** 

:k*  ** 

**  USER-DEFINED  FAILURE  CRITERION  TO  COMPUTE  ** 

**  MATERIAL  FAILURE  UNDER  APPLIED  LOAD  ** 

**  ** 

IMPLICIT  DOUBLE  PRECISION (A-H,0-Z) 

DIMENSION  SGPTS(27,6) ,EGTPS(27,6) ,STRESS(6) ,STRAIN(6) 
DIMENSION  PARAM (15) 


###  BEGIN  EXAMPLE  ### 

C 

C  «<  VALUE  IN  PARAM  (1)  IS  USED  TO  SELECT  THE  »> 

C  «<  FAILURE  CRITERION  TO  BE  APPLIED.  OTHER  »> 

C  «<  VALUES  STORED  IN  PARAM  ARRAY  AS  NEEDED.  »> 

C 

NFLAG  =  INT(PARAM(1)) 

C 

IF  (  NFLAG  .EQ.  1  )  THEN 
C 

C  «<  FAILURE  CRITERIA  AND  DAMAGE  LAWS  AS  PRESENTED  IN:  »> 

C  «<  FU-KUO  CHANG  AND  KUO-YEN  CHANG,  "A  PROGRESSIVE  DAMAGE  »> 

C  «<  MODEL  FOR  LAMINATED  COMPOSITES  CONTAINING  STRESS  »> 

C  «<  CONCENTRATION,"  J.  COMP.  MAT,  PP.  834-855,  1987.  »> 

C 

51  =  STRESS (1) 

52  =  STRESS (2) 

S12  =  STRESS (6) 

C 

C  «<  NONLINEAR  SHEAR  STRESS-STRAIN  FACTOR  ALPHA  »> 

C 

ALPHA  =  PARAM (2) 

C 

C  «<  ULTIMATE  STRENGTH  MEASURES  »> 

C 

XT  =  PARAMO) 

YT  =  PARAM (4) 

SC  =  PARAM (5) 

C 

EM  =  SQRT(  (S2/YT)**2  +  (S12**2/(2*G12)+3*ALPHA*S12**4/4)/ 
1  (SC*=t‘2  /(2*G12)+3’i‘ALPHA*  SC**4/4)  ) 

C 

C  «<  FIBER-MATRIX  SHEARING  AND  FIBER  BREAKAGE  »> 

C 

EF  =  SQRT(  (S1/XT)**2  +  (S12**2/(2*G12)+3*ALPHA*S12**4/4)/ 
1  (SC**2  /(2*G12)+3*ALPHA*  SC**4/4)  ) 

C 

ET  =  MAX(EM,EF) 
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EF  =  1.0/ET 
EM  =  EF 
C 

NFSTAT  =  NSTAT 
NACCS  =  0 
C 

C  «<  EMPLOY  PROPERTY  DEGRADATION  LAVS  »> 

C 

IF  (  EM  .GT.  1.0  )  THEN 
C 

C  «<  MATRIX  CRACKING  »> 

C 

V12  =  0.0 
V23  =  0.0 
V13  =  0.0 
E2  =  0.0 
E3  =  0.0 
G23  =  0.0 
C 

IF  (  NSTAT  .EQ.  0  )  THEN 
NFSTAT  =  2 
ELSE 

NFSTAT  =  3 
END  IF 
NACCS  =  1 
C 

END  IF 
C 

IF  (  EF  .GT.  1.0  )  THEN 
C 

C  <«  FIBER  BREAKAGE/FIBER-MATRIX  SHEARING  DAMAGE  CALCULATED  »> 

C  «<  USING  A  MICROMECHANICAL  FIBER  BUNDLE  FAILURE  APPROACH  »> 

C 

A  =  PTHK 
AO  =  PARAM(6) 

BETA  =  PARAM(7) 

«<  DAMAGE  IS  ASSUMED  TO  FOLLOW  A  WEIBUL  DISTRIBUTION  »> 

C 

El  =  El*  EXP (- (A/AO) **BETA) 

G12  =  G12*EXP(-(A/A0)**BETA) 

G31  =  G31*EXP(-(A/A0)**BETA) 

C 

IF  (  NSTAT  .EQ.  0  )  THEN 
NFSTAT  =  2 
ELSE 

NFSTAT  =  3 
END  IF 
NACCS  =  1 
C 

END  IF 
C 

END  IF 
C 
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c  ***  END  EXAMPLE  *** 
C 

RETURN 

END 
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Appendix  C: 

Sample  User-Defined  Subroutine  to  Assign  Initial  Material  Damage 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  USERID  (PTHK, THETA, El, E2,E3,G13,G23,Gi2,V13,V23, 

1  V12,P1,P2,P3,P4,P5,P6,P7,P8,P9,NELID, 

2  IPLY,NSTAT  ) 

>|c*  ** 

**  USER-DEFINED  SUBROUTINE  TO  ASSIGN  ** 

**  INITIAL  ELEMENT  MATERIAL  DAMAGE.  ** 

stesjc 

:JC*3<Cj|c3jCJiciC5»C****3|C***J|eSit**3|t***5t(**5|e**2t£****************** 

IMPLICIT  DOUBLE  PRECISION (A-H,0-Z) 

***  EXAMPLE  ROUTINE  *** 

USER  DEFINITION  FOR  THE  INITIAL  PROPERTIES  OF 
Kth  PLY  IN  ELEMENT  WITH  ID  LABEL  NELID: 

SELECTION  OF  INITIAL  DAMAGE  IS  BASED  ON  ELEMENT 
ID  SERIES  AND  SPECIALIZED  MATERIAL  PROPERTY 
DEGRADATION  RELATIONS  USING  INPUT  PARAMETERS 
PI  ->  P9 

IF  (  NELID  .GE.  2000  .AND.  NELID  .LT.  3000  )  THEN 

PARTIAL  FAILURE  I .  TRANSVERSE  MODULI  REDUCTION 

E2  =  P1*E2 
G12  =  P1+G12 
G23  =  P2*G23 

NSTAT  =  2 

ELSE  IF  (  NELID  .GE.  3000  .AND.  NELID  .LT.  4000  )  THEN 

PARTIAL  FAILURE  II.  TRANSVERSE  INPLANE  AND  NORMAL  ®@@ 

@@@  MODULI  REDUCTION  @@@ 

E2  =  P3*E2 
E3  =  P5*E3 
G12  =  P3*G12 
G23  =  P4*G23 
G13  =  P6*G13 

NSTAT  =  2 

ELSE  IF  (  NELID  .GE.  4000  .AND.  NELID  .LT.  5000  )  THEN 


PARTIAL  FAILURE  III.  TRANSVERSE  INPLANE  AND 
@0®  NORMAL  MODULI  REDUCTION 

E2  =  P3*E2 
E3  =  P5*E3 
G12  =  P3*G12 
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G23  =  P4*G23 
G13  =  P6*G13 
C 

NSTAT  =  2 
C 

ELSE  IF  (  NELID  .GE.  5000  .AND.  NELID  .LT.  7000  )  THEN 
C 

C  TOTAL  FAILURE.  REDUCE  ALL  MATERIAL  PROPERTIES  @@@ 

C  TO  ZERO 

C 

El  =  0.0 
E2  =  0.0 
E3  =  0.0 
G13  =  0.0 
G23  =0.0 
G12  =  0.0 
V13  =  0.0 
V23  =  0.0 
V12  =  0.0 
C 

NSTAT  =  3 
C 

END  IF 
C 

END  EXAMPLE 


RETURN 

END 
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Appendix  D: 

Sample  User-Defined  Subroutine  to  Interpret  Buckling  Modes 
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no 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 


SUBROUTINE  USERBl  (  EMODE,NELEM,NDSET,NXSET,NSETI ,NDID,EIGV, 

1  NUMNOD , NBPRC , LCNSL , ITERB , NUMEXC ,  NELNUM , 

2  NACC , NSURF , NINC , MDXO , MDXl , MAXB , MDX3  ) 

*******  sic**  :ie*:ic**:<e**:»c*****j|e***3}c***:<t*****  *  +  **:»:**** 

**  ** 

**  USER-DEFINED  SUBROUTINE  TO  EVALUATE 

**  BUCKLING  FAILURE  MODE  ** 

**  ** 

************************************************* 

IMPLICIT  DOUBLE  PRECISI0N(A-H,O-Z) 

DIMENSION  EMODE (MDXO, MAXB) ,NELEM(MDX3, 12) ,NDSET(MDX1) 
DIMENSION  NXSET(MDXl) 

###  BEGIN  EXAMPLE  ### 


@@@  DESCRIPTION:  FOR  DELAMINATION  NODE  SETS  10  AND  20,  ELEMENT  @@@ 
@@@  PAIRS  HAVE  BEEN  IDENTIFIED  ON  EITHER  DELAMINATION  FACE  WHICH 
Q@@  WILL  INDICTATE  BUCKLING  IF  THEY  SHOW  THE  SAME  SIGN  OF  @@@ 

000  RELATIVE  MOTION  AND  IF  THE  RELATIVE  DEFORMATION  IS  GREATER  000 

000  THAN  0.4.  THIS  IS  BASED  ON  THE  FACT  THAT  THE  MODE  SHAPE  IS  @00 

@00  NORMALIZED  AND  THAT  A  PARTICIPATION  OF  THE  SELECTED  ELEMENTS  000 

00©  GREATER  THAN  407.  OF  THE  MAXIMUM  DEFLECTION  IS  TAKEN  AS  A  @00 

000  SIGNIFICANT  PARTICIPATION  IN  THE  MODE.  THE  SIGNS  OF  THE  000 

@00  DISPLACEMENTS  INDICATING  AN  OPENING  OR  CLOSING  MODE  ARE  @00 

@00  ARBITRARY  SINCE  THE  MAGNITUDE  OF  THE  MODE  SHAPE  IS  @00 

@00  INDETERMINATE.  @00 


NACCS  =  0 

IF  (  NDID  .EQ.  10  )  THEN 

000  DELAMINATION  ID  10  IS  BEING  PROCESSED.  BUCKLING  WILL  BE  @00 
000  BASED  ON  THE  MAGNITUDE  OF  RELATIVE  DISPLACEMENT  BETWEEN  @00 
@00  OPPOSING  ELEMENTS  2050  AND  3050  WHICH  ARE  LOCATED  ABOVE  000 
000  AND  BELOW  THE  DELONANE.  @00 

@00  OBTAIN  THE  FOLLOWING  DISPLACEMENTS:  @00 
000  NODE  1  IN  ELEMENT  2050  @00 

@00  NODE  5  IN  ELEMENT  3050  @00 

@00  LOCATE  ELEMENTS  IN  ARRALEEQ.  2050  )  NL2050  =  I 
IF  (  NELEM(I,1)  .EQ.  3050  )  NL3050  =  I 
END  DO 

@00  CONVERT  INTERNAL  NODE  NUMBERS  IN  NELEM  ARRAY  TO  @00 
@00  INPUT  NODE  NUMBERS  00© 

NODEl  =  NSETI(NELEM(NL2050,2)) 

NODES  =  NSETI(NELEM(NL3050,6)) 

000  LOCATE  NODE  MODAL  DEFLECTIONS  IN  EMODE  ARRAY.  TEST  @00 
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C  @@@  ON  INPUT  NODE  NUMBER  AND  GENERATED  COINCIDENT  NODE  @@@ 

C  @@@  NUMBER  STORED  IN  SAME  RECORD  @@@ 

C 

NPNTl  =  0 
NPNT5  =  0 
DO  I  =  1,  NUMNOD 

IF  (  NODEl  .EQ.  INT(EMODE(I,l))  -OR. 

1  NODEl  .EQ.  INT(EM0DE(I,5))  )  NPNTl  =  I 

IF  (  NODES  .EQ.  INT(EMODE(I , 1) )  .OR. 

1  NODES  .EQ.  INT(EMODE(I,S))  )  NPNTS  =  I 

END  DO 
C 

C  OBTAIN  W-COMPONENT  FOR  COMPARISON 

C 

IF  (  NODEl  .EQ.  INT(EMODE(NPNTl,l))  )  W1  =  EMODE(NPNTl ,4) 

IF  (  NODEl  .EQ.  INT(EM0DE(NPNT1,5))  )  W1  =  EM0DE(NPNT1,8) 

IF  (  NODES  ,EQ.  INT (EMODE (NPNTS, 1) )  )  WS  =  EM0DE(NPNT5,4) 

IF  (  NODES  .EQ.  INT (EMODE (NPNTS, S) )  )  WS  =  EMODE (NPNTS, 8) 

C 

IF  (  ABS(Wl-WS)  .GT.  0.40  )  THEN 
C 

@0©  LOCAL  BUCKLING  IS  DETERMINED  @@@ 

NACC  =  1 
NSURF  =  1 
C 

ELSE 

C 

NACC  =  0 
C 

END  IF 
C 

END  IF 
C 

RETURN 

END 
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Appendix  E: 

Sample  User-Defined  Subroutine  to  Assign  Material  Degradation  Due  to  Buckling 
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SUBROUTINE  USERB2  (  NELEM,NODSET,NDELM,NDID,NELFL,ELFAC,NFAIL, 
1  NUMELF,NSURF,MDX1,MDX2,MDX3,MDX4  ) 

C 

Q  jfer+e***  ********  *******************  +  *************** 

c  **  ** 

C  **  USER-DEFINED  SUBROUTINE  TO  DETERMINE 

C  **  ELEMENTS  FAILED  DUE  TO  BUCKLING  AND  ** 

C  **  ASSIGNED  REDUCED  MATERIAL  PROPERTIES  ** 

C  **  ** 

C  *♦*♦**♦♦****♦*  :|t***!|!i(C***!(C****  ****♦****!(! 

c 

IMPLICIT  DOUBLE  PRECISION (A-H,0-Z) 

C 

DIMENSION  NELEM(MDX3, 12) ,N0DSET(MDX2) 

DIMENSION  NELFL(MDXl)  ,ELFAC(MDX1)  ,NDELM(MDX4) 

C 

NUMELF  =  0 
NELFL(l)  =  0 
ELFAC(l)  =  0.0 
C 

C  ###  BEGIN  EXAMPLE  ### 

C 

C  @@@  ASSIGN  DIFFERENT  LEVELS  OF  MATERIAL  FAILURE  DUE  @@@ 

C  @@@  BUCKLING  DEPENDING  ON  THE  DELAMINATION  LOCATION  (9©@ 

C  @©@  WHICH  IS  SIGNIFIED  BY  THE  DELAMINATION  ID  ©0© 

C 

IF  (  NDID  .EQ.  10  )  THEN 
C 

C  ©©©  PARTIAL  FAILURE  ©@© 

C 

C  ©©@  REDUCE  MODULI  BY  50%  IN  ELEMENTS  5000  THRU  5500.  ©©© 

C  ©©@  NUMBER  OF  FAILED  ELEMENTS  RETURNED  IN  NUMELF  ©©© 

C 

NUMELF  =501 
DO  I  =  1,  501 

NELFL(I)  =  4999+1 
ELFAC(l)  =0.50 
END  DO 
C 

END  IF 
C 

RETURN 

END 
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Appendix  F: 

Sample  Generated  MATHEMATICA  Input  File 
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(* 


{RGBColorC  1.0,  1.0,  1.0] ,Rectangle[{  0.0, 0.0}, {10.0,  6.0}]}, 
{RGBColorC  0.0,  0.5,  1.0]  , Rectangle [{  0.0, 3.0}, {10.0,  4.0}]}, 
Text["RESTRAN  GRAPHICAL  RESULTS" ,{5.0, 3. 5}] 

} 

Show [Graphics [  fcode  ]  ] 
model 1  =  { 

{RGBColorE  1.0,  1.0,  1.0] ,Rectangle[{  0.0, 0.5}, {10.0,  5.0}]}, 
{Text [FontFormE" ANALYSIS  CYCLE  1  ", 

{ " Courier-Bold" , 8}] , {5 . 0 , 4 . 0} , {0 , - 1}] } , 

{Text [FontForm["Linear  Buckling  Analysis", 

{ " Courier-Bold" , 8}] , {5 . 0 , 3 . 5}] } , 

Line [{{0.0, 2. 75}, {10. 0,2. 75}, {10. 0,4. 5}, {0.0, 4. 5}, {0.0, 2. 75}}] 

} 

Show [Graphics [modell , 

Frame  ->  False  ]  ] 

(* 

@00  UWDEFORMED  MODEL  GEOMETRY  000 
@00  SHOWING  MATERIAL  FAILURE  @00 

*) 

ppntsl  =  {{  0.0000,  0.0000,  0.9000},{  2.0000,  0.0000,  0.9000}, 

{  2.0000,  0.0000,  1.0000},{  0.0000,  0.0000,  1.0000}} 

ppnts2  =  {{  0.0000,  0.0000,  0.9000},{  0.0000,  0.0000,  1.0000}, 

{  0.0000,  1.0000,  1.0000},{  0.0000,  1.0000,  0.9000}} 

ppnts3  =  {{  0.0000,  0.0000,  0.9000},{  0.0000,  1.0000,  0.9000}, 

{  2.0000,  1.0000,  0.9000},{  2.0000,  0.0000,  0.9000}} 

ppnts4  =  {{  0.0000,  0.0000,  1.0000},{  2.0000,  0.0000,  1.0000}, 

{  2.0000,  1.0000,  1.0000},{  0.0000,  1.0000,  1.0000}} 


model2  =  { 

{FaceForm [RGBColorE. 00, .00, .99] , RGBColorE. 00, .00, .99]] , Polygon [ppntsl]}, 
{FaceForm [RGBColorE. 00, .00, .99] , RGBColorE. 00, .00, .99]] , Polygon [ppnts23}, 
{FaceForm [RGBColorE. 00, .00, .99] , RGBColorE. 00, .00, .99]] ,PolygonEppnts3]}, 


} 

ShowE  Graphics3DEmodel2, 
Boxed  ->  False, 
Axes  ->  True, 

Ticks  “>  None, 
Lighting  ~>  False, 
Render All  ->  True, 
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Viewpoint  ->  {1.3, -2. 4, 2.0}, 

ViewVertical  ->  {0,0,1}, 

DefaultFont  ->  {"Times-Italic" ,6} , 
PlotLabel  ->  FontForm ["Initial  Geometry", 
{"Helvetica-Bold" ,9}]  , 
AxesLabel  ->  {"X" , "Y" , "Z"}  ]  ] 


(* 


Clear All [ppntsl ,ppnts2 ,ppnts3 ,ppnts4 ,ppnts5] 
DEFORMED  MODEL  GEOMETRY 


*) 


THE  FOLLOWING  AMPLITUDE  FACTOR,  GAMMAi,  CAN  BE 
CHANGED  TO  GIVE  THE  DESIRED  DEFORMED  SHAPE 

gamma 1  =  -1.5000 


tpntsl  =  {0.0000,0.0000,0.9000}  +  gaimnal{0.0000,0. 0000,-1 .0000} 
tpnts2  =  {2.0000,0.0000,0.9000}  +  gainmal{0. 0188,0. 0000,-0. 6335} 
tpnts3  =  {2.0000,0.0000,1.0000}  +  gainmal{-.0182,0. 0000,-0. 6335} 
tpnts4  =  {0.0000,0.0000,1.0000}  +  gaimnal{0. 0000,0.0000,-1 .0000} 
ppnts223  =  {tpntsl  ,tpnts2  ,tpnts3  ,tpnts4  } 

(* 

@00  DEFORMED  MODEL  GEOMETRY 

*) 

model3  =  { 

Polygon [ppnts223  ], Polygon [ppnts 224  ], Polygon [ppnts225  ], 
Polygon [ppnts226  ], Polygon [ppnts227  ], Polygon [ppnts 228  ], 


} 

Show[  Graphics3D[model3, 

Boxed  ->  False, 

Axes  ->  True, 

Ticks  ->  None, 

Lighting  ->  True, 

RenderAll  ->  True, 

Viewpoint  ->  {1.3, -2. 4, 2.0}, 

ViewVertical  ->  {0,0,1}, 

DefaultFont  ->  {"Times-Italic" ,6}, 
PlotLabel  -> 

FontForm ["Buckling  Mode", 

{"Helvetica-Bold", 9}] , 

AxesLabel  ->  {"X" , "Y" , "Z"}  ]  ] 

Clear All  [tpntsl  ,tpnts2  ,tpnts3  ,tpnts4  ,tpnts5  ] 
ClearAll[tpnts6  ,tpnts7  ,tpnts8  ,tpnts9  ,tpntsl0  ] 
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Appendix  G 

Sample  User-Defined  Subroutine  to  Output  Model  Data  for  Graphical  Display 


121 


INTENTIONALLY  LEFT  BLANK. 


122 


SUBROUTINE  USRGRF  (NELEM, COORDS, U,NSETI,NELNUM,NUMNOD, 

1  NWTRN , NDFRM , GMFAC , SCALE , NPATH , NTG , 

2  NFLAG,NDEV,NTOTF,MDXl,MDX2,MDX3  ) 

C 

C  :^e*5|c***:ic;te3|cjfc5t::(e***:fc3|£jK*************** *********** 

C  **  ** 

C  **  USER-DEFINED  GRAPHICS  OUTPUT  ROUTINE  ** 

C  **  ** 

Q  jte**5jc*5!e*5|cj|c  +  5|e:^c**  ************************  ****** 

C 

IMPLICIT  DOUBLE  PRECISION (A-H,0-Z) 

C 

DIMENSION  NELEM (MDX3, 12) , COORDS (MDX2, 4) ,U(MDX1) 
DIMENSION  NSETI(MDXl) 

C 

C 

C  NODE  CONVENTION 

C 

C  Y.V 

C  Z,W  / 


c 

1 

/ 

SURFACE 

c 

1  8— 

- 7 

FACE 

NODES 

c 

1  /I 

/I 

FI: 

(1,2, 6, 5) 

c 

1  /  1 

/  1 

F2: 

(2, 3, 7, 6) 

c 

1/  1 

/  1 

F3: 

(3, 4, 8, 7) 

c 

5 - 

- 6  1 

F4: 

(1,5, 8, 4) 

c 

1  1 

1  1 

F5: 

(1,4, 3, 2) 

c 

1  4— 

- 1 - 3 

F6: 

(5, 6, 7, 8) 

c 

1  / 

1  / 

c 

1  / 

1  / 

c 

1/ 

1/ 

c 

1 - 

- 2 - 

—  X,U 

c 

c 

c 

«< 

PROGRAM  STATEMENTS  TO  OUTPUT  GRAPHICAL  »> 

c 

«< 

INFORMATION  IN 

\  USER-DEFINED  FORMAT  »> 

c 

<« 

ELEMENT  FAILURE  MEASURES: 

»> 

c 

«< 

»> 

c 

«< 

PRl:  PERCENT  OF 

PLIES  WITH 

FIBER  FAILURE  »> 

c 

«< 

PR2:  PERCENT  OF 

PLIES  WITH  MATRIX  FAILURE  »> 

C  «<  PR3:  PERCENT  MEASURE  OF  BUCKLING  FAILURE  »> 

C  «<  PR4:  PERCENT  MEASURE  OF  TOTAL  FAILURE  »> 

C 

C  «<  NDEV  =  15  FOR  UNIT  ASSOCIATED  WITH  FILE  restran.prp  »> 
C  «<  CONTAINING  GRAPHICS  GENERATED  DURING  PREPASS  »> 

C  «<  =  24  FOR  UNIT  ASSOCIATED  WITH  FILE  restran.grf  »> 

C  «<  FOR  GRAPHICS  GENERATED  DURING  ANALYSIS  »> 

C 

C  «<  IF  TOTAL  FAILURE  FLAG  IS  SET  (NTOTF  =  1) ,  THE  »> 

C  «<  DEFLECTION  OR  MODE  SHAPE  VECTOR  WILL  »> 

C  «<  GENERALLY  BE  NULL  »> 

C 

C  m®  SAMPLE  ROUTINE  FOR  PROCESSING  ELEMENT  GEOMETRY  ®@@ 

C  @@@  AND  NODE  DISPLACEMENTS  FOR  DEFORMATION  PLOT  @@@ 
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C 

C  @@@  DETERMINE  COORDINATES  OF  ELEMENT  NODES  FOR  A  @@@ 
C  @@@  PLOT  OF  UNDEFORMED  GEOMETRY  @@@ 
C 


DO  I  =  1,  NELNUM 
C 

NELID  =  NELEM(I,1) 

DO  J  =  1,  9 

M  =  NELEM(I,J+1) 

C 

C  @@@  SEARCH  COORDINATE  ARRAY  FOR  ELEMENT  NODE  @@@ 
C 

DO  K  =  1,  NUMNOD 

IF  (  M  .EQ.  INT(COORDS(K,l))  )  THEN 
N  =  K 
GO  TO  100 
END  IF 
END  DO 

100  CONTINUE 

C 

C  @@@  COORDINATES  OF  NODE  M  IN  ELEMENT  I  ARE:  @@@ 

C 

CX  =  COORDS (N, 2) 

CY  =  COORDS (N, 3) 

CZ  =  COORDS (N, 4) 

C 

END  DO 
C 

END  DO 


@@@  EXTRACT  DISPLACEMENTS  AT  NODES  TO  ADD  TO  NODE  @@@ 

C  @@@  COORDINATES  FOR  A  PLOT  OF  DEFORMED  GEOMETRY  @@@ 

C 

DO  I  =  1,  NELNUM 
C 

DO  K  =  1 ,  8 

N  =  NELEM(I,K+1) 

DO  IC  =  1,  NUMNOD 

IF  (  N  .EQ.  INT(C00RDS(IC,1))  )  THEN 
NL  =  IC 
GO  TO  200 
END  IF 
END  DO 

200  CONTINUE 

C 

C  @@@  COORDINATES  FOR  THE  K+1  NODE  IN  ELEMENT  I  DETERMINED  @@@ 
C 

CX  =  COORDS (NL, 2) 

CY  =  COORDS (NL, 3) 

CZ  =  COORDS (NL, 4) 

C 

C  @@@  DISPLACEMENTS  AT  THE  CURRENT  NODE  ARE  OBTAINED  AS 
C 


DX  =  U(3*N-2) 
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DY  =  U(3*N-1) 

DZ  =  U(3+N  ) 

@@@  THE  LOCATION  OF  THE  DISPLACED  NODE  IS  GIVEN  BY  @0® 

DNX  =  CX  +  DX 
DNY  =  CY  +  DY 
DNZ  =  CZ  +  DZ 

@@@  THESE  COORDINATES  ARE  THEN  USED  TO  DEVELOP  THE  @@@ 
@@@  THE  GEOMETRY  IN  A  USER-DEFINED  FORMAT  TO  PLOT  @0® 
C  @00  THE  DEFORMED  GEOMETRY  OF  THE  MODEL  000 

C 

END  DO 
C 

END  DO 
C 

RETURN 

END 
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INTENTIONALLY  LEFT  BLANK. 
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NO.  OF 
COPIES 

2 


1 

1 

1 

1 

1 

1 


ORGANIZATION 


NO.  OF 
COPIES 


DEFENSE  TECHNICAL  1 

INFORMATION  CENTER 
DTIC  OCA 

8725  JOHN  J  KINGMAN  RD 
STE  0944 

FT  BELVOIR  VA  22060-6218 

3 

HQDA 

DAMO  FDT 

400  ARMY  PENTAGON 

WASHINGTON  DC  20310-0460 

OSD  3 

OUSD(A&T)/ODDR&E(R) 

DRRJTREW 

3800  DEFENSE  PENTAGON 

WASHINGTON  DC  20301-3800 

COMMANDING  GENERAL 
US  ARMY  MATERIEL  CMD 
AMCRDATF 

500 1  EISENHOWER  AVE  2 

ALEXANDRIA  VA  22333-0001 

INST  FOR  ADVNCD  TCHNLGY 
THE  UNTV  OF  TEXAS  AT  AUSTIN 
3925  W  BRAKER  LN  STE  400 
AUSTIN  TX  78759-5316 

DARPA 

SPECIAL  PROJECTS  OFFICE 

JCARLINI 

3701  N  FAIRFAX  DR 

ARLINGTON  VA  22203-1714 

US  MILITARY  ACADEMY 

MATH  SCI  CTR  EXCELLENCE 

MADNMATH 

MAJ  HUBER 

THAYER  HALL 

WEST  POINT  NY  10996-1786 

DIRECTOR 

US  ARMY  RESEARCH  LAB 

AMSRLD 

DR  D  SMITH 

2800  POWDER  MILL  RD 

ADELPHI MD  20783-1 197 


ORGANIZATION 

DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  Cl  AI R 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  Cl  LL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1 197 

DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  Cl  1ST 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1 197 


ABERDEEN  PROVING  GROUND 
DIRUSARL 

AMSRL  Cl  LP  (BLDG  305) 
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NO.  OF 
COPIES 

1 


1 

3 

1 

1 

1 

1 

4 


ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  CP  CA 
D  SNIDER 

2800  POWDER  MILL  RD 
ADELPHI MD  20783-1 145 

DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  OP  SD  TA 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1 145 

DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  OP  SD  TL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1 145 

DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  Cl  1ST 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1 145 

DIRECTOR 
DA  OASARDA 
SARD  SO 

103  ARMY  PENTAGON 
WASHINGTON  DC  203 1 0-0103 

DPTY  ASST  SECY  FOR  R&T 
SARDTT 
THE PENTAGON 
RM  3EA79 

WASHINGTON  DC  20301-7100 

COMMANDER 

US  ARMY  MATERIEL  CMD 

AMXMIINT 

5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 

COMMANDER 
US  ARMY  ARDEC 
AMSTAARCC 
G PAYNE 
J  GEHBAUER 
C  BAULIEU 
HOPAT 

PICATINNY  ARSENAL  NJ 
07806-5000 


2  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  AE  WW 
E BAKER 
J  PEARSON 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  ARTD 
C  SPINELLI 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  ARFSE 
PICATINNY  ARSENAL  NJ 
07806-5000 

6  COMMANDER 

US  ARMY  ARDEC 
AMSTAARCCHA 
W  ANDREWS 
S  MUSALLI 
RCARR 
M  LUCIANO 
E  LOGSDEN 
TLOUZEIRO 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  ARCCHP 
JLUTZ 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  ARFSFT 
C  LIVECCHIA 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  ASF 

PICATINNY  ARSENAL  NJ 
07806-5000 
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NO.  OF 
COPIES 

1 


1 


3 


2 


1 


2 


NO.  OF 


ORGANIZATION 

COPIES 

ORGANIZATION 

COMMANDER 

1 

COMMANDER 

US  ARMY  ARDEC 

US  ARMY  ARDEC 

AMSTAARQACTC 

AMSTAARWET 

C  PATEL 

TSACHAR 

PICATINNY  ARSENAL  NJ 

BLDG  172 

07806-5000 

PICATINNY  ARSENAL  NJ 

07806-5000 

COMMANDER 

US  ARMY  ARDEC 

9 

COMMANDER 

AMSTAARM 

US  ARMY  ARDEC 

D  DEMELLA 

AMSTAARCCHB 

PICATINNY  ARSENAL  NJ 

P  DONADIA 

07806-5000 

F  DONLON 

P  VALENTI 

COMMANDER 

C  KNUTSON 

US  ARMY  ARDEC 

G  EUSTICE 

AMSTAARFSA 

S  PATEL 

AWARNASH 

GWAGNECZ 

BMACHAK 

RSAYER 

M  CHIEFA 

F CHANG 

PICATINNY  ARSENAL  NJ 

PICATINNY  ARSENAL  NJ 

07806-5000 

07806-5000 

COMMANDER 

6 

COMMAJTOER 

US  ARMY  ARDEC 

US  ARMY  ARDEC 

AMSTAARFSPG 

AMSTAARCCL 

M  SCHIKSNIS 

F  PUZYCKI 

D  CARLUCCI 

R  MCHUGH 

PICATINNY  ARSENAL  NJ 

D  CONWAY 

07806-5000 

E  JAROSZEWSKI 

R  SCHLENNER 

COMMANDER 

MCLUNE 

US  ARMY  ARDEC 

PICATINNY  ARSENAL  NJ 

AMSTAARFSP  A 

07806-5000 

P  KISATSKY 

PICATINNY  ARSENAL  NJ 

5 

PMSADARM 

07806-5000 

SFAE  GCSS  SD 

COL  B  ELLIS 

COMMANDER 

M  DEVINE 

US  ARMY  ARDEC 

W  DEMASSI 

AMSTAARCCHC 

J  PRITCHARD 

HCHANIN 

SHROWNAK 

S  CHICO 

PICATINNY  ARSENAL  NJ 

PICATINNY  ARSENAL  NJ 

07806-5000 

07806-5000 

1 

US  ARMY  ARDEC 

COMMANDER 

INTELLIGENCE  SPECIALIST 

US  ARMY  ARDEC 

AMSTAARWELF 

AMSTAARQACT 

M  GUERRIERE 

D  RIGOGLIOSO 

PICATINNY  ARSENAL  NJ 

PICATINNY  ARSENAL  NJ 

07806-5000 

07806-5000 
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NO.  OF 
COPIES 

2 


12 


1 


1 


1 


NO.  OF 

ORGANIZATION  COPIES 

PEO  FIELD  ARTILLERY  SYS  3 

SFAE  FAS  PM 

H  GOLDMAN 

T  MCWILLIAMS 

PICATINNY  ARSENAL  NJ 

07806-5000 

PMTMAS 
SFAE  GSSC  TMA 

R  MORRIS  1 

CKIMKER 

D  GUZIEWICZ 

E  KOPACZ 

RROESER 

R DARCY 

R  KOWALSKI 

RMCDANOLDS  1 

L  D  ULISSE 
C  ROLLER 
J  MCGREEN 
B  PATTER 

PICATINNY  ARSENAL  NJ 
07806-5000 

1 

COMMANDER 
US  ARMYARDEC 
AMSTAARWEA 
J  BRESCIA 

PICATINNY  ARSENAL  NJ 
07806-5000 

COMMANDER  2 

US  ARMYARDEC 

PRODUCTION  BASE 

MODERN  ACTY 

AMSMC  PBM  K 

PICATINNY  ARSENAL  NJ 

07806-5000 

COMMANDER 
US  ARMYTACOM 

PM  ABRAMS  1 

SFAE  ASM  AB 

6501  ELEVEN  MILE  RD 

WARREN  MI  48397-5000 

COMMANDER 
US  ARMYTACOM 
AMSTA  SF 

WARREN  MI  48397-5000 


ORGANIZATION 

COMMANDER 

US  ARMY  TACOM 

PM  TACTICAL  VEHICLES 

SFAETVL 

SFAE  TVM 

SFAETVH 

6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

COMMANDER 
US  ARMY  TACOM 
PM  BFVS 
SFAE  ASM  BV 
6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

COMMANDER 

US  ARMY  TACOM 

PMAFAS 

SFAE  ASM  AF 

6501  ELEVEN  MILE  RD 

WARREN  MI  48397-5000 

COMMANDER 
US  ARMY  TACOM 
PM  RDT&E 
SFAE  GCSS  W  AB 
J  GODELL 

6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

COMMANDER 
US  ARMY  TACOM 
PM  SURV  SYS 
SFAE  ASM  SS 
TDEAN 

SFAE  GCSS  W  GSI M 
D COCHRAN 
6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

USARMYCERL 

RLAMPO 

2902  NEWMARK  DR 
CHAMPAIGN  IL  61822 


130 


NO.  OF 
COPIES 
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1 


1 


1 


2 

3 


ORGAFnZATION 


NO.  OF 

COPIES  ORGANIZATION 


COMMANDER 
USARMYTACOM 
PM  SURVIVABLE  SYSTEMS 
SFAE  GCSS  W  GSI H 
MRYZYI 

6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

COMMANDER 

USARMYTACOM 

PMBFV 

SFAE  GCSS  W  BV 
S  DAVIS 

6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

COMMANDER 
USARMYTACOM 
CHIEF  ABRAMS  TESTING 
SFAE  GCSS  W  AB  QT 
TKRASKIEWICZ 
6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

COMMANDER 
WATERVLEET  ARSENAL 
SMCWV  QAE  Q 
B  VANINA 
BLDG  44 

WATERVLIET  NY  12189-4050 

TSM  ABRAMS 
ATZKTS 
SJABURG 
W  MEINSHAUSEN 
FT  KNOX  KY  40121 

ARMOR  SCHOOL 

ATZKTD 

RBAUEN 

JBERG 

APOMEY 

FT  KNOX  KY  40121 


14  COMMANDER 

USARMYTACOM 
AMSTATRR 
R  MCCLELLAND 
D  THOMAS 
J BENNETT 
D  HANSEN 
AMSTA  JSK 
S  GOODMAN 
J FLORENCE 
KIYER 

D  TEMPLETON 
A  SCHUMACHER 
AMSTA  TRD 
D  OSTBERG 
L  HINOJOSA 
BRAJU 
AMSTA  CS  SF 
H  HUTCHINSON 
F  <;rFrwAT?7 
WARREN  MI  48397-5000 

14  BENET  LABORATORIES 
AMSTA  ARCCB 
RFISCELLA 
MSOJA 
EKATHE 
M  SCAVULO 
G  SPENCER 
P  WHEELER 
SKRUPSKI 
JVASILAKIS 
G  FRIAR 
RHASENBEIN 
AMSTA  CCB  R 
S  SOPOK 
EHYLAND 
D  CRAYON 
R  DILLON 

WATERVLIET  NY  12189-4050 

2  HQ  IOC  TANK 

AMMUNITION  TEAM 
AMSIO  SMT 
R  CRAWFORD 
W  HARRIS 

ROCK  ISLAND  IL  61299-6000 
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NO.  OF 
COPIES 

2 

1 


2 


1 

1 

1 

2 


ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


COMMANDER 
US  ARMYAMCOM 
AVIATION  APPLIED  TECH  DIR 
J  SCHUCK 

FT  EUSTIS  VA  23604-5577 


2  OFC  OF  NAVAL  RESEARCH 
D  SIEGEL  CODE  351 
J KELLY 

800  N  QUINCY  ST 
ARLINGTON  VA  22217-5660 


DIRECTOR 
US  ARMYAMCOM 
SFAE  AV  RAM  TV 
D  CALDWELL 
BLDG  5300 

REDSTONE  ARSENAL  AL 
35898 


NAVAL  SURFACE  WARFARE  CTR 
DAHLGREN  DIV  CODE  G06 
DAHLGREN  VA  22448 

NAVAL  SURFACE  WARFARE  CTR 
TECH  LIBRARY  CODE  323 
17320  DAHLGREN  RD 
DAHLGREN  VA  22448 


US  ARMY  CORPS  OF  ENGINEERS 

CERDC 

TLIU 

CEWET 

TTAN 

20  MASSAVENW 
WASHINGTON  DC  203 14 

US  ARMY  COLD  REGIONS 
RSCH  &  ENGRNG  LAB 
PDUTTA 
72  LYME  RD 
HANOVER  NH  03755 

SYSTEM  MANAGER  ABRAMS 

ATZKTS 

LTCJHNUNN 

BLDG  1002  RM  110 

FT  KNOX  KY  40121 


1  NAVAL  SURFACE  WARFARE  CTR 

CRANE  DIVISION 
M  JOHNSON  CODE  20H4 
LOUISVILLE  KY  40214-5245 

8  DIRECTOR 

US  ARMY  NATIONAL  GROUND 
INTELLIGENCE  CTR 
D  LEITER 
M  HOLTUS 
M  WOLFE 
S  MINGLEDORF 
J  GASTON 

W  GSTATTENBAUER 
R  WARNER 
J  CRIDER 

220  SEVENTH  ST  NE 
CHARLOTTESVILLE  VA  22091 


USA  SBCCOM  PM  SOLDIER  SPT 
AMSSB  PM  RSS  A 
J  CONNORS 
KANSAS  ST 

NATICK  MA  01760-5057 


2  NAVAL  SURFACE  WARFARE  CTR 

U  SORATHIA 
C  WILLIAMS  CD  6551 
9500  MACARTHUR  BLVD 
WEST  BETHESDA  MD  20817 


USA  SBCCOM 
MATERIAL  SCIENCE  TEAM 
AMSSB  RSS 
J  HERBERT 
MSENNETT 
KANSAS  ST 

NATICK  MA  01760-5057 


2  COMMANDER 

NAVAL  SURFACE  WARFARE  CTR 
CARDEROCK  DIVISION 
R  PETERSON  CODE  2020 
M  CRITCHFIELD  CODE  1730 
BETHESDA  MD  20084 
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COPIES  ORGANIZATION 

8  US  ARMY  SBCCOM 
SOLDIER  SYSTEMS  CENTER 
BALLISTICS  TEAM 
JWARD 

WZUKAS 
P  CUNNIFF 
JSONG 

MARINE  CORPS  TEAM 

J  MACKIEWICZ 

BUS  AREA  ADVOCACY  TEAM 

W  HASKELL 

AMSSB  RCP  SS 

WNYKVIST 

S  BEAUDOIN 

KANSAS  ST 

NATICK  MA  01760-5019 

9  US  ARMY  RESEARCH  OFC 
A  CROWSON 

J CHANDRA 
H EVERETT 
J  PRATER 
R  SINGLETON 
G  ANDERSON 
D STEPP 
D  KISEROW 
J  CHANG 
PO  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
27709-2211 

8  NAVAL  SURFACE  WARFARE  CTR 

J  FRANCIS  CODE  G30 
D  WILSON  CODE  G32 
R  D  COOPER  CODE  G32 
JFRAYSSECODEG33 
E  ROWE  CODE  G33 
T  DURAN  CODE  G33 
L  DE  SIMONE  CODE  G33 
R  HUBBARD  CODE  G33 
DAHLGRENVA  22448 

1  NAVAL  SEA  SYSTEMS  CMD 

D  LIESE 

253 1  JEFFERSON  DAVIS  HWY 
ARLINGTON  VA  22242-5160 

1  NAVAL  SURFACE  WARFARE  CTR 

M  LACY  CODE  B02 
17320  DAHLGREN  RD 
DAHLGRENVA  22448 


NO.  OF 

COPIES  ORGANIZATION 

2  NAVAL  SURFACE  WARFARE  CTR 

CARDEROCK  DIVISION 
R  CRANE  CODE  2802 
C  WILLIAMS  CODE  6553 
3A  LEGGETT  CIR 
BETHESDA  MD  20054-5000 

1  EXPEDITIONARY  WARFARE 

DIVN85 
FSHOUP 

2000  NAVY  PENTAGON 
WASHINGTON  DC  20350-2000 

1  AFRLMLBC 

2941PSTRM  136 

WRIGHT  PATTERSON  AFB  OH 

45433-7750 

1  AFRLMLSS 
R THOMSON 

2179  12THSTRM122 
WRIGHT  PATTERSON  AFB  OH 
45433-7718 

2  AFRL 

F  ABRAMS 

J  BROWN 

BLDG  653 

2977  P  ST  STE  6 

WRIGHT  PATTERSON  AFB  OH 

45433-7739 

1  WATERWAYS  EXPERIMENT 

D SCOTT 

3909  HALLS  FERRY  RD  SC  C 
VICKSBURG  MS  39180 

5  DIRECTOR 

LLNL 

R  CHRISTENSEN 
S  DETERESA 
FMAGNESS 
M  FINGER  MS  313 
M  MURPHY  L  282 
PO  BOX  808 
LIVERMORE  CA  94550 

1  AFRL  MLS  OL 

L  COULTER 
7278  4TH  ST 
BLDG  100  BAY  D 
HILL  AFB  UT  84056-5205 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


I  OSD 

JOINT  CCD  TEST  FORCE 
OSD  JCCD 
R  WILLIAMS 
3909  HALLS  FERRY  RD 
VICKSBURG  MS  29180-6199 

1  DEFENSE  NUCLEAR  AGENCY 
INNOVATIVE  CONCEPTS  DIV 
6801  TELEGRAPH  RD 
ALEXANDRIA  VA  22310-3398 

3  DARPA 

MVANFOSSEN 

SWAX 

L  CHRISTODOULOU 
3701  N  FAIRFAX  DR 
ARLINGTON  VA  22203-1714 

2  SERDP  PROGRAM  OFC 
PMP2 

C  PELLERIN 
B  SMITH 

901  N  STUART  ST  STE  303 
ARLINGTON  VA  22203 

1  FAA 

MIL  HDBK  17  CHAIR 
L ILCEWICZ 
1601  LIND  AVE  SW 
ANM  115N 
RESTONVA  98055 

1  US  DEPT  OF  ENERGY 

OFC  OF  ENVIRONMENTAL 
MANAGEMENT 
P  RITZCOVAN 
19901  GERMANTOWN  RD 
GERMANTOWN  MD  20874-1928 

1  DIRECTOR 

LLNL 

F  ADDESSIO  MS  B216 

PO  BOX  1633 

LOS  ALAMOS  NM  87545 

1  OAK  RIDGE  NATIONAL 

LABORATORY 
RM  DAVIS 
PO  BOX  2008 

OAK  RIDGE  TN  37831-6195 


3  DIRECTOR 

SANDIA  NATIONAL  LABS 

APPLIED  MECHANICS  DEPT 

MS  9042 

J  HANDROCK 

YRKAN 

J  LAUFFER 

PO  BOX  969 

LIVERMORE  CA  94551-0969 

1  OAK  RIDGE  NATIONAL 

LABORATORY 
C  EBERLE  MS  8048 
PO  BOX  2008 
OAK  RIDGE  TN  37831 

1  OAK  RIDGE  NATIONAL 

LABORATORY 
CD  WARREN  MS  8039 
PO  BOX  2008 
OAK  RIDGE  TN  37831 

5  NIST 

JDUNKERS 

M  VANLANDINGHAM  MS  8621 
J  CHIN  MS  8621 
J  MARTIN  MS  8621 
DDUTHINHMS  8611 
100  BUREAU  DR 
GAITHERSBURG  MD  20899 

1  HYDROGEOLOGIC  INC 

SERDP  ESTCPSPTOFC 
S  WALSH 

1155  HERNDON  PKWY  STE  900 
HERNDON  VA  20170 

3  NASA  LANGLEY  RSCH  CTR 

AMSRLVS 
W  ELBER  MS  266 
F  BARTLETT  JR  MS  266 
G  FARLEY  MS  266 
HAMPTON  VA  23681-0001 

1  NASA  LANGLEY  RSCH  CTR 

T  GATES  MS  188E 
HAMPTON  VA  23661-3400 

1  FHWA 

EMUNLEY 

6300  GEORGETOWN  PIKE 
MCLEAN  VA  22101 
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1 

1 

1 


1 

2 

1 

1 

1 


ORGANIZATION 
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CYTEC  FIBERITE 

R DUNNE 

DKOHLI 

RMAYHEW 

1300  REVOLUTION  ST 

HAVRE  DE  GRACE  MD  21078 

USDOT  FEDERAL  RAILRD 
M  FATEH  RDV  31 
WASHINGTON  DC  20590 

MARINE  CORPS 
INTLLGNC  ACTVTY 
D  KOSITZKE 

3300  RUSSELL  RD  STE  250 
QUANTICO  VA  22134-501 1 

DIRECTOR 

NATIONAL  GRND  INTLLGNC  CTR 

lANGTMT 

220  SEVENTH  ST  NE 

CHARLOTTESVILLE  VA 

22902-5396 

SIOUX  MFG 
BKRIEL 
PO  BOX  400 
FT  TOTTEN  ND  58335 

3TEX  CORPORATION 
A  BOGDANOVICH 
J  SINGLETARY 
109MACKENANDR 
CARY  NC  275 11  . 

3M  CORPORATION 
J  SKILDUM 

3M  CENTER  BLDG  60  IN  01 
ST  PAUL  MN  55144-1000 

DIRECTOR 

DEFENSE  INTLLGNC  AGNCY 
TA5 

KCRELLING 
WASHINGTON  DC  20310 

ADVANCED  GLASS  FIBER  YARNS 
T  COLLINS 

281  SPRING  RUN  LANE  STE  A 
DOWNINGTON  PA  19335 


1  COMPOSITE  MATERIALS  INC 

D  SHORTT 
19105  63  AVENE 
PO  BOX  25 

ARLINGTON  WA  98223 

1  IPS  GLASS 

L  CARTER 
PO  BOX  260 
SLATER  RD 
SLATER  SC  29683 

1  COMPOSITE  MATERIALS  INC 

R  HOLLAND 
11  lEWELCT 
ORINDA  CA  94563 

1  COMPOSITE  MATERIALS  INC 
C  RILEY 

14530  S  ANSON  AVE 
SANTA  FE  SPRINGS  CA  90670 

2  SIMULA 
JCOLTMAN 
RHUYETT 
10016  S  5 1ST  ST 
PHOENIX  AZ  85044 

2  PROTECTION  MATERIALS  INC 

M  MILLER 
F  CRILLEY 
14000  NW  58  CT 
MIAMI  LAKES  FL  33014 

2  FOSTER  MILLER 

M  ROYLANCE 
WZUKAS 
195  BEAR  HILL  RD 
WALTHAM  MA  02354-1196 

1  ROM  DEVELOPMENT  CORP 

ROMEARA 
136  SWINEBURNE  ROW 
BRICK  MARKET  PLACE 
NEWPORT  RI 02840 
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2  TEXTRON  SYSTEMS 

T  FOLTZ 
M  TREASURE 
1449  MIDDLESEX  ST 
LOWELL  MA  01851 

1  O  GARA  HESS  &  EISENHARDT 
M  GILLESPIE 
9113LESAINTDR 
FAIRFIELD  OH  45014 

2  MILLIKEN  RSCH  CORE 
HKUHN 

M  MACLEOD 
PO  BOX  1926 
SPARTANBURG  SC  29303 

1  CONNEAUGHT  INDUSTRIES  INC 

J  SANTOS 
PO  BOX  1425 
COVENTRY  RI 02816 

1  BATTELLE  NATICK  OPNS 

BHALPIN 

209  W  CENTRAL  ST  STE  302 
NATICK  MA  01760 

1  ARMTEC  DEFENSE  PRODUCTS 

SDYER 
85  901  AVE  53 
PO  BOX  848 
COACHELLA  CA  92236 

1  NATIONAL  COMPOSITE  CENTER 
T  CORDELL 

2000  COMPOSITE  DR 
KETTERING  OH  45420 

3  PACIFIC  NORTHWEST  LAB 
M  SMITH 

G  VAN  ARSDALE 
R  SHIPPELL 
PO  BOX  999 
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